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Prof.  Tomlinson  Fort  of  the  University  of  Georgia  and  Dr. 
Alfred  Hind  were  engaged  on  this  project  from  March  1,  1952  to 
August  31,  1952.  Prof.  Fort  was  director  of  the  research  and 
Dr.  Hina  was  his  assistant.  Prof.  Fort  alone  has  been  engag- 
ed on  this  project  from  Sept.  1,  1952  to  Feb.  28,  1953. 

Two  papers  have  been  made  ready  fcr  publication.  The 
first  is  entitled} 

*Xhs  Lg&aafl-  Vibrating  Mat  sad  Baanltlas  flauadaixj&lMfl. 

Problems  for  a Partial  Difference  Equation  of  the 

Second  Order. ' 

This  paper  is  under  the  authorphip  of  Tomlinson  Fort. 

Two  copies  were  sent  to  the  Office  of  Ordnance  Research  at  Duke 
University  early  in  December  1952,  and  permission  asked  to 
submit  the  paper  for  publication  in  Proceedings  of  London 
Mathematical  Society. 

The  second  paper  is  under  the  authorship  of  Alfred  Hind. 

II  boars  »i.ao..tI»ily  the  following  title: 

'Convergence  of  approximate  solutions  of  the  vibrating 
membrane  equation  to  its  eaact  solutions- • 

Two  copies  have  been  sent  to  the  Office  cf  Ordnance  Re- 
search at  Duke  University  and  permission  asked  to  submit  the 
paper  to  the  Proceedings  of  the  American  Mathematical  Society. 

Brief  abstracts  only  will  be  given  of  these  two  papers, 
inasmuch  as  the  complete  papers  are  in  the  hands  of  the  Office 
of  Ordnance  Research  at  Duke  University.  These  abstracts 
will  be'  followed  by  an  equally  brief  resume  of  work  completed 
but  not  yet  readied  for  publication.  It  is  believed  that 
there  Is  ample  material  for  at  least  one  additional  paper. 
Certain  other  scattered  results  have  teen  obtained  and  many 
leads  followed  which  have  not  been  productive.  This  should 
be  helpful  in  further  research  in  this  field. 

Abstract  of  Paper  Number  One 

Let  there  be  given  a rectangular  net  composed  of  elastic 
cords  and  loaded  at  the  points  (i,J)  with  particles  of  mass 
mjj  moving  with  small  vibrations  in  lines  perpendicular  to  the 
plane  of  the  net  when  at  rest.  An  attempt  to  determine  the 
displacement  of  each  particle  at  time  t , leads  to  a set  of 
linear  differential  equations  of  the  second  order.  An  attempt 
to  solve  these  leads  one  to  a consideration  of  the  difference 
equation 


(1)  £b(i-lf  J^yU-l,  ^ £k(if  j-DAjyd,  j-lj 

+ It  y (i»  i)  = o 

subject  to  appropriate  boundary  conditions.  These  are  taken 
to  be 

(2)  y(o,j)  = y(5i+i.J)  = = y(i,n+i)  = o 

1 • 0 • f in | J - i,  £ f • # • f n* 


This  Is  tantamount  to  assuming  the  displacement  of  boundary  par- 
ticles to  bo  zero.  The  consideration  of  the  above  as  a char- 
acteristic-value problem  is  the  purpose  of  the  present  papor. 

The  characteristic  determinant  is  the  determinant  of  the  co- 
efficients of  the  linear  algebraic  equations  obtained  by  substi- 
tuting the  points  (i,J)  and  boundary  conditions  in  equation  (1). 
It  is  proved  that  the  characteristic  determinant  is  symmetrical 
about  its  principal  diagonal,  that  all  its  roots  are  real  and 
positive,  although  some  may  be  multiple,  as  has  boen  shown  by 
examples.  The  work  is  carried  forward  by  a kinematical  method. 
Among  other  things  it  is  proved  that  there  are  always  at  least 
m+n-1  distinct  roots  of  the  characteristic  equation.  These 
roots  are  denoted  of  ZQ,  •••»  ^nH-n-1*  Wh&n  <>  = a certain 
characteristic  function  which  is  aSnoted  by  V(i,J)  has  at  least 
^ nodal  lines  on  tho  rectangular  lattice. 

Abstract  of  naner  no.  2 

• • iu>  Friedrichs  and  Lewy  have  shown  that  under  appropri- 
c:-.  conultions  a solution  of 

(3)  i,  - 5 


A*  uf x, y-Ay,  t 



with  a mesh  ratio  of  rx  «■  r,  **  1 will  converge  to  a solution  of 
the  corresponding  differential  equation. 

Lautert  and  Obrien  proved  the  convergence  of  solutions  of 
the  equation 

A*u(x,  t-At)  _ u(x-Ax,  t)l 

(At)a  " ( IS?  J 

ba  > 0. 

to  the  corresponding  solutions  of  the  analagous  differential 
equation. 

The  present  paper  treats  equation  (3)  with  £ replaced  by 
b*>  0 . The  method  is  much  moro  oompact  than  that  employed  by 
Obrien  and  Lautert  in  their  one  dimensional  problem  and  permits 
of  immediate  extension  to  the  wave  equation  in  n-dimensions . 

It  is  proved  that  for  any  mesh  ratio  rx  = ra  = r r.  solution  of 
the  difference  equation  approaches  the  corresponding  solution 
of  the  differential  equation. 

Third  tonic  of  research  nqt  vet  submitted  for  publication. 

This  work  has  been  concerned  with  the  equation 

A|y(i-1,  1)  = Aa  y(i,  J-l) 


subjeot  to  the  condition  that  firstly,  y be  zero  on  the  boun- 
dary of  a square  lying  in  the  first  quadrant  and, secondly,  that 
y be  zero  along  the  boundary  of  a neoquadrant  of  a circle. 

Tho  method  of  attack  has  been  to  assume  a solution  in  the  form 
1(1)  J(j).  It  is  then  proved  that  any  solution  of  the  problem 
over  the  square  is  a linear  combination  of  products  of  solutions 


of  the  one  dimensional  Sturm-Lionville  problem.  It  follows 
that  all  characteristic  values  are  positive.  It  is  proved 
that  no  characteristic  value  i3  of  multiplicity  greater  than 
n . 


It  ia  prove^th&kg^glutions  exist  for  the  neo-quadrantal 
problem.  However,  ^no  soliition  car;  be  written  in  the  form 
I( 1 ) J(J)  where  I and  J are  Sturm-Lionvllie  functions  in 
one -dimension. 
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I.  Homogeneous  Nucleatlon  of  Supercooled  Vat «r  Drops 


Abstract  — The  experimentally  observed  crystal lltat Ion  of 
supercooled  water  near  -hOC  is  examined  In  terms  of  the  theory  of 
homogeneous  nucleatlon.  The  thermodynamic  and  molecular-kinetic 
nature  of  the  nucleatlon  process  Is  outlined  In  order  to  show  why 
supercooling  In  natural  clouds  can  occur  so  frequently.  Past  ef- 
forts to  explain  the  -hoc  transition  are  examined  critically  and 
are  found  to  contain  a number  of  significant  errors.  Because  the 
theoretical  nucleatlon  rates  are  extremely  sensitive  to  the  numeri- 
cal value  of  the  specific  surface  free  energy  of  a water-ice  Inter- 
face, particular  attention  lo  devoted  to  the  refinement  of  previ- 
ous estimates  of  this  parameter.  It  Is  shown  that  both  Krastanow's 
and  Mason's  estimates  were  Inaccurate,  and  that  In  the  latter's 
approach,  neglect  of  the  distortion  energy  of  the  surface  layer  of 
Ice  led  to  a marked  underestimate  of  the  nucleatlon  efficiency 
which  was  concealed  by  the  effects  of  several  counteracting  errors. 
Difficulties  lying  In  the  way  of  a direct  calculation  of  the  dis- 
tortion energy  for  Ice  are  examined  and  found  to  be  very  serious. 

A crude  correction  for  distortion  effects  leads  to  a theoretically 
predicted  temperature  of  -26C  for  the  threshold  of  spontaneous  nu- 
cleatlon of  drops  of  cloud -particle  size.  It  Is  concluded  that 
although  this  result  lies  well  above  the  experimentally  observed 
range  of  transition  temperatures,  It  Is  close  enough  to  that  range, 
considering  the  Inherent  difficulty  of  assessing  the  effect  of  dis- 
tortion, to  strengthen  the  belief  that  the  -hOC  transition  Is  due 
to  homogeneous  nucleatlon.  The  implications  of  this  conclusion 
for  the  theory  of  the  aircraft  icing  process  are  pointed  out 
briefly. 

1.  Introduction 

Although  a number  of  years  have  elapsed  since  the  first  clear  recognition 
(Cwllong,  19**5;  Schaefer,  19h6)  that  there  exists  a temperature  near  -hoc  at 
which  Ice  crystals  seem  to  form  abundantly  and  spontaneously  In  a cloud  of 
supercooled  water  drops,  the  exact  physical  nature  of  the  transition  has  not 
yet  been  ascertained.  Its  Importance  to  the  subjects  of  cloud  physics  and 
aircraft  Icing  makes  It  desirable  to  investigate  all  possible  mechanisms  which 
might  account  for  the  phenomenon.  The  principal  division  of  present  opinion 
seems  to  hinge  upon  the  question  of  whether  the  observed  formation  of  Ice 
crystals  near  -hoc  occurs  as  a result  of  direct  sublimation  from  the  vapor 
phase  (Cwllong,  19h7;  Bradley,  1951;  Schaefer,  1952)  or  as  a result  of  spon- 
taneous freezing  of  supercooled  water  drops  (Fisher,  Hollomon,  and  Turnbull, 
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19^9;  Lafargue,  1950;  Mason,  1952). 

The  sublloatlonal  hypothesis  encounters  serious  difficulties  on  both 
theoretical  and  experimental  grounds.  First,  Krastanow  (19**l)  has  shown 
thermodynamically  that  water  vapor  at  temperatures  below  0C  will  condense  to 
the  supercooled  liquid  state  In  an  energetically  easier  manner  than  It  will 
sublime  to  the  solid  state,  since  the  llquld-vapor  Interfacial  energy  Is 
certainly  smaller  than  the  solid -vapor  lnterfaclal  energy  In  the  temperature 
range  of  meteorological  Interest.  Second,  It  has  been  found  experimentally 
(e.g.,  Schaefer,  1952)  that  the  Ice  phase  does  not  make  Its  appearance  In 
vapor  saturated  Just  with  respect  to  Ice,  but  rather  In  vapor  saturated  with 
respect  to  liquid  water,  vhlch  Implies  that  formation  of  Ice  crystals  in  the 
atmosphere  proceeds  by  some  Intermediate  process  Involving  supercooled  liq- 
uid drops  (assuming  no  foreign  nuclei  to  be  present).  A very  complete  and 
critical  summary  of  past  studies  of  this  problem  has  been  given  by  Mason  ard 
Ludlam  (1951)*  sc  the  background  of  the  problem  need  not  be  elaborated  here. 

When  one  next  Inquires  as  to  how  a supercooled  drop  of  water  mleht 
freeze,  two  distinct  classes  of  processes  must  be  considered.  If  the  crys- 
tallization takes  place  In  entirely  pure  water.  It  Is  said  to  depend  upon 
"homogeneous  nucleatlon".  If  the  crystallization  Instead  depends  critical- 
ly upon  the  presence  of  trace  quantities  of  some  foreign  substances.  It  Is 
said  to  Involve  '“heterogeneous  nucleatlon".  It  has  been  clearly  demonstrated 
within  the  past  few  years  that  a large  number  of  Inorganic  salts  will  pro- 
mote heterogeneous  nucleatlon  when  added  to  suitably  supercooled  clouds,  and 
much  effoit  has  been  directed  towards  gaining  an  understanding  of  the  exact 
mechanisms  Involved.  In  contrast  to  these  experiments,  the  experiments  on 
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the  DC  transition*  seen  to  Indicate  that  heterogeneous  nucleatlon  Is  not 
Involved  In  this  phenomenon;  so  It  becomes  very  desirable  to  exploit  all 
available  experimental  and  theoretical  methods  for  gaining  Insight  Into  the 
problem  of  homogeneous  nucleatlon  of  supercooled  water.  The  present  paper 
summarizes  a theoretical  study  of  this  latter  phenomenon. 

It  may  be  noted  at  the  start  that  the  theoretical  expression  describ- 
ing the  rate  of  homogeneous  nucleatlon  Is  of  such  mathematical  form  as  to  be 
extremely  sensitive  to  the  numerical  values  assigned  to  certain  parameters 
entering  Into  It.  Consequently  a major  objective  of  the  present  study  has 
been  the  careful  examination  of  all  quantities  affecting  the  nucleatlon  ra:e, 
particular  effort  being  devoted  to  refining  such  numerical  estimates  of 
these  quantities  as  have  been  made  In  the  past.  In  addition,  several  thermo- 
dynamic aspects  of  the  nucleatlon  problem  have  been  Investigated  and  shown 
to  be  In  need  of  revision.  After  all  recognized  refinements  have  been  dis- 
cussed here,  calculations  of  the  nucleatlon  rate  will  be  carried  out  on  the 
basis  of  three  Independent  and  disparate  estimates  of  the  critically  Import- 
ant surface  free  energy  for  a water-ice  Interface. 

2.  Theory  of  homogeneous  nucleatlon 

One  of  the  principal  benefits  to  be  derived  from  even  a qualitative 
study  of  the  theory  of  homogeneous  nucleatlon  of  supercooled  liquids  Is  an 
appreciation  of  tbe  fact  that  as  the  temperature  of  a pure  liquid  substance 
le  lowered  to  end  then  below  the  melting  point  of  that  substances  solid 
phase,  a certain  degree  of  supercooling  Is  not  only  possible,  but  Is  In  fact 
difficult  to  avoid.  Thus  the  well-known  and  meteorologically  important  phe- 
nomenon of  supercooling  of  cloud  drops  ought  not  be  regarded  by  meteorolo- 

*Dlfferent  Investigators  have  observed  the  transition  In  question  at  slightly 
different  temperatures  lying  within  the  range  from  about  -38c  to  about  -4 1C. 
Here,  for  brevity,  It  will  be  referred  to  simply  as  the  "-40C  transition". 


-le- 
gists as  a paradox  demanding  some  special  explanation,  sieved  In  the  light 
of  the  theory  of  homogeneous  nucleatlon,  this  phenomenon  la  eeen  to  be  thermo- 
dynamically and  klnetlcally  Inevitable  In  pure  water  drops.  That  It  would 
really  be  the  opposite  caae  of  "we 11 -behaved"  free ting  at  OC  which  would 
demand  special  explanation  if  Invariably  observed  In  the  atmoephere  does 
not  appear  to  be  widely  appreciated!  If  one  Is  to  judge  from  consents  cm  the 
phenomenon  of  supercooling  to  be  found  In  many  meteorology  texts.  Therefore 
It  may  be  In  order  here  to  discuss  briefly  the  qualitative  nature  of  this 
problem  before  turning  to  a quantitative  examination. 

The  supercooling  of  a pure  liquid  Is  quite  closely  analogous  to  the 
supersaturAtion  of  a pure  vapor.  Since  the  latter  Is  of  conceptual  Interest 
In  cloud  physics,  but  particularly  since  certain  quantitative  aspects  of  the 
latter  are  likely  to  be  better  known  to  most  readers,  the  basic  physical  na- 
ture of  homogeneous  nucleatlon  will  be  described  first  In  terms  of  processes 
occurring  In  a supersaturated  vapor. 

In  this  latter  case,  one  observes  that  a vapor  which  contains  no  inter- 
mixed foreign  partlclee  and  which  Is  not  In  contact  with  a liquid  water  sur- 
face of  any  kind  may  be  cooled  far  below  Its  nominal  dev  point  without  any 
appearance  of  water  drops.  The  fundamental  reason  for  this  possibility  of 
supersaturatlon  of  pure  vapor  is  that  the  only  way  in  which  liquid  drops  can 
be  formed  under  such  homogeneous  conditions  Is  for  a chance  succession  of 
collisions  to  build  up  embryonic  water  droplets  one  molecule  at  a time.  But 
this  mode  of  formation  of  drops  Is  Inherently  Improbable,  even  at  tempera- 
tures well  below  the  nominal  dev  point  of  the  vapor  because  (to  describe  It 
In  simplified  terms)  the  total  binding  force  exerted  by  very  small  aggregates 
of  molecules  upon  one  of  their  surface  members  is  too  low  to  overcome  the 
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disgregatire  action  of  thermal  sot  Iona  of  the  individual  members.  Onlj  bj 
cooling  the  system  far  belov  Its  nominal  dev  point  can  one  suppress  this 
thermal  dlsgre gating  effect  on  randomly  forming  embryos  sufficiently  to  per- 
mit the  embryos  to  gain  member  molecules  faster  by  bombardment  from  the  va- 
por phase  than  the  embryos  lose  them  to  the  vapor  phase  by  evaporation. 

Sven  from  as  crude  a model  as  this,  one  can  sense  that  for  any  given  tem- 
perature belov  the  nominal  dev  point  there  ought  to  exist  a critical  embryo 
radius  belov  vhlch  an  embryo  is  unstab  La  and  vlll  tend  to  evaporate,  but 
above  vhlch  it  vlll  tend  to  grow  rapidly.  Furthermore,  the  same  model  sug- 
gests the  important  relation  that  the  magnitude  of  thlB  critical  radius  must 
surely  decrease  vlth  decreasing  temperature.  Both  of  these  conclusions  can 
be  given  firm  support  and,  still  better,  quantitative  expression  by  means 
of  the  corresponding  thermodynamic  arguments. 

One  of  the  main  results  of  the  thermodynamic  argument  is  a quite  famil- 
iar relation,  Kelvin's  equation.  Hovever,  the  form  into  vhlch  this  equation 
is  almost  invariably  put  in  meteorological  applications,  and  the  correspond- 
ing vay  in  vhlch  it  is  interpreted,  tend  to  conceal  its  interesting  impli- 
cations as  an  equation  from  the  theory  of  homogeneous  nucleatlon.  Kelvin's 
equation,  as  conventionally  vritten,  relates  the  vapor  tension  er  of  a pure 
vater  drop  of  radius  r and  temperature  T to  the  vapor  tension  e of  a plane 


cFor  aggregates  containing,  say,  only  a half  dozen  molecules,  the  binding 
force  for  any  given  member  molecule  may  be  understood  in  terms  of  the  small 
total  number  of  other  molecules  available  to  attract  the  given  molecule. 

For  larger  aggregates,  the  total  number  of  members  ceases  to  be  the  signif- 
icant para-  *r,  since  it  is  veil  established  that  the  range  of  lotermole- 
cuiar  forces  is  ao  small  as  to  become  generally  negligible  over  distances 
of  the  order  of  tvo  or  three  molecular  diameters.  Instead,  the  deficiency 
in  binding  energy  for  a surface  molecule  in  an  aggregate  of,  say,  a fev  hun- 
dred molecules  depends  almost  entirely  on  a purely  geometric  effect  due  to 
the  large  surface  curvature  of  such  tiny  aggregates.  The  curving  surface  of 
the  embryo  falls  avay  so  rapidly  on  all  sides  of  a given  surface  molecule 
that  there  exists  a slight  deficit  of  near -neighbors  vhose  radially  invard 
components  of  attraction  go  to  make  up  a portion  of  the  total  binding  force 
in  a perfectly  plane  liquid  surface . 


surface  of  pure  water  of  density  ^ , specific  surface  free  energy  (surface 
tension)  (fy  > and  temperature  T,  according  to 

.^er  — ^_^L.  (1) 

e ftR'Tr 

where  R'  Is  the  gas  constant  per  gram  for  water  vapor.  It  Is  usually,  »nd 
correctly,  stated  that  Kelvin’s  equation  shows  how  much  larger  the  vapor  ten- 
sion of  a tiny  drop  Is  than  that  of  a plane  surface  of  the  same  liquid  at  the 
same  temperature.  What  Is  too  frequently  Ignored  In  meteorological  discus- 
sions, Is  that  the  same  equation  is  equally  appropriately  regarded  as  an  equa- 
tion giving  the  critical  embryo  radius  cs  a function  of  the  degree  of  super- 
saturation  and  hence  of  the  degree  of  undercooling  below  the  nominal  dewpoint. 
TO  bring  out  the  latter  interpretation,  one  may  substitute  for  ln(er/e)  from 
the  Integrated  form  of  Clapeyron’s  equation  to  get,  on  expressing  it  as  an 


equation  to  determine  r, 

r = 


a az  lo 


where  Ly  Is  the  latent  heat  of  vaporization  of  the  liquid  and  T0  Is  the  tem- 
perature (greater  than  T)  at  which  er  corresponds  to  the  saturation  vapor 
pressure,  i.e.,  TQ  Is  simply  the  nominal  dew  point  of  the  given  sample  of 
supersaturated  vapor  which  Is  In  equilibrium  with  drops  of  radius  r.  Equa- 
tion (2),  though  still  essentially  Kelvin's  equation,  tells  rather  more  than 
(1),  for  It  defines  the  radius  r that  an  embryo  must  Just  attain  to  grow 
rather  than  to  evaporate  In  the  presence  of  water  vapor  which  has  been  under* 
cooled  Isobarically  from  Its  nominal  dew  point  T0  to  Its  actual  temperature 
T.  Kelvin's  familiar  but  amazing  derivation  of  (1)  from  consideration  of  a 
capillary  column  (Humphreys,  19^0)  completely  conceals  this  physical  Inter- 
pretation. A more  straightforward  thermodynamic  derivation  brings  It  out 
clearly  since  r In  (1)  or  (2)  Is  then  seen  to  represent  the  radius  for  which 
the  net  free  energy  change  due  to  addition  of  molecules  to  an  embryo  attains 
a maximum.  For  embryos  smaller  than  this  critical  size,  the  free  energy  in- 
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crease  due  to  increase  of  surface  area  accompany ing  the  addition  of  more  mole- 


cules more  than  cancels  the  free  energy  decrease  accompanying  the  phase  change. 
Hence  in  this  r-range,  further  growth  does  not  tend  to  occur  "spontaneously" 
in  the  thermodynamic  sense,  l.e.,  the  process  does  not  involve  a net  decrease 
of  free  energy.  For  embryos  greater  than  the  critical  site,  the  r2  depend- 
ence of  the  surface  term  is  overpowered  by  the  r5  dependence  of  the  bulk 
term  and  rapid  growth  ensues  spontaneously.  Viewed  in  this  light,  the  criti- 
cal radius  given  by  (2)  is  that  for  which  the  free  energy  has  a stationary 
value. 

To  gain  a firmer  feeling  for  the  homogeneous  nucleatlon  process  as  it 
operates  in  the  case  of  supersaturated  vapor,  one  may  use  (2)  to  compute  the 
actual  size  of  these  critical  embryos  in  a particular  case.  Thus,  a sample 
of  water  vapor  which  has  a vapor  pressure  of  17  mbs  is  known  to  be  nominally 
saturated  at  15C,  If  cooled  isubarlcally  to  5C  in  the  absence  of  condensa- 
tion nuclei  and  free  water  surfaces,  it  attains  a relative  humidity  of  about 
200  per  cent,  and  (2)  reveals  that  homogeneous  nucleatlon  cannot  then  occur 
unless  chance  collisions  build  up  embryos  with  radii  of  about  17  Angstrom 
units.  Since  this  radius  is  of  the  order  of  ten  times  the  molecular  radius 
of  water,  it  follows  that  spontaneous  condensation  will  not  begin  unless  ag- 
gregates of  some  10^  water  molecules  are  built  up.  This  turns  out  to  be  so 
very  Improbable  at  the  vapor  density  and  temperature  in  question  that  the 
chance  of  even  one  such  critical  embryo  appearing  in  a volume  of  many  liters 
of  vapor  in  a time  of  many  minutes  is  negligible. 

The  preceding  discussion  of  homogeneous  nucleatlon  in  a supersaturated 
vapor  raises  the  question  of  how  one  predicts  the  average  rate  at  which  em- 
bryos of  the  critical  size  may  form  by  random  molecular  processes  at  a given 
degree  of  undercooling.  This  question  is  of  equal  significance  in  the  prob- 
lem under  discussion  in  this  psper,  namely  nucleatlon  of  a supercooled  liquid; 


so,  having  outlined  the  basic  features  of  a typical  homogeneous  nucleatlon 
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procesa,  attention  will  now  be  returned  to  this  problem  of  primary  Interest 
here.  Just  as  was  true  for  the  case  of  nucleatlon  of  vapor,  one  can  only  ex- 
pect to  find  embryos  (now  crystalline  embryos)  forming  In  a supercooled  liq- 
uid as  a result  of  random  collision  processes;  and,  again,  survival  of  these 
Is  prejudiced  by  chance  until  they  can  grow  to  such  site  that  the  surface 
molecules  In  the  embryonic  crystal  lattice  are  bound  with  energies  rather 
greater  than  the  average  thermal  vibrational  energy  corresponding  to  the  tem- 
perature of  the  liquid.  An  equation  very  similar  to  (2),  namely  (5)  below, 
specifies  the  critical  crystallite  radius  for  any  given  degree  of  supercool- 
ing, and  the  problem  Is  to  find  a theoretical  expression  fer  the  rate  of 
formation  of  these  critical  embryos  per  unit  volume  of  supercooled  liquid  as 
a function  of  the  degree  of  supercooling. 

This  problem  has  been  solved  through  the  successive  efforts  of  Becker 
and  Boring  (1935)#  Turnbull  and  Fisher  (19^9).  and  others.  An  extensive  dis- 
cussion of  the  nucleatlon  rate  problem  (exclusive  of  the  recent  and  Important 
contribution  of  Turnbull  and  Fisher)  may  be  found  In  Frenkel  (19^6);  and  a 
more  recent  summary  has  been  given  by  Bradley  (1951).  In  a liquid  supercooled 
below  Its  nominal  freezing  point  TQ  to  an  actual  temperature  T,  the  rate  of 
formation  of  nuclei  (critical  embryos)  per  unit  volume  per  unit  time  Is 
given  by  Turnbull  and  Fisher  as  r\  . c 


where  J Is  the  specified  nucleatlon  rate,  n is  the  number  of  molecules  per 
unit  volume  In  the  liquid  phase,  k Is  Boltzmann's  constant,  u is  Planck's 
constant,  A is  the  free  energy  of  activation  for  aelf-d If fusion  of  the  liq- 
uid molecules,  and  Fc  is  the  free  energy  of  formation  of  a nucleus.  The  con- 
tribution to  this  theoretical  equation  made  by  Turnbull  and  Fisher  was  the  ap- 
proximate evaluation  of  the  pre factor  of  the  exponential  in  (5)  based  on  the 
quantum  statistical  theory  of  absolute  reaction  rates  (Glasstone,  Laidler, 
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and  Eyrlng,  19^1).  It  Is  to  be  noted  that  nkT/h  1b  not  simply  the  collision 
frequency,  as  vas  once  thought  to  be  true.  Very  loosely,  (3)  nay  be  said  to 
give  J as  the  product  of  a fundamental  molecular  (or  atomic)  rate  constant, 
nkT/h,  multiplied  by  two  fractional  probability  factors  rather  similar  to 
Boltzmann  factors.  The  first,  exp(-A/kT),  cuts  down  the  fundamental  rate  by 
a fraction  measuring  the  rather  low  probability  that  any  given  liquid  molecule 
will  possess  at  any  given  Instant  sufficient  energy  to  permit  It  to  "break 
loose"  from  the  liquid  structure  and  reach  an  energy  state  In  which  It  1b 
free  to  diffuse  from  the  liquid  region  onto  and  Into  the  embryonic  lattice 
(see  Section  9 below).  The  second  factor,  exp(-Fc/kT),  further  reduces  J by 
a factor  which  measures  the  generally  very  small  probability  that  random  pro- 
cesses can  succeed  In  amassing  an  embryo  of  critical  size  In  the  face  of  the 
associated  free  energy  Increase  of  the  system.  Clearly,  homogeneous  nuclea- 
tlon  should  be  thought  of  as  an  Intrinsically  Improbable  event,  as  has  been 
stressed  above  In  pointing  out  that  the  common  occurrence  of  supercooled 
cloud  drops  Is  not  really  paradoxical  at  all. 

The  activation  energy  for  self-d If fusion,  A,  cannot  readily  be  expressed 
as  a thermodynamic  function  of  more  elementary  parameters,  so  Its  evaluation 
will  be  reserved  for  Section  9«  On  the  other  hand,  the  free  energy  Increase, 
Fc,  associated  with  the  formation  of  an  embryo  of  critical  radius  rc  (or  of 
characteristic  length  rc  If  one  visualizes  formation  of  some  other  geometric 
shape  than  a sphere)  Is  found  from  thermodynamic  considerations  (see,  for 
example,  Frenkel,  19^6)  to  be 


where  Is  the  specific  surface  free  energy  of  the  solid-liquid  Interface 
and  g Is  a geometric  factor  such  that  grc^  Is  the  total  surface  area  of  the 
critical  embryo.  An  expression  of  the  form  of  (4)  holds  so  long  as  edge 
energies  may  be  Ignored.  The  critical  radius,  rc,  Is  given  by  an  equation 
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wnalogous  to  the  Kelvin  equation  vrltten  above  as  (2)  for  the  corresponding 
case  of  vapor-llquld  nucleatlon.  The  form  of  the  rc -equation  given  by 
Frenkel  (19**6)  la  r-p 

c ea;(x-T)  »> 

where  /O  Is  the  density  of  the  solid  phase  ttnd  Is  the  latent  heat  of 
fusion  of  the  substance  In  question.  It  has  been  argued  qualitatively  In  the 
fUacUrtSlon  "ucleatlon  above  that  the  rate  of  formation  of  nuclei 

(critical  embryos)  Increases  with  Increasing  degree  of  supercooling.  Super- 
ficial Inspection  of  (3)  might  seem  to  yield  a contradiction  to  this,  for 
the  nkT/h  prefactor  as  veil  as  the  factor  kT  In  the  exponent  both  tend  to 
lover  the  nucleatlon  rate  as  the  temperature  of  the  supercooled  liquid  falls. 
It  Is,  however,  In  spite  of  these  factors  and  because  of  the  sensitive  de- 
pendence of  Fc  on  T,  as  given  by  (U)  and  (5),  that  J Increases  rapidly  vlth 
decreasing  temperature.  Hence  spontaneous  crystallization  does  tend  to  be- 
come more  probable  the  greater  the  degree  of  supercooling.  The  problem  of 
cloud  physical  Interest  Is:  Hoes  this  probability  pass  through  some  sort  of 

threshold  value  near  -kOC  such  that  freezing  of  drops  becomes  nearly  Inevit- 
able within  periods  of  time  characteristic  of  cloud  processes?  Some  past  ef- 
forts to  answer  this  question  will  be  examined  next. 


3.  Previous  cloud  physical  applications  of  nucleatlon  theory 

The  present  position  of  the  theoretical  Investigation  of  homogeneous  nu- 
cleatlon of  supercooled  water  will  be  summarized  In  this  section  by  review lng 
tne  salient  features  of  the  principal  past  studies  of  this  problem.  In  the 
course  of  the  present  study,  significant  errors  have  been  found  5.n  each  of 
these.  Perhaps  the  major  contribution  of  the  present  investigation  has 
been  the  detection  of  end  at  least  partial  correction  for  these  errors. 

Krastanow  (191*!)  sought  to  apply  the  theoretical  work  of  Becker  and  Dor- 
Ing  (1955)  'ind  of  Volaer  (1959)  to  the  meteorological  problems  of  homogeneous 
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and  heterogeneous  nucleatlon  of  both  supersaturated  rapor  and  supercooled 
water.  Although  he  reached  conclusions  (e.g. , energetically  greater  proba- 
bility of  condensation  than  sublimation  in  rapor  below  OC)  whose  significance 
hare  not  always  recelred  the  attention  they  seem  to  deserve,  he  could  make 
no  deflnltlre  calculations  because  at  the  time  of  his  writing  the  rate -factor 
appearing  In  (5)  had  not  yet  been  adequately  eraluated.  Kraetanov  made  cer- 
tain estimates  of  A and  of  0^  which  will  be  discussed  below.  It  appears  not 
to  hare  been  prerloualy  noticed  that  Kras tan ow  failed  to  take  account  of  the 
appreciable  temperature  variation  of  the  latent  heat  of  fusion  of  water  (see 
Section  5 below)  In  his  calculations  of  . 

Fisher,  Hollomon,  and  Turnbull  (19b9),  recognizing  that  It  Is  re ry  dif- 
ficult to  determine  O5  by  experiment  or  theoretical  calculation,  used  nu- 
cleatlon theory  to  work  backwards  from  the  observed  -bOC  transition  temper- 
ature to  compute  0^  on  the  assumption  that  -bOC  is  in  fact  the  temperature 
for  which  J becomes  of  the  order  of  unity.  For  unstated  reasons,  they  omit- 
ted A In  the  nucleatlon  equation  which,  as  will  be  pointed  out  below  (Section 
9),  is  neither  qualitatively  nor  quantitatively  permissible.  Second,  they, 
like  Krastanow,  ignored  the  substantial  decrease  of  Lf  with  decreasing  tem- 
perature. Finally,  although  they  did  not  state  explicitly  what  mass  of 
supercooled  water  they  were  considering,  one  can  determine  this  by  solving 
for  their  N (the  n of  the  present  paper  multiplied  by  the  volume  considered) 
after  putting  their  computed  value  of  \j^  , 52.8  erg  cm”^, 
bark  Into  their  equation  along  with  all  other  numerical  values  they  specif^. 
One  finds  in  this  way  that  they  were  tacitly  considering  the  -hoc  (-580 
according  to  them)  transition  as  occurring  when  one  nucleus  formed  each  second 
somewhere  within  a sample  of  one  gram  of  supercooled  water.  This  was  a quite 
fallacious  basis  for  the  calculation  since  the  temperature  at  which  homogen- 
eous nucleatlon  occurs  will  clearly  vary  with  the  mee?  of  liquid  Involved, 
and  experimental  observations  yielding  the  -bOC  effect  concern  individual 
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_Q 

water  drops  vhose  masses  are  only  of  the  order  of  10  7 gnu  Because  of  these 
several  errors,  no  significance  can  be  attributed  to  their  estimate  of  (7^  • 
A corrected  calculation  of  the  same  type  a 111  be  made  belov  at  the  end  of 
Section  10. 

Lafargue  (1950)  made  a rather  curious  calculation  vhose  results  he  pre- 
sented as  proof  that  the  -40C  transition  Is  due  to  homogeneous  nucleatlon  of 
supercooled  drops.  He  began  by  pointing  out  that  x-ray  diffraction  studies 
of  the  structure  of  liquid  vater  have  revealed  that  each  oxygen  is  surround- 
ed tetrahedrally  by  four  other  nearest -neighbor  oxygens  at  2.76  Angstroms 
and  that  the  second -nearest  neighbors,  tvelve  In  number,  lie  at  a radial 
distance  of  14.53  Angstroms.  Next,  and  vltheut  further  comment,  he  Inserted 
rc  » U.53  Angstrom*  Into  the  Glbbs-Thomson  equation  (5)  along  vlth  numerical 
values  for  all  quantities  appearing  therein  except  T,  and  solved  for  T.  Els 
result  van  -4 1C,  In  excellent  numerical  agreement  vlth  Cvllong’s  ( 19*47)  ob- 
served value  for  the  transition  point.  He  concluded  that  this  argument  had 
demonstrated  that  the  transition  must  be  dependent  upon  homogeneous  nuclea- 
tlon. For  several  reasons  this  argument  of  Lafargue's  cannot  be  accepted. 
First,  although  the  fact  that  the  radius  of  the  Becond  sphere  of  coordina- 
tion In  liquid  vater  Is  about  4.5  Angstroms  Is  veil  established,  there  Is  no 
good  a priori  reason  vhy  one  should  Insert  this  particular  distance  rather 
than,  say  the  first  or  the  third,  or  even  some  higher  order  coordination  ra- 
dius Into  the  Cibbs-Thomson  equation  as  an  alleged  critical  radius.  It  Is  an 
epsent*H’.  feature  of  the  homogeneous  nucleatlon  process  that  rc  Is  not  a con- 
stant, txo  tuls  step  of  Lafargue's  assumes,  but  Is  instead  a quantity  vhlch 
decreases  vlth  decreasing  temperature . Second,  for  Lafargue  to  have  sought 
an  explanation  of  the  -40C  transition  In  terms  of  an  essentially  static  model 
rather  than  to  grapple  vlth  the  critically  Important  notion  of  the  rate  of 
nucleatlon  vas  to  Ignore  the  Intrinsically  kinetic  nature  of  the  nucleatlon 
process,  namely  that  at  all  times  there  exists  a spectrum  of  rates  of  cher.ce 
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formatlon  of  embryos  of  various  sices  from  two  molecules  on  up*  and  that 
spontaneous  crystalllratlon  only  occurs  vben  the  rate  of  formation  of  em- 
bryos of  critical  site  becomes  of  the  order  of  one  per  rolume  unit  -na  time 
unit  characteristic  of  the  experiment.  Lafargue  was  effectively  saying  that 
the  lce-llke  short-range  order  In  supercooled  water  waits  passively  for  the 
temperature  to  fall  far  enough  for  the  local  structure  to  simulate  an  Ice 
embryo  of  slie  sufficient  to  satisfy  the  Glbbs-Thomson  relation.  This  view 
Ignores  the  fundamental  role  played  by  statistical  fluctuations  In  the  nu- 
c lea t Ion  process.  Third,  Lafargue  used  the  value  of  Lf  corresponding  to  OC 
in  a calculation  referring  to  -41C.  Had  he  used  the  value  of  Lf  appropriate 
to  -4 1C  (see  Section  5)  a 1th  all  other  variables  unchanged,  he  would  have 
found  a transition  temperature  of  -58c  Instead  of  -4 1C,  and  the  agreement 
with  observation  would  have  disappeared.  It  would  seem  to  follow  from  these 
objections,  unfortunately,  that  the  problem  of  explaining  the  -40C  transition 
was  certainly  not  solved  by  Lafargue. 

A valuable  x-ray  diffraction  study  of  the  molecular  structure  of  super- 
cooled water  (Dorsch  and  Bcyd,  1951) » which  appears  to  have  been  stimulated 
In  part  by  Laf argue' s paper,  revealed  that  the  peak  In  the  angular  Intensity 
pattern  for  water  at  a scattering  angle  of  18  degrees  becomes  Increasingly 
better  resolved  as  the  degree  of  supercooling  Increases  dovn  to  those  Inves- 
tigators' lower  limit  of  observation  of  -16c.  Since  this  peak  Is  due  to  x- 
ray  scattering  by  the  second  sphere  of  coordination  at  about  4.5  Angstrom 
radius,  Dorsch  and  Boyd  concluded  that  the  structure  of  water  does  grow  In- 
creasingly ice -like  se  supercooling  proceeds,  as  suggested  a number  of  years 
ago  by  Bernal  and  Fowler  (1933).  They  pointed  out  that  this  trend  suggests 
that  at  the  still  greater  degree  of  supercooling  prevailing  near  the  -40C 
point,  the  structure  of  liquid  water  might  be  so  nearly  isomorphic  with  Ice 
as  tc  justify  Lafargue's  disregard  of  the  whole  concept  of  the  work  of  forma- 
tion of  an  Ice  embryo.  However,  the  work  of  formation  (properly  the  free 
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energy  of  formation)  of  a nucleus  goes  to  aero,  according  to  (4)  only  If 
goes  to  aero;  and  If  Isomorphism  developed  to  this  extent  near  -bOC,  then 
Lafargue's  calculation  vould  become  Invalid  for  still  another  reason,  since 
his  result  requires  a value  of  10.5  erg  cm"^  for  . There  seems  no  vay 
to  Justify  Lafargue's  approach. 

A still  more  recent  effort  to  account  for  the  -bOC  transition  In  terms 
of  homogeneous  nucleatlon  Is  due  to  Mason  (1952).  Unlike  Krastanow  vho  did 
not  have  available  the  quantum-statistical  expression  for  the  rate  factor  In 
the  nucleatlon  equation,  and  unlike  Lafargue  who  may  have  been  unauare  of  It, 
Mason  proceeded  from  Turnbull  and  Fisher's  equation,  (5)  above.  Finally,  un- 
like Fisher,  Hollomon,  and  Turnbull,  he  did  not  merely  use  the  equation  to 
make  an  estimate  of  0^  by  reasoning  backyard  from  the  observed  -bOC  transi- 
tion temperature  assuming  the  very  thing  vhlch  the  meteorologist  seeks  to 
confirm.  Instead,  Mason  made  a direct  attack  on  the  problem  by  estimating 
on  theoretical  grounds,  and  then  used  this  estimated  value  In  (3)  to 
study  the  behavior  of  J. 

In  several  of  the  follovlng  sections,  a number  of  refinements  of  detail 
In  Mason's  treatment  vlll  be  discussed.  Here  only  one  general  objection  must 
be  raised  to  the  basis  upon  vhlch  Mason  sought  to  demonstrate  agreement  be- 
tween theory  and  observation.  Avare  of  the  Inherent  uncertainty  In  his  the- 


oretical estimate  of  07  , Mason  found  that  values  of 


> T 


less  sensitive  to  uncertainties  In  than  were  values  of 


were 


lteelf 


and  n?  only  discussed  the  former  In  his  paper.  On  calculating  ILe  values 
of  that  derivative  sear  -hoc,  he  found  a nearly  t* n - t_h  in^r^sAe  in  .t  per 
Centigrade  degree  In  that  vicinity,  and  he  offered  the  rapidity  of  this  In- 
crease, per  ae,  as  the  explanation  of  the  experimental  observations  of  a sud- 


den transition  of  supercooled  clouds  In  the  neighborhood  of  -bOC. 
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Nov  a rapid  temperature  variation  in  the  nucleatlon  rate  near  -UOG  la 


certainly  necessary  to  any  statistical  explanation  of  the  observations,  but 
it  is  not  by  itself  sufficient.  One  must  ask  whether  the  nucleatlon  rate 
might  not  vary  Just  as  rapidly  at  some  temperatures  other  than  -1(0C;  and 
equation  (3)  yields  the  answer  that  it  does.  Indeed,  near  -300  the  rate  of 
formation  of  nuclei  is  found,  from  (3)  using  Mason's  figures,  to  change  by  a 
factor  much  greater,  250  per  degree;  and  at  -200  it  changes  by  a factor  of  one 
hundred  million  per  degree.'  Hence  the  mere  ten-fold  change  per  degree  near 
-hoc  can  certainly  not  by  itself  be  regarded  as  Implying  that  the  neighbor- 
hood of  -hoc  will  be  a preferred  region  for  transition.  Consequently,  it  be- 
comes logically  indispensable  to  consider  not  Just  the  steepness  of  the  tem- 
perature variation  in  the  nucleatlon  rate,  but  also  the  absolute  values  of 


that  rate. 

By  taking  this  latter  approach,  an  explanation  for  the  observed  phenomena 
might  have  been  presented  in  a rather  more  convincing  form  than  that  chosen  by 
Mason.  A further  Important  modification  in  approach  should  then  have  been  to 
evaluate  the  nucleatlon  rates  not  simply  in  terms  of  nuclei  per  cubic  centi- 
meter per  second  but  rather  in  terms  of  nuclei  per  drop  per  second/  The 
latter  is  the  pertinent  rate  here  since  each  drop  is  essentially  an  Isolated 
physical  system  in  the  process  in  question,  and  if  a given  drop  is  to  crys- 
tallize within  a reasonable  length  of  time,  t,  there  must  be  a high  probabil- 
ity of  at  least  one  nucleus  forming  somewhere  within  that  drop  during  the 
time  t.  Failure  o~  the  of  Fisher,  Hollomon,  and  Turnbull  ( 19^9 ) to 

j.ook  at  tuts  — *t GC  problem  in  this  sinner  has  already  been  cited.  A measure  of 


Mason  has  pointed  out  (private  communication)  that  his  original  approach  was 
exactly  this  one,  but  that  his  recognition  of  the  difference  in  sensitivity 
to  CT£  of  J and  of  } 1/^  T led  him  to  discuss  only  the  temperature  deriva- 
tives. In  so  doing  he  apparently  overlooked  the  fallacy  of  this  approach  oc- 
casioned by  the  even  steeper  J -variations  at  temperatures  warmer  than  -f*0C. 
And,  indeed,  a merely  large  value  of  ^ j/  ^ T occurs  at  -1*0C  for  almost  any 
helfway  plausible  numerical  values  one  may  insert  into  (3). 
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the  average  tine  between  nucleatlon  evente  within  each  drop  is  given  by  the 
reciprocal  of  the  nucleatlon  rate  per  drop.  Using  the  several  tentative  nu- 
merical values  suggested  by  Mason  for  the  parameters  entering  into  (3)>  one 
finds  for  this  "average  nucleatlon  time"  the  values  shown  in  Table  1 for 
drops  of  1,  5,  and  10  micron  radii.  It  is  In  considering  the  absolute  values 
of  these  tabulated  times  as  well  as  their  rapid  temperature  variation  that 
one  finds  quite  nice  apparent  agreement  between  observation  and  Mason's  work. 

Table  1.  Average  time  between  nucleatlon  events  within  water  drops  of 
radius  r at  various  temperatures.  Revised  from  Mason  (1952). 


Temperature 

r * 1 micron 

r ■ 5 microns 

r « 10  microns 

-36C 

187  days 

32  hours 

4.5  hours 

-37 

6.9  days 

80  minutes 

10  minutes 

-38 

11  hours 

5 minutes 

40  seconds 

-39 

42  minutes 

20  seconds 

2.5  seconds 

-40 

5 minutes 

2 second b 

0.3  seconds 

-41 

30  seconds 

0.2  seconds 

0.  Q3  seconds 

In  experiments  where  the  supercooling  is  produced  in  an  adiabatic  expan- 
sion, the  entire  nucleatlon  process  must  occur  in  times  of  the  order  of  sec- 
onds (unless  the  chamber  is  unusually  well  Insulated).  Hence  the  times  shown 
in  Table  1 suggest  that  no  crystals  should  be  observed  in  a cloud  of  drops 
whose  radii  are  of  the  order  of  microns  until  the  peak  expanslonal  cooling 
extends  down  to  near  -40C,  but  that  near  that  point  the  probability  of  nuclea- 
tlca  finally  rises  to  a level  at  which  every  drop  in  the  5-10  micron  range 
may  be  expected  to  experience  a nucleatlon  event  within  the  sensitive  time  of 
the  chamber.  In  a natural  cloud  where  the  period  during  which  a given  drop 
exists  at  a given  temperature  is  of  the  order  of  minutes,  one  would  expect, 
from  the  implications  of  Table  1 that  spontaneous  freezing  could  occur  at  a 
higher  temperature,  but  only  a degree  or  two  higher  because  of  the  extreme 
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temperature  sensitivity  of  nucleatlon.  Thus,  to  emphasize  the  basic  issue 
hers,  the  apparent  success  of  a statistical  theory  of  homogeneous  nucleatlon 
is  to  be  assessed  not  Just  on  the  basis  of  its  prediction  of  a rather  steep 
ascent  of  the  nucleatlon  probability  per  drop  in  the  -40C  region,  as  Mason 
tacitly  suggested,  but  much  more  on  the  basis  of  whether  it  predicts  that 
this  ascent  vlll  in  fact  carry  that  probability  to  values  Implying  that  a 
majority  of  drops  vlll  be  nucleated  at  least  once  in  the  time  Interval  char- 
acteristic of  the  cooling  process  involved.1* 

Judged  on  both  these  grounds,  the  version  of  Mason's  vork  shown  in  Table 
1 looks  nearly  perfect.  Even  the  slight  spread  of  the  experimentally  ob- 
served transition  temperatures  is  readily  understood  in  terms  of  Table  1 if 
different  observers  worked  with  clouds  of  slightly  different  drop-size  dis- 
tributions, and  particularly  if  the  sensitive  times  of  their  chambers  varied 
from  fractions  of  a second  to  several  seconds.  Both  types  of  dispersion  in 
experimental  conditions  are  entirely  plausible.  The  writer's  first  Interest 
in  nucleatlon  theory  was  aroused  by  these  implications  of  Table  1.  It  is 
with  some  regret,  then,  that  he  must  next  point  out  that  an  extensive  study 
of  the  basis  of  Mason's  calculations  has  uncovered  a number  of  defects  which, 
when  removed  in  the  best  way  recognized  by  the  writer,  destroy  the  neat 
agreement  between  Table  1 and  the  experimental  observations.  The  next  six 
Sections  will  be  devoted  to  a critical  examination  of  these  points  and  to  ef- 
forts to  improve,  in  every  way  possible,  the  basis  of  the  calculations  of  the 
dropwlse  nucleatlon  rates. 

She  calculated  rates  shown  in  Table  1 have  been  given  here  despite  their 
now  recognized  lack  of  validity,  partly  to  document  a logical  objection  to 
Mason's  mode  of  assessing  his  results,  but  much  more  to  illustrate,  for  the 
reader  not  familiar  with  nucleatlon  theory,  exactly  the  sort  of  numerical  be- 
havior of  the  dropwlse  nucleatlon  rates  which  one  would  hope  to  find  in 
searching  for  a statistical  explanation  of  the  -1*0C  transition. 
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4.  The  Olbbs-Thomaon  equation 

The  theoretical  nucleatlon  rate  Is  determined  Jointly  by  (3)*  (*♦  )>  and 
(5).  The  writer  has  been  unable  to  find  any  basis  for  modifying  (3)  or  (l*), 
though  (h)  vlll  be  Incomplete  If  ons  can  demonstrate  that  edge  energies  are 
being  Ignored  at  the  expense  of  numerical  accuracy.  The  Gibbs -Thomson  equa- 
tion, (5),  however,  appears  to  require  a slight  change.  Mason  followed 
Frenkel  (19^6,  p.  415)  In  using  the  form  (5);  but  an  examination  of  Frenkel's 
derivation  of  this  relation  reveals  that  It  is  only  approximately  correct  for 
such  large  values  of^T0  - T^as  one  encounters  in  the  present  problem.  From  a 
thermodynamic  argument  which  need  not  be  reproduced  here  (Frenkel,  1946,  p. 
368),  one  can  show  that 

( ▼ £ - )<Jp  * d(X/r) 

where  v^  and  v^  are  the  specific  volumes  of  the  liquid  and  solid  phases  re- 
spectively and  p Is  the  external  pressure  on  the  system.  Eliminating 
(v.  - v ) between  this  equation  and  the  differential  form  of  Clapcyron's 

V> 

equation,  and  Integrating  between  the  limits  r » rc  at  T and  r • 00  at  To  — 
273K,  one  gets  as  a more  generally  applicable  form  of  the  Olbbs-Thomaon  equa- 


tion. 


rt  = 





The  form  (5)  given  by  Frenkel  and  used  by  Mason  Is  an  approximation  to  (6) 
obtainable  by  expanding  ln(TQ/T)  Into  a series  and  using  only  the  first  term. 
For  only  very  small  degrees  of  supercooling,  negligible  error  is  Introduced 
by  using  (5);  but  when  as  here  - T becomes  40  Centigrade  degrees,  there  Is 
about  an  8 per  cent  error  In  (5),  which  moans  about  a 16  per  cent  error  in 
the  computed  value  of  Fc,  and  finally  a much  larger  error  In  J.  It  may  be 
well  to  state  here  that  though  this  elimination  of  an  error  of  less  than  10 
per  cent  In  rc  represents  a rather  slight  Improvement  when  viewed  against  the 
large  uncertainties  which  will  be  shown  to  exist  in  0^  in  Section  7,  still 
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lt  enters  (6)  factortaliy  and  not  addltlvely,  and  so  the  revision  Is  meaning- 
ful, A more  concrete  measure  of  the  significance  of  the  replacement  of  (5) 
by  (6)  is  found  in  the  fact  that,  keeping  all  other  parameters  unchanged  in 
Mason's  calculations,  this  modification  raises  J by  a factor  of  lO1*  for  drops 
at  -tOC,  and  thereby  raises  the  predicted  "spontaneous  freer lng  point"  by 
about  h Centigrade  degrees,  from  near  -hOC  as  shown  In  Table  1 to  about  -36c. 
That  Is,  the  change  from  (5)  to  (6)  already  alters  appreciably  the  sort  of 
close  agreement  with  observation  which  made  Table  1 seem  so  promising  an  ex- 
planation of  the  -tOC  transition. 


5.  Latent  heat  of  fusion  of  Ice 

It  has  been  pointed  out  above  that  previous  Investigators  of  the  problem 
of  homogeneous  nucleatlon  of  supercooled  water  have  failed  to  take  account  of 
the  temperature  dependence  of  the  latent  heat  of  fusion  of  Ice,  It  will  next 
be  shown  that  this  oversight  has  Introduced  a quite  large  error  Into  this 
previous  work,  Including  Mason's  work  as  modified  In  Table  1. 

A general  thermodynamic  relation,  sometimes  referred  to  as  Klrchoff's 
equation  (Olasstone,  19^7) , relates  the  lsobarlc  temperature  variation  of  the 
latent  heat  of  fusion  to  the  difference  In  specific  heats  at  constant  pres- 
sure for  the  liquid  phase,  Cpi,  and  the  solid  phase,  Cp8,  according  to 

/BU\  _ c _ c 

Tr  — n (7) 

'■  j I h .11  . 

For  water  substance  at  OC,  cp^  ^ 1 cal  gm  deg-1,  while  Cp8  ^ 0.5  cal  gm 
deg"l,  so  the  latent  heat  of  fusion  of  Ice  decreases  with  decreasing  tempera- 
ture at  a rata  of  about  0,5  cal  gm"*deg-*  Just  below  the  melting  point.  If 
this  rate  held  constant  down  to  -^OC,  would  decrease  to  about  60  cal  gm-1 
at  that  temperature.  Actually  the  situation  la  slightly  worse,  since  Cp8  Is 
known  to  decrease  with  decreasing  temperature  (Glaque  and  Stout,  1956),  and 
at  the  same  time  Cp^  Increases  with  decreasing  temperature  according  to  exper- 
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lmental  results  extending  to  -15C  (Scheel  and  He use,  1909)  plus  theoretical 
extrapolation  to  -50C  (Goff,  19*9)  based  on  thermodynamic  relations  between 
Cpi  and  de/dT. 

Values  of  Lf(T)  are  shown  in  Figure  1.  The  dashed  portion  of  the  curve 
in  Figure  1 below  -50C  is  the  writer's  own  simple  linear  extrapolation  be- 
yond Goff's  values.  Although  this  part  of  the  curve  is  a gross  extrapolation, 
it  will  be  used  in  certain  estimates  made  in  Section  10.  The  very  careful 
nature  of  Goff's  extrapolation  to  -50C  from  experimental  data  extending  to- 
-13C  warrants  considerable  faith  in  the  values  down  to  well  below  the  -40C 
transition  point,  which  is  fortunate  here.  It  is  clear  from  Figure  1 that 
previous  analyses  of  the  homogeneous  nucleatlon  of  supercooled  water,  in  which 
a value  of  80  cal  gm"1  has  repeatedly  been  used  for  Lf  in  the  Glbbs-Thomaon 
equation,  have  thereby  been  very  seriously  in  error  inasmuch  as  this  roughly 
^0  per  cent  overestimate  of  Lj  is  raised  to  the  second  power  and  then  used 
in  an  exponential.  Ftor  example,  if  just  this  one  correction  is  made  in 
Mason's  work  it  shifts  the  predicted  transition  point  from  near  -4oc  as  im- 
plied in  Table  1 to  the  very  much  lower  temperature  of  -68c.  Comment  has  al- 
ready been  made  on  what  the  same  type  of  error  did  to  Lafargue's  work.  In 
8ectlon  7,  below,  it  will  be  shown  that  failure  to  recognize  that  Lf  decreases 
with  decreasing  temperature  produced  still  another  type  of  error  in  Kr a sta- 
nce ' s (1949)  estimate  of  QT^  , 

There  is  a second,  less  obvious  correction  of  that  must  be  considered. 
This  second  correction  has  also  been  overlooked  in  previous  studies,  though 
fortunately  without  numerical  error  since  still  another  oversight  closely  com- 
pensates this  one,  as  will  be  shewn  later.  Despite  this  near-cancellation  of 
one  error  against  another,  it  seems  advisable  to  examine  each  here  to  avoid 
possibility  of  one  of  them  being  detected  and  corrected  and  not  also  the  other. 
The  present  point  hinges  upon  the  fact  that  the  ice  whose  fusion  is  involved 


Figure  1.  Temperature  dependence  of  the  latent 
heat  of  fusion  of  ice,  L^.  Values  from  OC  to  -50C 
after  Gcff.  Values  from  -50C  to  -80C  obtained  by 
linear  extrapolation. 


-22- 


in  homogeneous  nuc  lest  Ion  occurs  In  the  form  of  ex  trend/  snail  embryos  vlth 
radii  of  the  order  of  10  - 20  Angstroms.  for  particles  of  such  snail  radii 
of  curvature,  there  appears  the  sane  sort  of  geometric  effect  on  binding  en- 
ergy as  was  referred  to  above  In  a footnote  In  Section  2 for  the  case  of  liq- 
uid embryos.  The  binding  energy  (and  hence  the  latent  heat  of  fusion)  must 
decrease  vlth  decreasing  enbryo  radius,  and  although  this  decrease  Is  en- 
tirely negligible  for  radii  above,  say,  10~^cm,  It  becomes  numerically  sig- 
nificant for  particles  of  the  site  of  those  here  In  question. 

To  obtain  thermodynamically  a measure  of  the  reduction  of  Lf  due  to  sur- 
face curvature  effects,  consider  a spherical  Ice  embryo  of  radius  r from 
vhlch  a small  mass  dm  of  vater  molecules  Is  removed  by  melting.  Let  Lf  nov 
denote  the  latent  heat  of  fusion  for  r » oO at  the  temperature  In  question. 
Then  the  amount  of  energy  dQ  that  must  be  supplied  to  melt  off  dm  grama  of 
Ice  from  the  sphere  Is  less  than  Lfdm  by  the  amount  of  the  decrease  of  sur- 
face energy,  dA,  where  dA  is  the  change  In  surface  area  of  the  embryo  ac- 
companying the  spherically  symmetric  stripping  of  the  dm  grams.  Thus,  Inas- 
much as 

dm  - 4 7T  ^ r2dr 
and 

dA  = 8 7f  rdr, 

and  rince 

dQ  ■ Lfdm  - (j^  dA, 

It  follows  that 

ifr=  aa  - L,  -2 or  , (8> 

- am 

where  Lfr  Is  the  latent  heat  of  fusion  per  gram  of  Ice  at  a radius  of  curva- 
ture r.  It  may  be  noted  that  O5  Itself  Is  site  dependent,  but  allowance  for 
this  would  only  Introduce  second-order  correction  Into  (8), 
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The  numerical  magnitude  of  this  slit  correction  of  Lf  la  subject  to  tbe 
**m  relative  error  as  Is  the  value  assigned  to  # so  It  will  turn  out 
that  Its  precise  value  cannot  here  be  determined  since  0^  vlll  finally  turn 
out  to  be  uncertain.  However,  an  eatlaate  of  Its  effect  can  be  made  by  con- 
sidering a nuclear  radius  of  10  Angstroms  (approximate  size  in  water  nuclea- 
tlon  at  -40C)  and  assuming  20  erg  cm~^  (rough  average  of  three  estimates 
of  0^  given  In  Section  7 below).  In  this  case  the  correction  term  amounts 
to  about  13  per  cent  of  Lf.  This  correction  Is  thus  about  half  again  as  large 
as  the  correction  due  to  the  revision  In  form  of  the  Gibbs -Thomson  equation, 
which  was  shown  In  Section  4 to  shift  the  theoretically  predicted  transition 
temperature  by  about  4c. 

It  Is  Interesting  to  compare  (8),  which  has  been  obtained  here  thermo- 
dynamically, with  the  results  of  an  approx laate  molecular  argument  given  by 
Benson  and  Shuttlevorth  (19t>l).  They  calculated  the  ratio  of  the  molecular 
heat  of  vaporisation  for  a small  cluster  of  molecules  to  that  of  a plane  sur- 
face on  the  assumption  that  all  pairs  of  molecules  Interact  according  to  an 
Inverse  twelfth  power  repulsive  potential  plus  an  Inverse  sixth  power  attrac- 
tive potential.  For  the  s Imp lest  case  of  a central  molecule  surrounded  by 
twelve  close -packed  neighbors,  they  computed  the  ratio  to  be  about  0.4.  Since 
Benson  and  Shuttlevorth 's  model  Implies  that  this  ratio  may  be  expected  to 
vary  approximately  as  (1  - cH”1/^),  where  c Is  a constant  and  N Is  the  number 
of  molecules  In  the  cluster,  their  value  of  0.4  for  thirteen  molecules  may  be 
raised  to  about  0.6  for  a cluster  of  the  site  referred  to  In  the  preceding 
paragraph,  since  such  n cluster  would  comprise  23  to  30  water  molecules.  The 
corresponding  ratio  for  latent  heats  of  fusion  predicted  by  (8)  was  0.67* 

This  order  of  magnitude  agreement  with  Benson  and  Sbuttleworth's  result  Is 
surprisingly  good  In  view  of  the  present  uncertainty  In  and  In  view  of  the 
other  approximations  Involved  in  both  approaches  to  this  ratio.  It  seems  ap- 
propriate to  regard  this  agreement  In  size  effect  estimates  as  providing  mu- 
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tual  support  for  these  Quite  Independent  methods  of  estimating  latent  heat* 
of  small  aggregates  of  molecules. 

6.  Density  of  Ice 

The  coefficient  of  volume  expansion  of  Ice  near  OC  la  (Dorsey,  19*»0) 

1.5  x lO^deg."*.  It  follows  that  even  at  -50C  thermal  contraction  has  In- 
creased the  density  of  Ice  by  only  about  one  per  cent  of  Its  value  of  0.92  at 
OC.  This  change  Is  too  small  to  be  considered  here.  It  will,  furthermore, 
be  opposed  by  a decrease  of  density  associated  with  surface  distortion  of  the 
Ice  lattice  (Brown,  19**7)>  Since  nuclei  have  relatively  large  surfacc/vcluae 
ratios,  this  latter  type  of  density  change  may  become  sensible.  Shuttleworth 
(19^9)  calculates  that  It  amounts  to  several  per  cent  for  homopolar  crystals, 
but  this  density  reduction  Is  localized  almost  entirely  to  a surface  monolayer. 
Pending  further  developments  In  the  theory  of  dipole -bonded  crystals  like  Ice, 
the  effect  must  be  Ignored. 

7.  Surface  free  energy  of  a water-ice  Interface 

No  one  has  yet  devised  an  experimental  technique  for  directly  measuring 

for  an  Interface  between  supercooled  water  and  Ice,  nor  even  for  water 
and  Ice  at  the  triple  point,  to  the  writer’s  knowledge.  Nothing  could  be  more 
unfortunate  for  the  problem  at  hand.  Inspection  of  (3),  (M  and  (6)  reveals 
that  has  a controlling  Influence  on  the  homogeneous  nucleatlon  rate  since 
It  enters  to  the  third  power  In  an  exponential  In  the  J-equatlon.  It  follows 
that  major  effort  should  be  devoted  to  getting  Improved  estimates  of  by 
any  available  theoretical  means.  In  the  present  section,  attests  will  be  made 
to  refine  Krastanow's  ( 104 1)  and  Mason's  (1952)  values.  The  present  study  has 
not  led  to  any  essentially  new  basis  for  estimating  O5  > but  It  has  revealed 
that  numerically  significant  changes  have  to  be  made  In  both  of  these  existing 


methods 
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Krastanow's  basis  for  estimating  0^  depended  upon  a general  relationship 
rather  casually  suggested  by  Yolmer  (1939j  P»  181).  The  latter  proposed  that 


one  night  get  a rough  value  of  this  quantity  from  the  proportion 

0?  _ Gi 


L 


U 


(9) 


The  validity  of  such  a proportion  depends  on  hov  surface  free  energies  and 
latent  heats  are  related,  respectively,  to  Intermolecular  binding  energies, 
and  depends  particularly  on  whether  that  relationship  Is  Identically  the  same 
for  the  solid  and  the  liquid  states.  It  must  not  be  thought  that  this  Iden- 
tity has  been  proved  by  either  Volmer  or  Krastanow;  In  fact  neither  has  given 
any  comment  of  justification  for  (9)*  Since  Krastanow 's  estimates  of 
as  derived  from  (9)  have  been  used  without  apparent  question  In  subsequent 
studies  by  other  Investigators  (Lafargue,  1951;  Welckmann,  1951  )>  ud  will 
also  be  used  below,  It  becomes  a matter  of  some  interest  to  try  to  assess  this 
equation's  validity,  even  If  only  to  first  approximation. 

In  the  latter  spirit,  the  writer  offers  the  following  crude  derivation  of 
Volmer'8  equation  based  upon  a result  of  a theoretical  analysis  of  the  rela- 
tion between  surface  energies  and  latent  heats  of  sublimation  of  crystals 
(Shuttleworth,  19**9)«  From  direct  molecular  calculations  of  the  surface  en- 
ergy, Shuttleworth  found  that  an  expression  of  the  form  aB2  C £ holds 
for  homopolar  crystals,  where  a8  Is  the  Intermolecular  (or  Interatomic)  dis- 
tance In  the  solid,  Is  the  surface  free  energy  of  the  solid  against  Its 
vapor,  and  I Is  a constant  depending  upon  the  Interatomic  potential  used  and 
on  the  lattice  geometry  Involved.  For  a (111)  face  In  a face-centered  cubic 
homopolar  crystal,  K Is  about  0.2  for  one  plausible  potential.  Shuttleworth 
discussed  quantitatively  the  effect  on  0^  due  to  surface  distortion  of  the 
lattice  of  van  der  Waale  crystals,  but  could  only  point  out  qualitatively 


that  for  crystals  In  which  permanent  dipoles  Influence  the  bonding 
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will  be  appreciably  lower  due  to  greater  contribution  of  distortion  energy. 
Since  water  la  a substance  of  particularly  marked  dipole-dipole  bonding 
(hydrogen  bonding),  one  would  expect  a 1C -value  somewhat  below  0.2  in  the  case 
of  interest  here. 

If  one  next  considers  the  liquid  state,  he  observes  that  an  interaction 
picture  very  much  like  that  assumed  by  Shuttleworth  exists  as  far  as  concerns 
the  short-range  order  important  in  homopolar  or  dipole  bonding,  and  particu- 
larly so  in  the  case  of  supercooled  water,  which  appears  to  have  a very  ice- 
like  structure,  at  least  out  as  far  as  second -nearest  neighbors  (Dorsch  and 
Boyd,  1951)*  Hence  one  might  expect  that  supercooled  water  should  also  be 
characterlted  by  a Shuttleworth  relation  of  the  form  e^  /L^K',  where 
*JL  18  the  lnterm°l-ecular  distance  in  the  liquid  structure  and  K'  is  the  new 
constant  for  liquid  water.  But  since  to  within  the  limits  of  accuracy  of  the 
present  discussion  aB  and  since  structural  similarities  and  interaction 

similarities  should  imply  K££K',  it  should  follow  that  a proportion  of  the 


form 


j 1 1°) 

might  hold  approximately.  Furthermore,  each  ratio  in  (10)  would  be  expected, 
as  Shuttleworth  noted,  to  be  somewhat  smaller  than  0.2.  Actually,  the  right 
member,  which  is  readily  evaluated,  equals  0.08  at  OC^.  From  the  OC  value  of 

, , 85  , 

Lb,  one  finds  from  (10)  that  is  about^erg  cm  . Then  invoking  Antonov's 
rule  (Adam,  1951)  to  estimate  , one  finds 

0^  ~ — C7£  ^ (B5  - 75)erg  cm’2  ■ 10  erg  cm-2 

Finally,  putting  this  estimated  OC  value  of  0^  back  into  Volmer's  equation  (9), 
one  finds  that  the  left  and  right  members  have  numerical  values  of  2.93  x 10“^ 
and  5.01  x 10”9,  respectively,  when  the  latent  heats  are  once  again  expressed 
in  ergs  per  gram.  This  agreement  of  the  two  members  of  (9)  to  better  than  ore 


'One  must  use  the  lateut  heats  per  molecule  in  Shuttleworth' a relation. 
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per  cent  la  almost  certainly  fortuitous  because  of  the  cany  approx lms  t Ions 
Involved  In  the  preceding  discussion.  What  can  be  said,  however,  Is  that 
Yolmer's  equation  has  now  been  provided  with  at  least  some  theoretical  sup- 
port, which  has  not  been  true  previously,  to  the  writer’s  knowledge.  Briefly, 
one  can  regard  (10)  as  following  from  Shutt lew orth's  relation  plus  the  sim- 
ilarity in  structures  of  ice  and  supercooled  water;  and  Yolmer's  equation  (9) 
then  appears  to  be  related  to  (10)  through  Antonov's  rule.  It  should  be 
strongly  emphasised,  however,  that  (9)  has  been  justified  here  in  only  a very 
approximate  manner;  so  computations  based  upon  It  simply  cannot  yet  be  re- 
garded as  definitive,  a point  which  has  not  been  given  deserved  emphasis  in 
previous  applications. 

Since  (9)  has  now  been  rendered  at  least  plausible  as  one  relation  for 
making  approximate  calculations  of  (X,  the  next  step  Is  to  examine  whether 
Krastanow's  (1941)  applications  of  (9)  can  merely  be  taken  over  here  for 
later  use  In  estimating  J,  Although  Krastanow  did  not  state  the  numerical 
values  he  used  to  compute  his  0^  values,  the  writer  finds  that  Krastanow's 
results  can  be  reproduced  to  within  one  per  cent  by  Introducing  Into  (9) 
temperature -dependent  values  ofQ^  extrapolate*  below  OC  from  data  such  as 
those  given  by  Zemansky  (1943)  and  temperature -dependent  values  of  Ly  (List, 
1951)*  but  by  using  a constant  value  of  Lf  equal  to  80  cal  gm"*.  This  last 
feature  of  Krastanow's  approach  was  erroneous  because  of  the  temperature  de- 
pendence of  Lf . His  original  computed  values  of  are  shown  in  the  first 
line  of  Table  2.  By  using  Lf  values  taken  from  Figure  1,  one  obtains  the  set 
of  values  shown  In  the  second  line  of  Table  2. 


Table  2.  Water-Ice  Interfacial  free  energy  as  computed 

p 

from  Volraer's  equation.  Units  of  7^  are  erg  cm  . 


Temperature 

(°c) 

0 

-10 

-20 

• 

V* 

0 

-40 

-50 

Krastanow's 

10.0 

10.2 

10.4 

10.6 

10.8 

11.0 

Corrected 

,J2 

10.0 

9.6 

9.1 

8.5 

7.7 

6.8 
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Wbereaa  Krastanow's  values  Increase  slowly  w 1th  decreasing  temperature, 
the  temperature -corrected  values  decrease  appreciably.  That  the  latter  aort 
of  temperature  dependence  la  the  more  reasonable  would  seem  to  follow  from 
Dorach  and  Boyd's  (1951)  x-ray  diffraction  studies,  which  Indicate  that  as 
the  degree  of  supercooling  of  water  Increases,  the  liquid  becomes  Increasing- 
ly Ice -like  in  structure.  And  as  the  structures  of  the  two  phases  grow  In- 
creasingly more  similar,  It  should  follow  that  the  surface  free  energy  of  the 
Interface  between  the  two  phases  should  decrease  towards  the  tero  value  It 
must  exhibit  In  toe  limit  of  complete  isomorphism.  Consequently  the  present 
revision  of  Krastanow's  computation  must  be  regarded  as  more  nearly  correct 
than  the  original  values  and  will  be  used  later  In  Section  10.  Note  that  in 
the  neighborhood  of  the  -40C  point,  Krastanow's  value  Is  about  40  per  cent 
too  large,  If  the  present  revision  Is  taken  as  the  reference  value. 

Although  Krastanow's  estimates  of  were  substantially  in  error,  It 
must  next  be  noted  that  since  (j ^ and  Lf  enter  the  Gibbs  -Thomson  equation  only 
In  the  form  of  the  ratio  '^/l f,  and  since  exactly  this  ratio  forms  the  left 
side  of  (9),  Krastanow's  errors  lnO^  and  Lf  Inevitably  cancel  each  other  ex- 
actly as  far  as  **ny  calculations  of  rc  from  (5)  la  concerned.  Nevertheless, 

Fc  varies  as  CT$  rc 2>  80  Krastanow's  roughly  40  per  cent  overestimate  of  C3V 
Itself  near  -40C  would  still  remain  to  throw  off  any  nucleatlon  rate  calcula- 
tions based  on  his  results.  It  can  be  seen  that  some  extenuating  remarks 
might  be  made  along  the  same  line  about  Lafargue's  (1951)  calculation;  but 
Lafargue  mixed  a OC  value  of  Lf  with  a value  of  which  appears  to  have  beer 
Interpolated  linearly  between  Krastanow's  OC  value  and  his  -40C  value.  This 
somewhat  undoes  the  work  of  one  error  canceling  the  other,  so  the  point  will 
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not  be  amplified  here,  especially  since  Lafargue's  whole  approach  is  funda- 
mentally unsound. 

Having  examined  In  some  detell,  and  refined  In  one  way,  the  Volmer-Kras- 
tsncv  method  for  estimating  (5^*  attention  will  next  be  turned  to  a second 
method,  first  used  In  a meteorological  application  by  Mason  (1952).  Mason's 
approach  was  as  follows:  He  computed,  using  crystallographic  data  for  Ice, 

the  total  number  of  hydrogen  bonds  per  unit  area  of  a (0001)  plane  In  Ice, 
combined  this  with  an  estimate  of  the  energy  required  to  break  one  hydrogen 
bond,  and  therefrom  determined  the  cleavage  work  per  unit  area  required  to 
pull  an  Ice  crystal  apart  perpendicular  to  the  basal  plane.  Since  this  cleav- 
age process  creates  two  units  of  new  surface  area  per  unit  area  of  basal  plane, 
he  Identified  one -half  of  the  computed  cleavage  work  (one-half  of  204  erg 
cm”2  or  100  erg  cm"2)  with  l5^\  the  specific  surface  free  energy  of  Ice 
against  water  vapor.  Next  he  Implicitly  Introduced  Antonov's  rule  to  obtain 
( J ' by  subtracting  the  specific  surface  free  energy  of  water  against  vapor 
from  the  100  erg  cm"2.  Using  a value  of  80  erg  cm"2  for  the  surface  tension 
of  water  at  -40C,  he  obtained  22  erg  cm"2  forQ^*  Th*8  value,  then,  under- 
lies the  results  shown  In  Table  1 above.  In  view  of  toe  extreme  sensitivity 
of  calculated  J values  to  the  value  used  for  (5c, » It  Is  lramed late ly  disturb- 
ing to  see  that  Mason's  estimate  Is  about  three  times  as  large  as  that  which 
has  been  obtained  above  for  -40C  from  a refinement  of  the  Vclmer-Krastanow 
method.  This  discrepancy  demands  that  a careful  evaluation  of  Mason’s  method 
of  estimating  3^  he  carried  out,  for  although  It  has  been  emphasized  above 
that  the  Volmer-Kraetanow  method  remains  open  to  serious  question,  choice  be- 
tween It  and  Mason's  method  cannot  be  made  until  the  latter's  validity  Is  also 
assessed. 

A first,  and  rather  obvious  question  concerne  the  value  which  Mason  took 
for  (3^  at  -40C,  namely  80  erg  cm"2.  Perhaps,  although  he  did  not  so  state, 
Mason  may  have  been  Influenced  In  this  choice  by  a suggestion  due  to  Sander 
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and  Dsmk'ohler  (1943)  to  the  effect  that  a value  of  about  70  erg  cm*®  la  re- 
quired for  near  -40C  to  make  certain  experimental  data  on  water  vapor  nu- 
de at  on  agree  with  theory.  Bradley  (1951)  aa  well  aa  LaMer  and  Pound  (1951) 
have  suggested  that  the  decrease  In  below  values  extrapolated  from  exper- 
imental data  above  0C  may  actually  have  been  due  to  the  effect  of  high  sur- 
face curvature  of  the  water  nuclei  Involved.  This  explanation  seems  quite 
plausible  to  the  writer  (see  comments  at  the  end  of  this  Section  on  a similar 
phenomenon  In  Ice  nuclei).  If  true,  then  It  Is  definitely  Incorrect  to  make 
size  correction  In  (3^ and  then  use  this  In  Antonov's  rule  as  has  to  be  done 
In  Mason's  approach.  Fortunately,  some  observational  data  exist  to  throw 
light  on  this  question.  Hacker  (1951)  has  recently  carried  out  an  excellent 
experimental  study  of  for  water  supercooled  down  to  as  low  as  -22C.  His 
results  are  In  very  good  agreement  with  the  International  Critical  Tables  to 
as  far  as  the  latter's  data  extend,  -IOC,  and  confirm  the  existence  of  the 
slight  Inflection  point  In  (J^(T)  near  0C  that  can  be  discerned  In  the  I.C.T. 
values.  Hacker's  results  definitely  contradict  the  deductions  of  Sander  and 
Damkohler  (to  the  extent  that  the  latter  Investigator^’  conclusions  are  ap- 
plied to  plane  surfaces  of  water  rather  than  to  highly  curved  surfaces),  for 
near  -20C,  Hacker's  curve  Is  curving  slightly  upward  while  Sander  and  Dam- 


kohler ' s Is  curving  downward  and  lies  already  well  below  Hacker's  absolute 
values.  It  may  be  mentioned  that  there  Is  little  ground  for  questioning 
Hacker's  results  on  the  basis  that  surface  contamination  might  have  affected 
his  results,  since  contamination  can  easily  lower  surface  tensions,  but  al- 
most never  raises  them. 

An  extrapolation  of  Hacker's  experimental  curve  to  -40C,  preserving  the 
slight  curvature  of  the  OC  to  -22C  range,  yields  a value  of  about  84  erg  cm"  ; 
so  It  appears  necessary  to  alter  Mason's  value  of  80  erg  cm“^,  and  hence  In 
turn  to  revise  immediately  his  value  of  (j^  from  102  - 80  - 22  erg  cm  ’**  to 
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102  - 84  * 18  erg  cm"2.  This  rather  alight  revision  in  , it  must  be 
noted,  has  marked  effects  on  calculated  values  of  J and  raises  the  theoret- 
ically predicted  transition  temperature  by  about  6 Centigrade  degrees.  It  Is 
clear  that  rather  slight  changes  in  one's  basis  for  assigning  numerical  val- 
ues to  either  0^  or  0^  will  have  rather  large  effects  on  fltoce  the  lat- 
ter Is,  unfortunately,  obtained  from  Antonov ' s rule  as  a small  difference, 

Of ' -($£  in  two  large  terms. 

The  above  revision  has  sent  the  present  estimate  of  downward  to 

Improve  somewhat  the  agreement  between  tbe  Volmer-Krastanow  estimate  and  that 
obtained  from  the  cleavage-work  approach.  The  sensitivity  of  (5*,  to  the  cal- 
culated value  of  requires  that  the  latter  quantity  be  scrutinized  next. 

Hie re  appears  to  be  no  possibility  of  revising  Mason's  purely  crystallographic 
arguments  concerning  the  density  of  hydrogen  bonds  per  unit  area  of  a (0001) 
plane,  a density  which  can  be  expressed  inversely  as  an  area  per  bond  of  1.77 
x 10-1*  cm?;  but  two  objections  must  be  raised  against  his  method  of  calcu- 
lating the  the  energy  per  bond. 

Mason  follows  Pauling  (1940,  p.  304)  in  assigning  a value  of  4500  cal 
mole"l  to  the  hydrogen  bond  energy  at  OC,  but  then  states  with  no  explanation 
that  the  corresponding  value  at  -40C  should  be  about  51°0  cal  rnole”^-.  Now 
the  temperature  variation  of  the  heat  of  sublimation  (quantity  from  which  the 
hydrogen  bond  energy  is  determined  by  Pauling's  argument)  can  be  estimated 
from  Klrchoff’s  equation  applied  now  to  the  solid-vnpor  transition, 

(I r)e  = cr'r  ' 

where  cpy  and  cp£J  are  the  specific  heats  at  constant  pressure  for  the  vapor 
and  solid  phases  respectively.  Near  OC,  Cpy  - 0.443  cal  gm'^deg"*  and  Cpfl  = 
0.468  cal  gm’^deg"1,  so  I^/^T)pCi -0.025  cal  gm'^deg"1.  Both  cpy  and  CpB 
decrease  with  T,  but  in  such  a way  that  their  difference  remains  nearly  con- 
stant. Hence,  extrapolating  to  -40C  from  the  OC  value,  La  * 677  cal  gm"1, 
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one  finds  that  near  -40C  Lfl  678  eal  g*"1.  nils  change  represents  an  In- 
crease of  less  than  0.2  per  cent  over  the  OC  Talus,  in  serious  disagreement 
with  Mason's  roughly  15  per  cent  boost  In  the  hydrogen  bond  energy  orer  this 
te literature  interval. 

Since  the  discrepancy  between  the  above  thermodynamic  estimate  of  tem- 
perature effect  on  L_  and  Mason's  unexplained  adjustment  is  so  large  and  has 
a marked  effect  on  the  calculated  value  of  (5^  » It  is  fortunate  that  at 
least  a crude  independent  check  can  be  effected  as  follows*  The  coefficient 
of  expansion  of  ice  is  about  5 x 10"^  deg"*,  so  the  40  degree  temperature 
change  from  OC  to  -40C  will  decrease  all  lattice  distances  by  about  two  parts 
per  thousand  If  the  effects  of  the  slight  anisotropy  of  ice  are  ignored.  Now 
by  Badger's  rule  (Pauling,  1940,  p.  171),  the  force  constant  k0  for  a bond 
varies  about  as  d~3,  vhere  d is  the  pertinent  Interatomic  distance.  Since 
the  range  of  forces  here  Involved  is  only  of  the  order  of  d Itself,  the  work 
of  breaking  a bond  Is  of  the  order  of  l^d2,  or  d~3ds,  or  g-l#  Badger's  rule 
Is  only  an  empirical  relationship,  but  that  it  yields  tolerably  good  approxi- 
mations can  be  shown  by  checking  it  against  heats  of  dissociation  of  two  dif- 
ferent hydrogen-bonded  substances  having  known  bond  distances  (e.g.,  water 
and  formic  acid  dimer).  Combining  the  thermal  expansion  datum  with  Badger's 
rule,  one  predicts  that  the  bond  energy  in  ice  should  Increase  by  about  0.2 
per  cent,  Just  the  value  estimated  on  purely  thermodynamic  grounds  above. 

Since  Mason  gave  no  Justification  for  his  13  per  cent  increase  in  the 
bond  energy  associated  with  the  temperature  change  from  OC  to  -40C,  and  since 
two  independent  estimates  assign  It  a value  of  only  about  0.2  per  cent,  it 
appears  necessary  to  reject  Mason's  figure  of  5100  cal  mole*1  in  favor  of  the 
0C  Vilue  of  4500  cal  mole"*  for  the  hydrogen  bond  energy  per  mole  (simply  ig- 
noring the  here  predicted  0.2  per  cent  Increase  at  -40C).  This  revision  of 
Mason's  estimate  of  (5^  yields  a value  of  102(4500/5100)^90  erg  cm"2.  Ap- 
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p lying  Antonov ' a rule  and  using  (j£—  ^ «rg  cm"^  at  -1*0C,  one  now  gate  only 
6 erg  cb"2  for  a second  -revised  value  of  Q ^ . This  revision  has  sent  0^ 
below  the  revised  Volmer-Erastancw  value  of  7.7  erg  cm-2. 

However,  a further  error  appears  to  have  entered  Mason's  calculation  of 
, and  this  error  more  than  cancels  the  two  errors  now  cited  In  his 
calculation.  In  computing  the  work  that  must  be  expended  In  cleaving  an  Ice 
crystal  parallel  to  the  basal  plane,  Mason  took  account  of  only  a part  of  the 
total  bond  energy,  that  ascrlbable  to  hydrogen  bonding.  That  is,  the  value 
of  l 5 00  cal  mole*1  which  he  quotes  from  Pauling  was  obtained  In  turn  by  Paul* 
Ing  from  the  12,200  cal  male"*  heat  of  sublimation  of  Ice  by  observing  that 
only  about  three ‘fourths  of  the  total  bond  energy  In  an  Ice  crystal  Is  due 
to  hydrogen  bonding.  The  remaining  one-fourth  la  due  to  van  der  Waals  forces. 
Now  In  the  conceptual  experiment  of  cleaving  an  Ice  crystal,  work  would  be 
done  both  In  breaking  hydrogen  bonds  and  In  separating  neighboring  molecules 
against  the  van  der  Waals  attractions,  so  Mason  was  Incorrect  in  ignoring  the 
latter.  The  van  der  Waals  attractions  here  In  question  are  mainly  of  the  fc 
type  known  as  London  dlsperson  forces  and  arise  in  Ice  because  of  dlpole-dl- 
pole  Interactions  between  adjacent  oxygens.  Fluctuations  In  the  Instanta- 
neous dipole  moment  of  one  oxygen  nucleus  and  Its  surrounding  electron  cloud 
Induce  Instantaneous  dipoles  In  neighboring  oxygens  and  the  time  average  of 
these  fluctuating  dipole -Induced -dipole  Interactions  constitutes  an  attract- 
ive effect.  The  exact  nature  of  these  forces  Is  quite  pertinent  here  because 
all  of  the  available  electrons  In  water  may  be  regarded  as  surrounding  the 
oxygen  nuclei  so  that  It  becomes  clear  that  the  effective  areal  density  of 
van  der  Waals  bonds  on  a (0001)  plane  will  be  Identical  with  that  of  the  hy- 
drogen bonds  (since  we  may  Ignore  all  but  nearest-neighbor  Interactions). 

The  two  typc3  of  bonding  will,  in  fact,  act  coaxially  along  the  oxygen-oxy- 
gen axes.  Hence  Mason's  figure  of  1*500  cal  mole"^  has  only  to  be  replaced 
by  6100  cal  mole"l  g^ve  a estimate  for  of  >0(6l00/U500)=122  erg  cn"^. 
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Applylng  Antonov's  rule  as  before,  one  gets  (5^-38  erg  cm"2  ae  a final  re- 
vised value  resulting  from  Mason's  approach. 

Heglect  of  the  war  der  Waals  bonding  vas  the  moat  numerically  serious 
error  In  Mason's  work  and  la  seen  to  lead  to  a value  of  which  Is  about 
five  times  greater  than  the  revised  Yolmer-Krastanow  value  for  the  same  tem- 
perature. This  large  discrepancy  in  valuea  of  predicted  by  the  two  avail- 
able methods  leads  to  vastly  different  values  of  J near  -40C,  as  vlll  be 
shown  below,  and  demands  that  search  be  made  for  still  further  refinements  In 
the  two  calculatlonal  approaches  so  that  the  discrepancy  may  be  reduced.  It 
would,  In  particular,  be  very  helpful  even  to  know  If  the  revised  Yolmer- 
Krastanov  method  somehow  tends  to  underestimate  (5^  or  If  the  revised  Mason 
method  tends  to  overestimate  the  same  quantity.  A comment  will  next  be  made 
concerning  the  latter  possibility. 

In  estimating  the  cleavage  work  In  the  way  that  gave  Mason  a value  of 
102  erg  cm"2  for  Gj^and  which  has  here  yielded  122  erg  cm"2,  no  explicit 
allowance  has  been  made  for  the  reduction  In  the  work  of  separation  occasioned 
by  the  distortion  of  the  lattice  In  the  vicinity  of  the  two  newly  created 
(0001)  faces.  The  local  density  of  atoms  In  the  surface  layer  of  a crystal 
Is  slightly  lower  than  the  value  characteristic  of  the  deep  Interior,  and  as 
a result  a newly  created  surface  of  the  sort  envisioned  In  the  cleavage  opera- 
tion will  relax  Into  a state  of  lower  potential  energy  during  the  separation 
following  cleavage.  The  result  Is  (Brown,  1947;  Shuttleworth,  1949)  that  the 
net  cleavage  work  Is  reduced  by  an  amount  dependent  upon  the  degree  of  sur- 
face distortion  peculiar  to  the  crystalline  substance  In  question.  Shuttle- 
worth  has  succeeded  In  calculating  this  reduction  term  for  the  case  of  inert- 
gas  crystals  near  absolute  zero,  and  finds  that  It  is  of  the  order  of  only  a 
few  per  cent  of  the  gross  cleavage  work.  However,  Shuttleworth  observes  that 
for  a crystal  In  which  permanent  dipoles  contribute  to  the  total  bonding, 
this  reduction  term  will  become  appreciably  larger,  though  by  an  amount  which 
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he  Indicates  is  sot  easily  computed*  Since  the  lea  lattice  is  chiefly  held 

together  by  hydrogen  bonds  of  the  peraanent  dipole  type,  it  would  seem  to 

... ! 

folios  that  one  must  obtain  a too -large  value  of  0^  if  no  reduction  Is  ap- 
plied to  the  cleawage  work  to  allow  for  surface  distortion.  This  would  mean 
that  36  erg  cm"2  would,  in  turn,  be  an  overestimate  of  O5  , 

Although  the  writer  has  not  succeeded  In  making  any  precise  determination 
of  the  distortion  correction,  the  following  crude  argument  Is  offered  as  an 
Indication  of  the  order  of  magnitude  of  the  effect}  One  can  make  a theoret- 
ical calculation  of  CT^  by  the  cleavage-work  method  used  above;  but  to  get 
one  must,  of  course  base  the  calculation  on  Ly  rather  than  L8<  The  re- 

P JO 

suit  for  OC  Is  107  erg  cm  , which  Is  32  erg  cm  greater  than  the  known  val- 
ue of  75  erg  cm"2.  At  -Uoc  one  gets  112  erg  cm"2,  which  exceeds  by  28  erg 
cm"2  the  -1»0C  value  used  here  on  the  basis  of  Hacker's  study.  Since  the  sort 
of  surface  distortion  which  occurs  In  a crystal  lattice  also  occurs  In  a llq- 
uld  (see,  for  example.  Brown,  19^7),  one  may  take  30  erg  cm  as  an  average 
distortion  energy  for  liquid  water  over  the  temperature  range  of  Interest 
here.  That  this  Is  a percentually  very  much  larger  correction  than  the  dis- 
tortion correction  Shuttleworth  found  for  inert-gas  crystals  Is  not  surpris- 
ing In  view  of  the  very  different  force  laws  Involved  In  the  bonding  of  the 
two  types  of  substances;  but  whether  this  figure  may  be  regarded  as  a close 
estimate  of  the  distortion  correction  cannot  be  said  to  be  certain.  It  will 
be  a good  estimate  only  If  the  cleavage  calculation  upon  which  the  uncorrected 
values  of  were  based  do  take  appropriate  account  of  all  other  factors.  The 
writer  proceeds  on  the  assumption  that  the  latter  Is  so. 

Since  the  dipole  and  higher  multipole  binding  In  water  substance  Involves 
significant  Interactions  out  to  and  possibly  beyond  next -nearest  neighbors 
(Campbell,  1952),  one  must  expect  that  the  surface  structure  In  water  would 
be  somewhat  more  open  than  that  of  ice  since  In  water  only  short-range  order 
exists,  while  In  Ice  long-range  order  prevails.  Consequently  the  distortion 


energy  In  Ice  must  be  slightly  less  than  that  In  water,  and  roughly  in  In- 
Terse  proportion  to  the  latent  heats  of  vaporization  of  the  two  phases. 
Hence,  the  distortion  correction  in  ice  nay  be  taken  here  as  about  30(LvAs^ 
% 26  erg  cm-^  for  the  temperature  range  from  OC  to  -40C.  Applying  this  cor- 
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rectlon  to  the  previously  obtained  value  * 122  erg  cm  one  gets  9 C 
erg  cm-2  as  a crudely  distortion-corrected  estimate  of  the  specific  free 
energy  of  a vapor- Ice  Interface, 

Two  distinct  values  have  now  been  found  for  (tJ  from  the  cleavage-work 
approach,  122  erg  cm  as  a value  uncorrectsd  for  distortion,  and  96  erg 
cm  as  a value  very  roughly  corrected  for  distortion.  Using  Antonov's  rule 
with  values  of  (5^  taken  from  Hacker  ( 1951)  end  from  a smooth  extrapolation 
to  -50C  of  Hacker's  data,  one  gets  the  two  sets  of  0^  values  shown  In 
Table  3.  For  reference,  tho  first  line  contains  the  (T^  values  used  In  the 
calculations. 

— - 

Table  3.  Values  of  Os  In  erg  cm"  calculated  from  the  work  of  cleaving 
an  Ice  crystal  elong  a (0001)  plane. 


Temperature  (°C)  0 

-10 

-20 

-30 

-40 

-50 

Value 

of  C)£  used  (erg  cm’^)  75 

77 

80 

82 

84 

86 

uncorrected  for  distortion  energy  49 

45 

43 

4o 

38 

36 

"corrected"  for  distortion  energy  21 

19 

16 

l4 

12 

10 

A last  point  to  be  considered  under  the  heading  of  surface  free  energy 
of  a water-ice  Interface  concerns  the  question  of  the  size -dependence  of 
It  has  been  argued  theoretically  (e.g.,  To loan,  1948)  that  surface  energy 


must  In  general  be  a function  of  surface  curvature.  Indeed,  the  sort  of 
simple  molecular  model  In  terms  of  which  the  r-dependcnce  of  Lf  was  obtained 
earlier  here  demands  r-dependence  also  In  since  both  of  these  parameters 
are  largely  determined  by  surface  binding  forces  and  the  net  effect  of  the 
latter  Is  dependent  upon  surface  curvature.  LaMer  and  Pound  (1951)  have  con- 
sidered the  experimental  evidence  for  a decrease  of  surface  tension  of  liquid 
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water  with  decreasing  radius  of  surface  curvature  and  hare  discerned  what 
may  be  such  an  effect  In  the  experimental  data  of  Sander  and  Damkohler  (19U3). 
Bradley  (1940)  aade  a calculation  based  on  an  Inverse -seventh  power  law  of 
attraction  Inserted  Into  the  Lap  laden  theory  of  capillarity  and  concluded 
that  the  surface  free  energy  might  drop  by  slightly  less  than  10  per  cent 
at  radii  equal  to  ten  molecular  diameters!  and  by  perhaps  13  per  cent  for 
drops  with  radii  equal  to  only  five  molecular  diameters  (order  of  magnitude 
of  critical  Ice  embryos).  More  recently,  Benson  and  Shuttleworth  (1931) 
have  made  a similar  calculation  except  that  they  Included  repulsive  Interac- 
tions. They  found  that  for  clusters  of  as  few  aa  thirteen  molecules  the 
size -reduction  amounts  to  perhaps  13  per  cent  of  the  plane-surface  value. 

The  fairly  close  numerical  agreement  between  Bradley's  and  Benson  and 

Shuttleworth' e estimates  of  10-15  par  cent  for  the  particle  site  range  here 

Involved,  coupled  with  the  further  agreement  between  this  estimate  and  the 

writer's  a lie  correction  for  the  parameter  (Section  5 above)  has  led  him 

to  treat  0£/ as  size -Independent  In  the  Glbbs-Thomson  equation  (6).  This 

leaves  Just  the  explicit  factor  of  0^  In  (4)  to  produce  net  size  dependence 

In  the  nucleatlon  rates,  and  the  writer  chooses  to  employ  a uniform  10  per 

cent  reduction  In  whenever  Tc  Is  being  computed,  since  the  embryo  radii 

iYl 

are^Just  the  size  range  where  the  above  evidence  suggests  a correction  of 
this  order  magnitude.  Clearly  there  Is  room  for  improvement  here,  but  It 
does  not  seem  likely  that  much  will  be  possible  until  a more  complete  theo- 
retical treatment  of  surface  energy  of  small  aggregates  Is  forthcoming. 

8.  The  nuclear  shape  factor  g 

Once  all  parameters  entering  the  Glbbs-Thooson  equation  are  specified, 
the  value  of  Fc  Is  still  dependent  upon  the  shape  factor  g In  (4).  Mason 
chose  g 5S  23,  the  vaLue  appropriate  to  an  hexagonal  prism  whose  height  equals 
three  times  Its  apical  semidiameter,  but  he  did  not  discuss  the  basis  for 
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thlB  particular  choice.  Since  all  of  the  thermodynamic  arguments  that  under- 
lie the  Olbbs-lbomson  equation  assume  spherical  nuclei,  and  since  g for  a 
sphere  is  only  l47t 12.6,  the  use  of  any  other  value  calls  for  son*  Justifi- 
cation. As  soon,  hove rer,  as  one  gives  explicit  attention  to  this  matter  he 
recognizes  that  here  is  still  another  of  the  numerous  details  of  nucleatlon 
theory  vhlch  pose  fairly  subtle  problems.  In  considering  nuclear  shapes, 
one  must  cope  vlth  the  constraint  Imposed  by  the  lattice  geometry  character- 
istic of  ice;  and  for  such  small  crystallites  as  those  comprising  nuclei  in 
supercooled  vater  near  -bOC  it  is  not  clear  that  one  can  choose  arbitrarily 
either  a sphere  or  any  other  shape.  One  guiding  principle  here  is  that  the 
most  probable  nuclear  shape  should  be  that  vhlch  is  both  representable  by  a 
mlcrolattlce  and  is  consistent  vlth  the  requirement  that  the  total  free  en- 
ergy shall  be  a minimum  for  the  given  number  of  member  molecules.  By  anal- 
ogy vlth  gross  crystal  morphology  one  vould  expect  that  the  crystal  faces  of 
high  specific  surface  free  energy  vould  rapidly  grow  out,  leaving  best  devel- 
oped the  faces  of  low  specific  free  energy.  Since  Welckmann  (19*»7)  has 
clearly  established  the  dominance  of  elongated  hexagonal  prisms  near  the  -ItOC 
point,  one  concludes  that  near  that  temperature  the  basal  planes  must  be 
charactlzed  by  somewhat  larger  specific  surface  free  energy  than  are  the 
prism  faces.  It  follows  (really  by  direct  deduction  from  Welckmann' s obser- 
vations) that  the  preferred  shape  for  nuclei  near  *40C  ought  to  be  elongated 
prisms,  unless  curvature  effects  somehow  alter  this  problem  seriously.  Since 
Mason's  choice  of  g corresponds  to  a shape  that  might  be  built  up  from  the 
basic  lattice  geometry  of  ice  while  a sphere  does  not  (in  the  limit  of  aggre- 
gates of  nuclear  size),  and  because  furthermore  the  same  g seems  consistent 
with  Welckmann 's  observations  on  macroscopic  crystals,  the  writer  chooses  to 
follow  Mason  in  using  g£  23. 

As  for  the  implications  of  Inconsistency  this  choice  seems  to  introduce 
in  view  of  the  spherical  shapes  assumed  in  all  underlying  thermodynamic  argu- 
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tenta,  the  writer  would  make  the  follow lng  observationt  It  has  been  pointed 
out  earlier  that  (b)  is  correct  only  as  long  as  edge  energies  are  un important j 
and  the  ratio  of  the  contribution  of  edge  energies  to  that  of  surface  energies 
grows  large  for  erectly  the  size  limit  near  which  one  works  in  nucleatlon  the- 
ory. it  is  of  double  interest,  then,  to  note  that  in  Benson  and  Shuttle- 
worth's  (1951)  analysis  of  the  surface  energy  of  small  nuclei,  their  term, 
which  represents  the  contribution  of  edge  energy  to  the  total  surface  energy 
of  e polyhedral  cluster  car  be  shown  to  comprise  a term  of  identical  function- 
al form  and  approximately  equal  numerical  magnitude  as  that  which  has  been 
discussed  above  at  the  end  of  Section  7 in  connection  with  the  else -dependence 
of  G5  • That  is,  if  one  usee  a value  of  Q$  corrected  for  size -dependence, 
it  appears  that  he  is  essentially  making  an  edge-energy  correction  for  what 
is  in  microscopic  reality  a polyhedral  crystallite  rather  than  a true  sphere. 

9,  Activation  energy  for  ee If -d If fusion 

The  most  essential  difference  between  the  solid  phase  and  the  liquid 
phase  of  a given  substance  la  that  in  the  former,  long-range  order  exists  in 
the  form  of  a truly  crystalline  structure,  while  in  the  latter  only  e short- 
range  order  exists.  Yet  in  many  liquids,  of  which  water  is  an  excellent  ex- 
ample, this  short-range  order  is  well  enough  developed  that  one  seems  forced 
to  regard  each  molecule  as  being  rather  well  locked  into  a local  structure 
which  in  the  case  of  water  is  known  to  be  te  trailed  rally  bonded.  In  view  of 
this  tendency  towards  a crystal-like  local  structure,  molecules  cannot  move 
in  a truly  gas -like  fashion  among  each  other  in  a liquid.  Hence  the  phenom- 
ena of  self -diffusion  and  Viscosity  come  to  depend  upon  the  probability  with 
which  an  Individual  molecule  may  break  one  or  more  of  the  bonds  which  hold  it 
to  its  nearest  neighbors  preparatory  to  moving  relative  to  the  local  structure. 
The  energy  A required  to  thus  attain  the  more  mobile  state  of  higher  potential 
energy  la  drawn  from  the  thermal  energy  distribution  by  the  oechanlem  of  ran- 
dom vibrational  collisions,  end  the  fractional  probability  that  a given  liquid 


molecule  may  have  acquired  this  so-called  "activation  energy  for  self -diffu- 
sion" Is  given  by  the  factor  exp(-A/kT).  This  probability  factor  enters  Into 
the  nucleatlon  equation  because  each  molecule  (including  the  last  one  needed 
to  Just  attain  critical  size)  must  first  free  Itself  from  the  liquid  struc- 
ture near  the  water-ice  Interface  and  then  diffuse  from  the  liquid  region  over 
onto  the  crystal  lattice.  The  activation  energy  barrier  may  be  thought  of 
loosely  as  being  made  up  of  the  work  of  breaking  one  tntermolecular  bond  pri- 
or to  the  molecule's  chance  rotation  into  an  orientation  compatible  with  entry 
Into  the  local  lattice  plus  the  work  done  during  the  very  small  linear  trans- 
lation then  needed  to  lift  the  molecule  from  the  bottom  of  a potential  veil 
close  to  Its  former  liquid  ne ighbors  up  to  the  plateau  between  that  well  and 
the  potential  well  close  to  its  future  crystalline  neighbors.  In  the  limit 
of  very  low  temperatures,  where  Fc  becomes  negligibly  small,  the  rate  of 
growth  of  embryos  becomes  limited  almost  entirely  by  the  dlffualonal  barrier, 
and  a true  ice  lattice  may  have  to  give  way  to  a vitreous  form  of  the  solid 
phase . 

Since  exactly  the  same  activation  energy  enters,  In  an  Inverse  manner, 
Into  viscous  processes  It  becomes  possible  to  determine  A from  empirical  vis- 
cosity data.  Although  there  have  been  proposed  many  different  expressions  for 
describing  the  temperature  dependence  of  viscosity  (Partington,  1951)>  the 
most  widely  accepted  on  both  theoretical  and  observational  grounds  Is 


where  ^ Is  the  coefficient  of  viscosity  and  C Is  a constant. 

Mason  (1952)  made  use  of  (11)  to  determine  A using  viscosity  data  of  White  and 
Twining  (1913)*  He  obtained  the  value  A«3.3*10"^ergs,and  It  Is  of  Immediate 
Interest  to  note  that  this  Is  of  the  order  of  magnitude  of  the  energy  per 
bond  In  <at^r,  as  would  be  expected  on  the  basis  of  the  physical  picture  sug- 


gested above 
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It  should  be  mentioned  here  that  Erastanow  (1941)  gave  an  estimate  for  A 
of  2xi010  erg  gm~^  or  about  6x10"^  ergs  per  molecule,  vhlch  he  obtained  from 
a suggestion  due  to  Voloer  (1939)  that  A must  be  somevhat  less  than  the  mole- 
cular heat  of  vaporization  of  the  substance  in  question.  That  Erastanow 's 
value  Is  about  twice  as  large  as  that  found  by  Mason  from  (11)  Is  due  to  the 
fact  that  an  estimate  based  upon  (Erastanov  actually  used  300  cal  gm”*) 
Inevitably  measures  the  energy  needed  to  break  tvo  lntermolecular  bonds, 
whereas  diffusion  actually  appears  to  proceed  primarily  by  a sequence  of 
single -bond  ruptures  followed  by  rotations  about  the  re  mainly  bond,  even  In 
Ice  (Owston,  195 1)« 

To  complete  this  brief  examination  of  previous  treatments  of  A In  nucle- 
atlon  studies,  It  may  be  recalled  that  one  of  the  objections  raised  in  Section 
3 tc  the  calculations  of  FlBher,  Hollomon,  and  Turnbull  (1949)  concerned  their 
conqjlete  omission  of  a term  of  the  form  exp(-A/kT)  In  their  nucleatlon  equa- 
tion. In  terms  of  the  physical  picture  suggested  above  It  Beems  clear  that 
this  omission  left  out  an  effect  which,  especially  for  large  degrees  of  su- 
percooling, assumes  appreciable  importance  In  Inhibiting  embryo  growth. 
Finally,  It  may  be  noted  in  passing  that  a dlffuslonal  barrier  does  not  af- 
fect nucleatlon  of  either  liquid  or  solid  embryos  forming  In  supersaturated 
vapor  because  the  nearest -neighbor  distance  In  a vapor  is  generally  so  great 
as  to  preclude  the  sort  of  interactions  which  inhibit  diffusion  in  a liquid. 

In  estimating  A,  Mason's  approach  is  followed  here.  However,  Mason  did 
not  take  cognizance  of  the  fact  that  A Increases  with  decreasing  temperature 
due  to  the  Increasing  degree  of  order  In  the  liquid  structure  at  lower  tem- 
peratures. From  viscosity  data  given  by  Dorsey  (1940,  Tables  82,  85);  the 
writer  has  computed  A(T)  in  the  temperature  range  from  -10C  to  70C.  The  re- 
suita  are  shown  In  Figure  2 as  the  solid  portion  of  the  cunre.  Mason's  value 
of  3»3*10”13  erg  Is  seen  to  correspond  to  a temperature  of  about  5C,  while  In 
the  temperature  range  of  Interest  In  the  problem,  A Is  clearly  going  to  be 


Figure  2.  Temperature  dependence  of  the  activation  energy  for 
self-diffusion  in  water.  A la  given  In  units  of  ergs  per  molecule 
Values  from  60C  down  to  -IOC  calculated  from  empirical  viscosity 
data  (Dorsey).  Values  from  -IOC  down  to  -50C  obtained  by  linear 
eitrapolation. 
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largar.  In  Tien  of  the  fact  that  the  functional  form  of  (11)  Is  not  Indis- 
putably established,  the  writer  has  only  felt  Justified  In  making  a simple 
linear  extrapolation  dovn  to  the  temperature  llmltB  of  Interest  here.  This 
procedure  probably  tends  to  underestimate  slightly  the  activation  energy  bar- 
rier, but  Is  clearly  an  Improvement  over  the  use  of  a constant  value  corres- 
ponding to  a temperature  near  OC. 

10,  Calculations 

The  primary  objective  of  the  present  study  has  been  to  determine,  If 
possible,  whether  the  >^0C  transition  can  be  explained  theoretically  In  terms 
of  the  homogeneous  nucleatlon  process*  The  test  of  this  hypothesis  consists 
In  calculating  values  of  J from  (3),  ( ** ) , and  (6)  for  water  drops  of  speci- 
fied size  and  for  a number  of  different  temperatures  to  see  what  degree  of 
supercooling  Is  required  to  raise  the  nucleatlon  rate  per  drop  per  second  to 
the  order  of  unity  (cf.  discussion  of  Table  1 In  Section  3 abovej.  This  will 
now  be  done,  using  the  revised  numerical  values  of  the  several  parameters  dls 
cussed  In  Sections  3 through  9 above. 

In  the  calculations,  the  value  of  n In  (3)  will  be  taken  to  be  the  num- 
ber of  water  molecules  In  an  Isolated  drop  of  10  micron  radius.  From  the  val 
ues  of  J here  calculated,  one  may  readily  obtain  the  corresponding  values  for 
smaller  or  larger  radii  to  form  a table  of  the  form  of  Table  1,  since  J var- 
ies linearly  with  n.  A radius  of  10  microns  Is  selected  here  as  being  typi- 
cal of  the  drop  size  reported  In  laboratory  studies  of  the  -U OC  transition 
(e.g.,  Schaefer,  19^9).  It  also  serves  reasonably  well  to  represent  the  nat- 
ural cloud  physical  situation. 

Since  (3)  Is  transcendental  In  T,  one  cannot  simply  put  I°S^qJ  * 0 and 
solve  directly  for  T to  determine  the  theoretical  transition  temperature.  In 
stead,  one  may  compute  and  plot  log^jT)  and  then  read  off  the  value  of  T 
where  the  curve  crosses  the  axis  of  abscissas.  This  hae  been  done  here  In 


Figure  3.  Temperature  dependence  of  the  nucieatlon  rate  J for  super 
cooled  water  drops  of  10  micron  radius.  Curve  I was  computed  using 
0"c  ®B  estimated  by  the  revised  Volmer-Krastanow  method;  Curve  II 
was  computed  using  as  estimated  by  the  revised  Mason  method  with 
an  approximate  correction  for  surface  distortion  energy. 
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three  separate  series  of  calculations,  one  corresponding  to  each  one  of  the 
methods  of  estimating  C5^  • Th«  results,  displayed  in  part  In  Figure  5,  are 
as  follows: 

(1)  Uclr.g  (JT  *>i  obtained  by  the  revised  Vo  laser  -Era  sta  now  method  (Table 
2),  one  gets  curve  I of  Figure  3*  It  Is  seen  that  the  temperature  for  spon- 
taneous nucleatlon  of  10  micron  drops  Is  predicted  to  be  about  -IOC  according 
to  this  calculation. 

(2)  Next,  using  values  of  (5^  obtained  from  the  present  revision  of 
Mason's  approach  (Table  3)  one  obtains  values  of  J which  are  so  minute  as  to 
preclude  spontaneous  nucleatlon  anywhere  near  40C.  The  lowest  temperature 
to  which  the  writer  extended  this  part  of  the  calculations  was  -70C,  and  even 
for  that  extreme  degree  of  supercooling,  a drop  of  10  micron  radius  would  ex- 
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perlence  a nucleatlon  only  once  In  about  10TO  seconds,  if  this  second  set  of 
CTj  values  were  correct.  The  age  of  the  universe  is  believed  to  be  of  the 
order  of  lO1^  seconds.  Such  low  nucleatlon  rates  cannot  conveniently  be 
shown  In  Figure  3,  and  would  clearly  be  of  no  meteorological  Interest  cnyvay. 

(3)  Finally,  using  the  values  of  based  on  the  present  revision  of 
Mason's  approach  but  corrected  very  roughly  for  distortion  effects  In  the  man- 
ner indicated  earlier  (Table  3),  one  gets  the  curve  shown  as  II  In  Figure  3, 
from  which  the  theoretically  predicted  transition  temperature  Is  found  to  be 
about  -26c. 

Before  discussing  these  results,  one  additional  calculation  will  be  ex- 
amined. Fisher,  Hollomon,  and  Turnbull  (19^9)  have  used  (3)  to  compute  the 
value  which  (J^  must  have  If  the  -UOC  transition  really  Is  due  to  homogeneous 
nucleatlon.  In  Section  3 above,  several  objections  to  their  calculation  have 
been  noted.  It  Is  of  Interest  here  to  repeat  that  calculation  on  the  basis  of 
such  revisions  as  have  been  made  In  the  course  of  the  present  study  (excepting, 
of  course,  revisions  In  methods  of  calculating  Itself).  Assuming  that  the 
laboratory  observations  of  the  -40C  transition  apply  to  drops  of  10  micron  ra- 
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dlus,  on#  finds  by  this  sort  of  lnrerse  ealoulstlon  that  a Talus  of  O5  of 
16  erg  cm"2  would  bring  the  present  theory  Into  accord  with  laboratory  obser- 
vations. This  figure  Is  to  be  compared  with  the  -J*0C  values  of  7«7>  38>  and 
12  erg  cm"2  obtained  In  this  study  by  three  different  methods,  and  with  the 
value  of  33  ®rg  cm-2  deduced  by  Fisher,  Hollomon,  and  Turnbull. 

11.  Discussion. 

On  the  basis  of  a rather  large  number  of  revisions  of  previous  efforts 
to  examine  the  -40C  problem,  three  separate  calculations  have  been  made  hers 
employing  three  different  sets  of  estimated  values  of  , the  variable  hav- 
ing dominant  numerical  Influence  on  the  theoretical  nucleatlon  rates. 

The  first  of  these  calculations  (Curve  I,  Figure  3)  Implies  a far  too 
efficient  nucleatlon  process.  If  all  of  the  other  parameters  controlling  J 
can  be  trusted  (and  they  can  certainly  be  trusted  to  far  greater  extent  than 
can  (5^),  then  one  can  conclude  that  the  values  of  O5  estimated  by  the  7ol- 
mer-Krastanov  method  must  surely  be  lover  than  the  correct  values  because  It 
Is  veil  known  that  liquid  water  drops  often  exist  In  clouds  at  temperatures 
substantially  below  -IOC.  It  must  be  recalled  that  no  really  firm  basis  for 
accepting  the  Volmer-Krastanov  method  for  estimating  has  yet  been  given, 
although  the  writer  has  suggested  earlier  here  one  way  In  which  it  can  at 
least  be  rendered  qualitatively  plausible. 

The  second  calculation  predicts  an  effectively  zero  nucleatlon  rate  down 
to  degrees  of  supercooling  well  beyond  any  of  meteorological  Interest.  That 
Is,  when  the  writer  calculates  (y ' after  the  manner  suggested  by  Mason,  but 
with  Mason’s  apparent  errors  corrected  In  the  best  way  recognized  by  the  writ- 
er, the  Implied  nucleatlon  rates  are  found  to  fall  completely  to  account  for  a 
-UOC  transition.  This  plus  the  results  of  the  first  calculation  seems  to 
place  the  burden  of  the  argument  squarely  on  tbe  correction  for  distortion  en- 
ergy In  the  cleavage  computations. 
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In  the  third  calculation,  where  a crude,  though  at  least  not  arbitrary 
correction  for  the  lattice  distortion  effect  has  been  Included,  the  predicted 
transition  temperature  still  falls  outside  the  range  of  observed  transition 
points,  though  by  an  amount  which  Is  only  about  half  as  great  as  the  discrep- 
ancy for  the  first  of  the  three  calculations.  The  phenomenal  Increase  In 
computed  nucleatlon  efficiency  that  accompanies  the  application  of  the  dis- 
tortion energy  correction  Is  very  disquieting  In  view  of  the  uncertain  grounds 
on  which  a distortion  correction  was  made  here.  Consequently  It  becomes  Im- 
portant to  examine  the  possibility  that  one  could  make  some  more  straightfor- 
ward attack  on  the  problem  of  determining  the  distortion  energy. 

If  one  had  experimental  data  concerning  the  lattice  expansion  in  the  out- 
er two  or  three  molecular  layers  near  an  ice  surface,  be  might  calculate  the 
distortion  energy  fairly  accurately  from  the  known  compressibility  coefficient 
of  Ice;  but  unfortunately  x-ray  data  Inevitably  provide  only  a picture  of  the 
average  structure  down  to  many  tens  of  molecular  or  atomic  distances.  Per- 
haps electron  diffraction  methods  hold  somewhat  more  promise,  but  the  diffi- 
culty In  adequately  preparing  a sample  surface  in  the  case  of  Ice  would  prob- 
ably be  a limiting  factor  here. 

On  the  theoretical  side,  it  might  be  hoped  that  a direct  calculation  mod- 
eled upon  that  made  by  Shuttleworth  for  Inert-gas  crystals  would  be  the  answer, 
and  this  Is  Indeed  so  In  principle.  In  fact,  however,  such  a calculation 
would  be  extremely  difficult.  Following  Shuttleworth,  one  would  seek  a gener- 
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interplane  o '.stance,  solve  for  the  implied  equilibrium  interplane  distance, 
and  then  Insert  this  back  Into  the  general  energy  expression  to  compute  the 
reduction  In  energy  due  to  surface  distortion.  But  whereas  the  Interatomic 
attractive  potential  In  an  inert-gas  crystal  falls  off  as  r"^,  that  for  di- 
polar ice  would  contain  terms  falling  off  only  as  r"^,  so  one  would  have  to 


iinlmize  Jclrtly  the  potential  energy  of  at  least  the  tvo  uppermost  planes 
of  water  molecules,  and  this  might  be  only  a first  approximation.  This  com- 
plication vould  not  be  so  serious  by  Itself,  but  each  of  these  potentials  In- 
volved would  ha  7C  tC  be  cxprs  ased  aa  a sum  over  the  lattice  of  a Taylor  ser- 
ies expansion  (multipole  expansion)  of  the  charge  distribution  around  each 
HgO  molecule.  Campbell  (1952)  has  examined  the  problem  of  calculating  the 
lattice  energy  of  Ice  (for  fixed  lntermoleculer  distances)  in  terms  of  mul- 
tlpoles,  and  has  found  that  even  when  one  includes  up  to  fifth-order  terms 
(e.g.,  octupole-octupole  Interactions)  for  nearest -neighbors  the  agreement 
vlth  the  corresponding  thermodynamic  data  Is  unsatisfactory.  Furthermore, 
Campbell  found  that  Interactions  between  next-nearest  neighbors  are  still 
significant  out  to  beyond  third -order  terms.  From  Campbell's  Table  XI , one 
can  see  that  this  means  that  the  potential  energy  for  each  molecule  In  a 
plane  near  the  surface  must  contain  some  fifteen  terms  In  the  variable  Inter- 
plane  distance;  and  since  one  must  find  a Joint  minimum  for  at  least  the 
first  tvo  Interplane  distances,  It  follows  that  at  least  thirty  terms  are  In- 
volved In  the  function  whose  minimum  vould  be  sought.  Finally  a very  serious 
further  complication  enters  by  virtue  of  the  fact  that  these  multipole  Inter- 
actions, unlike  the  dispersion  forces  (London,  1937)  vlth  which  Shuttlevorth 
dealt,  are  not  simply  additive,  so  the  type  of  lattice  sums  which  were  usable 
In  Shuttlevorth' s calculation  are  not  applicable  here. 

In  all,  it  seems  questionable  whether  the  heroic  efforts  that  vould  have 
to  be  made  to  effect  this  lattice  calculation  of  the  distortion  energy  are 
meteorologically  Justifiable.  This  last  step  required  to  complete  the  theo- 
retical exploration  of  the  problem  of  homogeneous  nucleatlon  of  supercooled 
water  drops  appears  to  be  of  an  order  of  difficulty  far  exceeding  that  of  any 
solid-state  calculations  that  have  yet  been  carried  out  for  Ice  or  other  dipole 
crystals.  Consequently  the  writer  feels  that  the  fact  that  even  a crude  esti- 
mate of  the  distortion  has  here  yielded  a theoretically  predicted  temperature 


of  -26c  may  for  the  present  be  taken  as  a strong  indication  that  the  -40C 
transition  Is  truly  an  effect  of  homogeneous  nucleation,  an  effect  whose  pre- 
cise nature  cannot  be  specified  quantitatively  in  the  present  state  of  know- 
ledge of  the  solid-state  physics  of  ice.  Careful  scrutiny  of  the  present  re- 
visions of  previous  investigators'  work  will  have  to  be  made  by  others  before 
this  conclusion  can  be  accepted;  but  if  the  present  calculations  can  be  re- 
garded as  essentially  correct  up  to  the  last  step  of  making  the  distortion 
correction,  and  if  no  experimental  evidence  clearly  contradicting  the  homo- 
geneous nucleation  hypothesis  is  forthcoming,  then  there  would  seem  to  be 
strong  enough  grounds  for  concluding  that  -40C  is  the  temperature  to  which 
water  droplets  of  cloud -particle  size  must  be  supercooled  in  order  that  there 
shall  be  nearly  unit  probability  of  the  formation  of  an  ice  embryo  of  critical 
size  somewhere  within  each  drop,  that  is,  that  the  cloud  there  undergoes  spon- 
taneous nucleation. 

13.  A remark  on  the  icing  of  aircraft 

The  theory  of  homogeneous  nucleation  may  shed  some  light  on  one  important 
aspect  of  the  aircraft  icing  process.  In  those  dlccue«ions  of  icing  with 
which  the  writer  happens  to  be  familiar,  no  clear  explanation  seems  to  be  of- 
fered for  the  reason  why  the  accreted  supercooled  water  freezes  after  deposi- 
tion even  though  it  may  have  remained  liquid  for  a long  period  prior  to  en- 
trance of  the  aircraft  into  the  given  cloud.  Or,  if  an  explanation  is  given, 
some  allusion  is  made  to  the  shock  of  the  Impact  and  the  reader  is  reminded  of 
the  allege!  role  of  mechanical  disturbance  in  initiating  freezing  in  bulk  su- 
percooled water. 

If,  as  the  writer  has  here  been  attempting  to  show,  the  existence  of  the 
supercooled  cloud  is  primarily  due  to  the  negligibly  small  rate  of  homogeneous 
nucleation  above  -4f)C  then  the  reason  why  water  freezes  after  deposition  though 
not  before  is  that  the  free  energy  barrier  to  the  formation  of  a critical  era- 
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bryo  on  vhich  the  crystallization  may  begin  Is  no  longer  an  obstacle  to  freez- 
ing vhen  the  drop  Impinges  on  the  surface  of  the  airfoil.  If  a sheet  of  Ice 
Is  already  present  on  the  airfoil,  the  Impinging  supercooled  vater  merely 
*~inc  tbnt  Ice  lattice  In  the  ordinary  vay  In  vhlcb  vater  freezes  on  lattices 
of  email  curvature.  If  Ice  does  not  preexist,  as  at  the  monaat  of  first  pen- 
etration Into  the  Icing  region,  then  the  microscopic  roughness  of  the  airfoil 
surface,  or  perhaps  some  adsorbed  forelgh  material,  Is  almost  certain  to  Ini- 
tiate freezing  by  heterogeneous  nucleatlon  aomevhere  v 1th In  the  vater  film 
covering  the  leading  edge.  Once  a microscopic  lattice  Is  created  anywhere 
along  the  vetted  surface  of  the  airfoil,  freezing  vlll  spread  epidemically 
throughout  all  portions  of  the  vater  film  continuous  vlth  the  locus  of  nucle- 
atlon. After  that,  the  freezing  of  Impinging  drops  occurs  by  the  sort  of 
process  first  outlined. 

The  essential  difference  betveen  the  film  of  vater  on  the  airfoil  and 
the  same  mass  of  vater  In  the  form  of  many  cloud  drops  Is  that  In  the  former 
case  only  one  successful  nucleatlon  event  need  occur  vlthln  its  total  volume 
to  produce  complete  crystallization,  vhlle  In  the  latter  case  some  supercool- 
ing vlll  persist  until  there  has  occurred  one  nucleatlon  event  vlthln  each  of 
the  Individual  drops.  That  la,  the  difference  Is  essentially  topological, 
having  to  do  vlth  the  connectedness  of  the  vater.  If  It  la  Indeed  true  that 
average  cloud  drops  cannot  be  expected  to  be  nucleated  once  vlthln  their  full 
lifetime  of  the  order  of  minutes  at  temperatures  much  above  -bOC  (recall  gen- 
eral Implications  of  Table  1),  then  one  can  understand  why  aircraft  Icing  Is 
observed  above  that  temperature  but  Is  almost  unknown  below  that  temperature, 

13.  Summary 

Previous  Investigations  of  the  -40C  transition  In  supercooled  water  drops 
have  been  shown  to  contain  a variety  of  inaccuracies  which  have  had  marked  ef- 
fect on  predicted  nucleatlon  rates.  After  a number  of  modifications  were  car- 
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rled  cut  In  existing  approaches  to  the  homogeneous  nucleatlon  problem,  three 
different  estimates  of  the  specific  surface  free  energy  of  a nateT-ice  Inter- 
face were  used  to  estimate  the  temperature  dependence  of  the  dropwlse  nuclea- 
tlon rate.  Of  the  three  results,  the  one  In  which  the  most  confidence  can  be 
placed  yielded  a predicted  transition  temperature  substantially  too  high, 

-26C,  but  was  based  on  a correction  for  surface  distortion  effects  that  Is 
recognised  to  be  quite  crude.  There  seems  little  Immediate  hope  for  gaining 
Improved  precision  In  the  estimate  of  the  dlatortlon  correction  on  theoreti- 
cal grounds,  since  the  required  calculation  poses  very  formidable  difficulties. 

The  fact  that  one  rough  estimate  of  the  distortion  energy  has  led  to  a 
predicted  transition  threshold  even  as  close  to  the  observed  value  as  -26c  Is 
tentatively  taken  to  Indicate  that  the  -40C  transition  Is  an  effect  of  homo- 
geneous nucleatlon  whose  precise  explanation  will  have  to  await  further  devel- 
opments In  the  solid-state  physics  of  Ice. 


14.  Suggestion  for  future  research 

It  would  be  of  the  greatest  Interest  to  have  some  sort  of  experimental 
determination  of  the  surface  free  energy  of  a water-ice  Interface.  Although 
a direct  evaluation  of  this  parameter  for  supercooled  water  seems  quite  out 
of  question,  It  may  not  be  entirely  Impossible  to  determine  this  quantity  at 
the  triple  point  of  water  by  employing  sufficiently  elaborate  thermostatic 
control.  It  should  be  clear  from  the  previous  discussions  of  the  distortion 
energy  correction  that  Immeasurable  Improvement  In  one's  understanding  of  the 
nucleatlon  problem  would  result  even  from  a reasonably  accurate  determination 
oi  at  the  Iplc  point. 
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II.  On  the  Electrical  Conductivity  of  the  Lover  Stratosphere 

Abstract  --The  reality  of  the  decrease  of  conductivity  vlth 
height  observed  at  the  top  of  the  Explorer  H flight  Is  supported 
by  some  apparently  overlooked  measurements  made  earlier  by  Idrac. 

The  cause  of  this  decrease  is  discussed  and  an  hypothesis  of  con- 
vective updraft  of  Altken  nuclei  from  troposphere  to  stratosphere 
Is  examined.  The  hypothesis  appears  Incapable  of  accounting  for  a 
steady-state  vorldvlde  population  of  stratospheric  nuclei  suffici- 
ent to  satisfy  existing  observations.  It  Is  urged  that  Idrac's 
measurements  be  repeated  on  a more  extensive  basis  using  more  mod- 
ern sounding  techniques. 

1.  Introduction 

One  of  the  many  results  of  the  1933  stratosphere  balloon  flight  of  the 
Explorer  II  vas  the  observation  of  a surprising  decrease  of  atmospheric  elec- 
trical conductivity  vlth  height  through  the  top  fev  kilometers  (19-22  km)  of 
that  flight.  Gish  and  Sherman  (1936)  have  discussed  this  feature  of  the 
sounding  and  Gish  (1939)  has  made  some  suggest  Iona  as  to  Its  possible  ori- 
gin; but  no  further  attention  seems  to  have  been  given  this  matter  until 
Holzer  and  Saxon  (1952)  recently  examined,  on  theoretical  grounds,  the  cur- 
rent distribution  that  may  be  expected  to  exist  above  and  around  an  active 
thunderstorm.  They  approached  this  problem  In  an  effort  to  check  the  signif- 
icance of  the  important  thunderstorm  electrical  measurements  made  by  Gish  and 
Walt  (1950). 

Holzer  and  Saxon  employed  the  assumption  that  atmospheric  conductivity 
Increases  exponentially  and  hence  monotonlcally  vlth  Increasing  height  through 
the  troposphere  and  stratosphere  all  of  the  vay  up  to  some  conducting  layer  In 
the  lover  Ionosphere.  In  making  this  assumption,  these  authors  did  not  over- 
look the  Gish  and  8herman  observations  of  a conductivity  minimum,  but  they  did 
choose  to  omit  this  feature  from  their  analysis  on  the  basis  that  It  Is  not 
knovn  vhether  this  represents  a commonly  occurring  condition  of  the  lover 
stratosphere  or  vhether  It  vaa  an  anomaly  peculiar  to  the  Explorer  II  sound- 
ing. Holzer  and  Saxon  point  out  that  a shallov  layer  of  lov  conductivity 


would  not  alter  their  principal  conclusions  concerning  the  lonopberlc  destin- 
ation of  the  currents  measured  by  Gish  and  WAlt  ( 195 0 ) , but  emphasize  that 
such  a layer  could  have  a marked  effect  on  certain  other  atmospheric  electric- 
al phenomena,  notably  the  vertical  field  fluctuations  observable  at  the  earth's 
surface  at  distances  of  the  order  of  many  tens  of  kilometers  from  active  thun- 
derstorms. 

2.  Earlier  evidence  for  a stratospheric  conductivity  decrease 

Holzer  and  Saxon's  discussion  has  pointed  up  the  uncertainty  as  to  whether 
a stratum  of  low  conductivity  Is  a cowaon,  or  Indeed  even  a real  feature  of 
the  stratosphere  and  has  prompted  the  present  writer  to  call  attention  to  an 
earlier  and  apparently  forgotten  study  of  the  electrical  state  of  the  lower 
stratosphere  In  which  there  was  found  evidence  for  a decrease  of  conductivity 
similar  to  that  found  In  the  Explorer  II  sounding.  Idrac  (1926)  made  a num- 
ber of  balloon  soundings  of  the  vertical  electric  field  Intensity  over  Trappes, 
France  during  a single  day  in  June,  1926.  Three  of  his  releasee  led  to  sound- 
ings extending  above  13  km  find  In  these  three  Idrac  found  that  the  field  In- 
tensity, after  decreasing  In  the  characteristic  manner  through  the  tropos- 
phere, started  Increasing  above  the  tropopause.  The  average  field  strength  at 
8 km  for  all  of  his  flights  for  that  day  wus  only  2.3  v/m,  while  In  the  region 
above  13  km  It  reached  values  as  high  as  1»0  v/m.  Assuming  a uniform  vertical 
current  density  for  all  heights  reached  by  Idrac 's  balloons,  one  finds  that 
these  Intensity  values  Imply  a conductivity  decrease  by  a factor  of  almost 
twenty  In  going  from  8 km  up  to  "his  level  of  maximum  local  field  strength 
above  13  km.  Fortunately,  one  of  Idrac's  flights  extended  to  20  km  and  re- 
vealed that  the  neia  intensity  decreased  again  above  16  km,  falling  to 
1.2  v/m  at  19  km,  thus  shoving  that  the  conductivity  did  not  remain  low  that 
day  throughout  the  vertical  extent  of  the  stratosphere  over  France. 

These  observations  by  Idrac  seem  not  to  have  been  known  to  Gish  and  Sher- 
man, with  the  result  that  there  has  probably  been  less  significance  attached 
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to  the  gyp  lore  r H results  In  the  19-22  km  region  than  might  otherwise  have 
been  the  case.  Idrac's  findings  seem  particularly  valuable  In  that  they  re- 
fer to  a region  remote  from  that  In  which  the  Explorer  H flight  was  made 
and  to  a time  several  years  earlier,  thus  casting  doubt  on  any  supposition 
that  the  sort  of  low  conductivity  stratum  found  by  Gish  and  Sherman  was  an 
anomaly  peculiar  to  that  sounding.  Furthermore,  Idrac's  measurements  suggest 
that  there  may  be  only  a rather  thin  layer  of  low  conductivity  In  the  strat- 
osphere, while  the  Explorer  II  observations  left  this  Important  point  Inde- 
terminate. This  latter  contribution  of  Idrac's  work  strengthens  the  position 
taken  by  Holier  and  Saxon  (1952)  with  respect  to  the  slight  Importance  of  any 
regions  of  low  conductivity  In  altering  the  upward  flow  of  positive  current 
from  thunderstorm  to  Ionosphere,  while  the  former  contribution  (indication  of 
worldwide  extent  of  the  low  conductivity  layer)  points  to  the  need  for  further 
study  of  the  suggestion  made  by  Holier  and  Saxon  that  such  a layer  may  strong- 
ly Influence  surface  field -strength  fluctuations  far  from  active  thunder- 
storms and  squall  lines. 

5.  Possible  causes  of  the  conductivity  minimum 

Despite  the  lack  of  agreement  between  Idrac's  and  Gish  and  Sherman's  ob- 
servations of  the  altitude  of  the  base  of  the  region  of  low  conductivity, 
their  agreement  as  to  the  presence  of  such  a region  in  the  lower  stratosphere 
would  seem  to  Justify  some  attempt  to  find  an  explanation  for  Its  existence. 
Atmospheric  conductivity  Is  almost  entirely  controlled  by  the  small-lon  dens- 
ity of  the  air,  and  this  density  Is  In  turn  controlled  Jointly  by  the  rate  of 
Ion  formation  (by  cosmic  ray  Ionizations  followed  by  molecular  attachment) 
and  by  the  rate  of  destruction  (by  recombination  processes  and  by  attachment 
to  Altken  nuclei,  forming  relatively  immobile  large  Ions);  hence  one  must 
search  for  some  phenomenon  capable  of  locally  altering  one  or  both  of  these 
rates  In  the  stratum  under  consideration.  There  would  appear  to  be  no  basis 
for  believing  that  there  might  be  any  local  decrease  of  cosmic  ray  bombardment 
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here,  nor  any  local  anomaly  in  the  molecular  attachment  rates  or  in  the  re- 
cotnbl native  processes  for  small  ions,  so  one  is  led  to  seek  an  explanation 
for  the  low  conductivity  in  terms  of  the  effect  of  some  local  concentration 
of  nuclei  in  the  lover  stratosphere.  Gish  (1939);  in  first  discussing  this 
type  of  explanation,  suggested  the  possible  role  of  nitrogen  pentoxide  which 
has  been  detected  spectroscopically  in  the  region  from  16  to  UO  km.  Gish 
(1951)  has  also  noted  that  ozone  was  found  to  be  unusually  abundant  near  the 
top  of  the  Explorer  II  sounding;  but  he  has  not  indicated  how  either  of  these 
substances  might  ever  appear  in  the  form  of  particles  large  enough  to  serve 
as  large-ion  nuclei. 

1*.  Convective  transport  of  nuclei 

The  writer  has  been  led  to  consider  quite  a different  hypothesis  which 
seems,  at  first  Inspection,  qualitatively  more  probable  than  those  advanced 
by  Gish,  This  hypothesis  would  account  for  the  presence  of  a stratum  of  Alt- 
ken  nuclei  in  the  lowest  portion  of  the  stratosphere  in  terms  of  injection  of 
nucleus-rich  air  into  the  base  of  the  stratosphere  by  thunderstorm  updrafts 
that  sweep  the  air  up  through  the  troposphere  from  the  lower  levels  of  higher 
nuclear  density  and  then  expel  this  air  into  the  stable  base  of  the  isothermal 
region.  The  nuclei  so  added  from  time  to  time  by  thunderstorms  around  the 
world  would  not  remain  forever  in  the  lower  stratosphere  but  would  slowly 
leave  this  region  by  virtue  of  the  Joint  action  of  fall-out  and  turbulent 
diffusion.  The  critical  test  of  the  hypothesis  thus  becomes  that  of  inquiring 
whether  the  processes  of  addition  and  removal  might  reasonably  be  expected  to 
come  to  balance  with  a steady -state  nuclei  count  at  the  base  of  the  strato- 
sphere sufficient  to  explain  the  sort  of  decreased  electrical  conductivity 
found  by  Idrac  and  by  Gish  and  Sherman, 

5.  Thunderstorm  heights 

First  it  may  be  noted  that  the  heights  to  which  thunderstorms  extend  are 
in  reasonable  agreement  with  the  convective  hypothesis,  at  least  for  the  Idrac 


observations  of  a low -cons activity  layer  with  its  base  near  13  km.  The  aver- 
age heights  of  the  tops  of  thunderstorms  observed  by  radar  during  the  Thunder- 
storm Project  (Byers  and  Braham,  19^9)  was  slightly  over  11  km,  and  40  per 
cent  of  all  observed  storms  built  up  to  13  km  or  above.  It  is,  of  course,  no 
mere  accident  that  the  base  of  the  stratosphere  coincides  closely  vith  the 
maximum  altitude  reached  by  thunderstorm  updrafts  since  the  stable  density 
distribution  above  the  tropopause  precludes  appreciable  growth  into  the  strat- 
osphere. 

The  location  of  the  low -conductivity  layer  in  the  Explorer  II  sounding 
is  higher  than  can  be  accounted  for  In  terms  of  average  thunderstorms  of  mid- 
dle-latitudes, and  particularly  so  for  the  November  date  on  which  the  flight 
was  made.  Examination  of  the  temperature  distribution  prevailing  during  the 
flight  (Broobacher,  1936)  reveals  that  a double  tropopause  existed  over 
South  Dakota  on  flight  day.  The  lower  inversion  began  at  11.5  ks  and  the 
upper  began  at  about  18.7  km.  This  plus  the  fact  that  the  winds  near  the  top 
of  the  flight  were  southwesterly  suggests  that  the  balloon  may  have  been  in 
air  that  had  recently  come  from  lower  latitudes  of  the  Pacific  area  where 
thunderstorm  convection  is  better  able  to  transport  air  to  heights  approach- 
ing those  at  which  the  Explorer  II  encountered  the  decrease  of  conductivity. 
However  this  necessity  of  an  appeal  to  s tropical  origin  of  the  nuclei  over 
South  Dakota  must  be  regarded  as  a weakness  of  the  thunderstorm  hypothesis 
and  if  further  measurements  of  stratospheric  conductivity  should  reveal  that 
the  average  level  of  the  minimum  agrees  more  closely  with  that  found  by  Gish 
and  Sherman  than  with  that  indicated  by  Idrac's  work,  one  could  not  even  con- 
sider the  conveotlve  explanation  here  proposed. 

6.  Rate  of  convective  transport  of  nuclei 

The  first  step  in  a quantitative  check  of  the  convective  hypothesis  con- 
sists in  estimating  the  average  worldwide  rate  of  thunderstorm  transport  of 


nuclei  up  to  the  base  of  the  stratosphere.  Using  data  on  the  average  vertical 
distribution  of  Altken  nuclei  based  on  twenty -eight  balloon  flights  in  the 
troposphere  (Land a berg,  1938),  and  combining  these  with  some  recent  estimates 
of  the  vertical  distribution  of  thunderstorm  inflow  rates  (Braham,  1952)*  one 
finds  that,  during  the  entire  lifetime  of  an  average  thunderstorm  cell  of  the 
mid  die -latitude  type  considered  by  Braham,  about  3 x 1020  nuclei  may  be  ex- 
pected to  enter  the  updraft.  Of  this  total,  almost  three-fourths  of  the  nu- 
clei are  found  to  enter  the  cell  in  the  0-1  km  layer.  At  greater  heights, 
where  the  mass  of  air  entrained  is  larger  than  in  the  surface  layer,  the  nuclei 
count  has  fallen  off  so  much  that  the  weighted  average  Influx  of  nuclei  is 
much  less  than  in  the  0-1  km  interval. 

Not  all  of  these  3 x 1020  nuclei  are  to  be  regarded  as  reaching  the  out- 
flow region  at  the  very  top  of  the  storm,  however.  A certain  number  will 
serve  as  condensation  nuclei  and  will  thus  be  largely  removed  by  the  precip- 
itation process;  but  this  number  is  so  small  compared  to  the  total  number  of 
Altken  nuclei  (most  of  which  are  too  small  to  be  activated  for  growth)  that 
it  may  be  ignored  here.  Second,  relative  motions  of  nuclei  and  cloud  drops 
will  remove  some  nuclei  by  accretion,  but  this  mechanism  will  also  be  Ignored 
here  on  the  ground  that  the  collection  efficiency  for  this  capture  process 
will  be  very  low  in  view  of  the  small  size  of  the  nuclei.  The  third  process, 
which  cannot  be  ignored,  is  that  of  horizontal  outflow  of  updraft  air  prior 
to  its  reaching  the  tropopause.  Again  Braham* s data  on  mass-exchange  in 
thunderstorms  provides  a basis  for  an  estimate.  Braham  (1952)  finds  that  of 
the  total  of  9.0  x 10*®  kg  of  air  entering  an  average  storm  throughout  its 
duration,  only  1.8  by  10*®  kg  flows  out  at  the  200  mb  level  (about  12  km). 

Thus,  one  may  regard  only  1.8/9. 0 of  the  total  of  3 * 102®  nuclei,  l.e,, 

19 

about  7 x 10  nuclei,  as  being  expelled  from  the  storm  top  into  the  lower 
layers  of  the  stratosphere.  It  must  be  admitted  that  present  ignorance  of 
the  details  of  the  kinematics  of  the  outflow  pattern  at  these  levels  leaves 
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doubt  as  to  whether  even  this  latter  number  of  nuclei  may  safely  be  assumed 
to  remain  at  the  outflow  level  rather  than  to  subside  back  Into  the  upper 
troposphere,  but  this  assumption  will  be  made  here. 

Taking  the  average  lifetime  of  a thunderstorm  cell  as  about  one  hcur 
(Byers  and  Bra  ham,  19^9),  and  using  Gish  and  Walt’s  (1950)  estimate  of  3 z 
105  storms  as  the  average  Instantaneous  rate  of  occurrence  of  thunderstorms 
over  the  entire  globe,  one  finds  an  average  rate  of  stratospheric  addition 
of  6 z 10  ' nuclei  per  second  for  the  whole  world.  Overlooking  the  fact  that 
this  rate  of  transport  must  certainly  decrease  rapidly  with  Increasing  lati- 
tude to  nearly  zero  values  In  both  polar  regions,  one  finds  that  for  the  en- 

l8  ? 

tire  area  of  the  earth,  5 2 10  cm  , the  average  rate  of  thunderstorm  In- 
jection of  nuclei  Into  the  base  of  the  stratosphere  may  be  of  the  order  of 
2 

10  nuclel/cm  sec. 

7.  Rate  of  fall-out  of  nuclei 

Having  estimated  the  rate  of  addition  of  nuclei,  the  nezt  step  Is  to  est- 
imate the  rate  of  removal  In  order  to  compare  these  two  rates  as  a test  of  the 
convective  hypothesis.  Gish  and  Sherman  (1936)  have  given  estimates  of  the 
density  of  nuclei  required  to  account  for  the  low  conductivity  In  the  19-22 
km  Interval  of  their  Explorer  II  measurements,  so  In  spite  of  the  fact  that 
this  Interval  lies  several  kilometers  above  the  level  to  which  one  may  ezpect 
thunderstorms  to  penetrate  In  middle  latitudes,  and  In  spite  of  the  author's 
warning  that  not  too  much  quantitative  significance  should  be  attached  to 
their  nuclei  estimates,  these  values  will  be  usee  here  as  the  only  available 
estimate  of  the  nuclear  densities  In  the  lower  stratosphere.  Certain  addi- 
tional deductions  can  be  made  from  Idrac’s  data,  and  note  will  be  taken  of 
these  later. 

Most  Altken  nuclei  are  less  than  about  2 z 10"-’  cm  In  diameter  (Junge, 
1951),  ao  Stokes’  lav  may  be  applied  with  reasonable  accuracy.  Assuming  a 
mean  density  of  2 gm/cm?  for  the  nuclear  substances,  the  particles  may  be  ex- 
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pected  to  bare  a maximum  fell  velocity  of  only  about  10“^  cm/sec,  l.e.,  less 
than  a meter  per  day  for  even  the  largest.  The  average  Instantaneous  rate  of 
fall-out  per  unit  horizontal  area  of  the  stratosphere  Is  then  given  by  the 
product  of  this  velocity  and  the  prevailing  nuclear  density,  which  we  way 
take  as  2 x 10^  nuclel/cnr*  (Gish  and  Sherman,  1936).  Thus  the  downward  grav- 
itational flux  of  nuclei  Is  only  about  2 nuclel/cm2sec.  This  Is  almost  an 
order  of  magnitude  less  than  the  estimated  rate  of  addition  of  nuclei  by  up- 
drafts, so  If  fall-out  were  the  sole  mechanism  capable  of  removing  the  nuclei 
from  the  stratosphere,  one  could  conclude  that  the  convective  hypothesis  was 
confirmed  by  the  above  estimates;  but  In  addition,  downward  turbulent  diffu- 
sion must  be  considered. 

8,  Turbulent  diffusion  of  nuclei 

The  rate  of  vertical  turbulent  diffusion  of  nuclei  in  a layer  depends  on 
the  eddy  diffusion  coefficient  D,  and  on  the  vertical  density  gradient  of  the 
nuclei,  dn/dz.  Lettau  (1951,  Fig.  2)  gives  1()5  cm  2/sec  for  the  diffusion  co- 
efficient at  the  15  km  level,  and  from  Gish  and  Sherman  (1936,  Fig.  7)  one 

_2 

finds  the  density  gradient  In  the  layer  of  low  conductivity  to  be  about  10 
nuclel/cni4.  Hence  for  this  combination  of  data  the  turbulent  flux  is 
D u (lcPcm^Bec”*)  (lO^cm"*4)  « 10^cm"2sec~^ 

downward.  This  rate  is  three  orders  of  magnitude  greater  than  the  rate  of 
fall-out,  so  It  appears  that  one  may  quite  safely  neglect  fall-out  as  compered 
with  diffusion.  But  more  pertinent  to  the  present  discussion  Is  the  fact  that 
the  estimated  rate  of  downward  diffusion  of  nuclei  Is  some  sixty  times  greater 
than  the  estimated  rate  of  addition  of  nuclei  by  thunderstorms  (15  nuclel/cm? 
sec.).  One  seems  forced  to  conclude  that  the  convective  transport  hypothesis 
Is  quantitatively  inadequate  for  accounting  for  a uniform,  worldwide  stratum 
of  high  enough  nuclear  density  to  fit  the  Explorer  II  conductivity  measure- 


ments 
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9*  Discuss  Ion 

Haring  found  this  negative  result  In  the  effort  to  check  the  corrective 
hypothesis.  It  Is  Interesting  to  note  that  If  one  Beeka  points  In  the  compu- 
tations where  some  nod  If lea t Ion  might  be  made,  these  appear  to  be  almost  en- 
tirely changes  that  only  strengthen  the  evidence  against  the  hypothesis. 

First,  the  data  on  nuclear  densities  In  the  troposphere  (Landsberg, 

1923)  were  obtained  from  balloon  flights  made  over  weLl  settled  areas  where 
industrial  pollution  tends  to  give  counts  unrepresentative ly  high  for  the 
world  as  a whole,  so  any  revisions  here  would  certainly  lower  the  estimated 
convective  transport  rate. 

Second,  the  assumption  that  all  of  the  air  diverging  from  the  thunder- 
storms at  the  200  mb  level  remained,  along  with  ltB  suspended  nuclei,  at  the 
level  of  outflow  cannot  be  defended  too  well.  A thunderstorm  that  builds  up 
to  the  trope pa use  probably  succeeds  In  locally  pushing  up  the  stable  overly- 
ing stratospheric  air  but  complete  Intermixing  of  the  outflow  with  BtratOB- 
pherlc  air,  as  assumed  above,  Is  a rather  unlikely  extreme.  If  any  correc- 
tions could  be  made  here  they  would  undoubtedly  lower  the  effective  rate  of 
convective  addition  of  nuclei. 

Third,  an  attempt  to  Incorporate  Idrac's  findings  Into  the  estimate  of 
the  downward  diffusion  of  nuclei  from  the  underside  of  the  stratum  of  low  con- 
ductivity yields  an  even  higher  rate  of  removal  than  was  found  above  from  the 
data  of  Gish  and  Sheraan.  Idrac  reports  a nearly  twenty-fold  Increase  of 
field  strength  between  8 km  and  about  14  km,  which  Implies  a roughly  equal 
factorial  decrease  of  conductivity  In  this  Interval.  Gish  and  Sherman  (1936, 
Fig.  4),  on  the  other  hand,  found  a conductivity  decrease  of  a factor  of  only 
two  In  the  interval  from  19  km  to  22  km.  ThlB  difference  in  Implied  nuclear 
gradients  Is  thus  seen  to  amount  to  about  a factor  of  five,  making  It  corres- 
pondingly more  unlikely  that  thunderstorm  updrafts  are  steadily  counterbal- 
ancing downward  diffusion  of  nuclei. 
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yourth,  one  might  choose  to  use  the  eddy  dlffualvlty  value  corresponding 
to  the  interval  in  which  the  conductivity  decreased  in  the  Explorer  II  sound- 
ing (19-22  fen)  rather  than  to  use  the  value  for  the  15  km  level  as  was  done 
above  In  an  effort  to  simulate  conditions  prevailing  just  above  the  tropo- 
pause.  Lettau  (1951)  gives  lO^cm^/sec  for  D at  this  level  (down  a hundred- 
fold from  D at  15  km),  so  combining  this  vlth  the  previously  considered  val- 
ue of  the  density  gradient  one  obtains  an  estimated  diffusion  rate  of  10  nu- 
clel/cm-sec . This  is  Just  the  estimated  transport  rate,  so  this  fourth  re- 
vision is  the  first  one  to  favor  the  convective  hypothesis.  But  since  this 
revision  requires  that  one  deal  vlth  altitudes  too  great  to  match  observed 
thunderstorm  heights  for  any  but  tropical  latitudes,  it  provides  no  real  sup- 
port for  the  convective  hypothesis  anyvay . It  is,  however,  interesting  to 
note  that  Lettau  (1951)  has  suggested  that  the  very  rapid  decrease  of  D vlth 
height  Just  above  15  km  must  tend  to  produce  vhat  he  terms  a "dust  horizon" 
at  this  level,  and  cites  some  light -scattering  observations  in  support  of  that 
contention.  This  theoretical  and  observational  evidence  for  some  Bort  of  zone 
of  accumulation  Just  above  the  tropopauee  does  give  further  support  to  the 
▼ lev  that  a layer  of  high  nuclear  density  and  hence  low  conductivity  is  the 
rule  rather  than  the  exception  in  the  lover  stratosphere  but,  dees  not  clari- 
fy its  ultimate  origin. 

In  all,  it  would  seem  to  have  been  shown  here  that  the  convective  hy- 
pothesis for  transport  of  nuclei  to  the  stratosphere,  though  qualitatively 
quite  plausible,  is  not  quantitatively  compatible  with  such  observations  of 
the  conductivity  minimum  as  exist  at  present. 

10.  Concluding  remarks 

That  Idrac's  observations  of  the  vertical  variation  of  the  electric 
field  intensity  should  constitute  the  only  check  on  the  decrease  of  conduct- 
ivity with  height  in  the  lever  ctrutcspherc  found  in  the  Explorer  II  flight, 
and  that  even  this  check  should  have  gone  so  long  unnoticed  seen  regrettable. 
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Present-day  balloon-sounding  techniques  should  be  readily  capable  of  provid- 
ing data  on  the  behavior  of  the  field  up  to  almost  30  km.  Bence  the  writer 
wishes  to  recommend  that  repetitions  of  Idrac's  measurements  be  carried  out 
at  enough  different  localities  and  times  to  determine  whether  a layer  of  lew 
conductivity  Is  In  fact  always  present  Just  above  the  tropopause,  and  if  so 
to  determine  at  what  heights  It  lies  antf  whether  It  la  uniform  or  patchy  in 
nature.  A series  of  such  soundings  might  clarify  many  of  the  questions 
raised  In  the  present  examination  of  the  convective  hypothesis  of  the  origin 
of  such  a layer.  Such  measurements  would  also  shed  light  on  the  Interesting 
suggestion  (Holier  and  Saxon,  1932)  that  a layer  of  minimum  conductivity 
could  be  responsible  for  exaggerated  surface  field  fluctuations  far  from 
active  thunderstorms. 
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m.  A Note  on  Erroneous  Cloud  Physical  Applications  of  Raoult's  Lav 

Abstract— It  is  pointed  out  briefly  that  an  error  has  appeared 
persistently  in  the  form  of  Raoult's  lav  used  in  analyses  of  the 
lowering  of  the  vapor  tension  of  cloud  drops  due  to  the  presence  of 
dissolved  solutes.  The  magnitude  of  this  error  is  shown  to  become 
as  great  as  191  per  cent  when  the  solute  is  sodium  chloride.  At- 
tention is  called  to  the  way  in  which  this  error  vas  incorporated 
into  an  analysis  of  condensational  growth  of  cloud  drops  by  Howell. 

1.  Introduction 

The  purpose  of  this  note  is  to  call  attention  to  an  error  vhlch  has  ap- 
peared repeatedly  in  statements  of  Raoult's  lav  in  the  meteorological  liter- 
ature and  vhlch  has  led  to  certain  inaccuracies  in  at  least  one  recent  paper 
on  cloud  physics. 

2.  Raoult's  lav 

When  m'  moles  of  a non -electrolyte  are  dissolved  in  m moles  of  water, 
the  relationship  between  the  vapor  tension  e'  of  the  resulting  solution  and 
the  vapor  tension  e of  pure  water  is  given  by  Raoult's  lav  as 

e*  - e - aJ , 

e m'  ♦ m (1) 

which  may  also  be  rewritten  as 

_e_i  « m , (2) 

e m'  + m 

Equation  (2)  states  that  the  ratio  of  the  vapor  tension  of  the  solution  to 
the  vapor  tension  of  the  pure  solvent  equals  the  mole -fraction  of  the  solvent 
present  in  the  solution.  In  sc-called  ideal  solutions,  the  lav  is  (by  defin- 
ition) exact  at  all  concentrations,  and  in  real  solutions  of  non-electrolytes 
it  holds  to  a good  degree  of  approximation  at  low  and  moderate  concentrations 
(Daniels,  1948). 

When  lowering  of  the  vapor  tension  of  a solution  is  due,  on  the  other 
hand,  to  an  electrolyte,  equations  (1)  or  (2)  no  longer  apply  because  of  dis- 
sociation of  the  solute.  There  Is  not  at  present  unanimous  agreement  as  to 
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the  exact  nature  of  the  solutes  present  in  hygroscopic  nuclei  in  the  atmos- 
phere, but  the  substances  most  frequently  suggested  as  Important  (sea  salts, 
nitrous  acid,  sulfuric  acid)  all  hare  in  common  the  property  of  dissociating 
Into  ions  in  solution  and  hence  oust  give  vapor  tension  reductions  different 
from  those  predicted  by  (1)  or  (2),  which  apply  only  to  non-electrolytes. 
Despite  this  fact,  Raoult's  lav  is  given  only  as  one  or  the  other  of  the 
above  equations  by  Haurvltt  (19^1) , Lowell  (19**5)>  and  again  quite  recently 
by  Neuberger  (1951).  (It  might  also  be  noted  that,  in  the  last  reference, 
the  symbols  corresponding  to  m and  m'  in  (1)  and  (2)  are  incorrectly  identi- 
fied vlth  the  masses  rather  than  vlth  the  numbers  of  moles  of  solute  and  sol- 
vent. ) The  same  form  of  Raoult's  lav  appears  in  the  recently  revised  Smith- 
sonian Meteorological  Tables  (List,  1951),  and  forms  the  Incorrect  basis  for 
all  of  the  equilibrium  supersaturetlons  over  solution  droplets  tabulated  on 
pp.  375-579  of  those  Tables. 

Tot  electrolytes,  Raoult's  lav  must  be  modified  to  the  form 

e'  - e - - lm'  , 

e im'  + m (3) 

vhere  1 is  a factor,  often  called  the  van't  Hoff  factor,  vhlch  varies  both 
vlth  the  chemical  nature  of  the  electrolyte  and  vith  the  concentration  of  the 
solution.  In  the  limit  of  infinitesimal  concentration  of  the  solute,  1 be- 
comes simply  the  number  of  lens  comprising  one  molecule  of  the  solute  (e.g., 
tvo  for  RaCl,  three  for  MgClg). 

3.  The  van't  Hoff  1 factor 

As  the  solute  concentration  Increases  from  rero,  the  value  of  the  van’t 
Hoff  factor  first  decreases,  but  then  begins  to  rise  again  for  moderate  con- 
centrations and,  for  most  strong  electrolytes,  attains  values  in  excess  of 
the  number  of  ions  per  molecule  at  concentrations  near  the  saturation  value. 
The  latter  effect  13  particularly  marked  In  the  cases  of  salts  which  charac- 
teristically form  hydrates  (e.g.,  MgCig.SH^C).  The  Debye-Euckel  theory  of 


/ 
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lnter  Ionic  attraction  gives  a fairly  good  explanation  of  the  behavior  of  1 
near  zero  concentration  vhere  1 decreases  with  Increasing  concentration 
(Daniels,  19^8) . At  high  concentrations  It  Is  believed  that  the  clustering 
of  the  highly  polar  vater  molecules  about  the  Ions  (partlculerly  about  the 
cations)  ties  up  a large  fraction  of  the  total  population  of  vater  molecules 
to  decrease  appreciably  the  number  of  vater  molecules  escaping  per  second 
per  unit  area  of  surface  of  the  solution,  the  latter  thereby  behaving  as  If 
there  vere  an  apparent  Ionic  concentration  greater  than  the  actual  value. 

Since  the  variation  of  1 vlth  concentration  Is  dependent  upon  the  chem- 
ical nature  of  the  solute  In  question,  and  since  there  Is  not  yet  complete 
agreement  as  to  the  nuclear  substances  operative  under  natural  conditions,  It 
vouid  scarcely  be  In  order. here  to  undertake  an  exhaustive  examination  of  the 
physical  chemistry  of  any  one  nuclear  substance  or  mixture.  However,  In  or- 
der to  show  how  Raoult’s  lav  should  be  treated  for  vhatever  solutes  ultimately 
prove  to  be  of  chief  Importance  In  atmospheric  condensation,  the  case  of  HaCl 
nuclei  Is  considered  here  In  somewhat  more  detail  than  has  been  done  In  the 
meteorological  literature  before.  If  sea  salt  nuclei  should  prove  to  be  the 
main  atmospheric  nuclei,  then  since  NaCl  comprises  some  77  per  cent  oy  velght 
of  the  mixture  of  salts  present  In  sea  vater,  the  Raoult  effect  of  that  salt 
vill  be  of  primary  interest,  though  MgClg,  present  to  the  extent  of  about  ll 
per  cent,  vlll  have  to  be  considered  carefully  because  Its  property  of  form- 
ing a hexahydrate  vlll  make  It  a quite  significant  factor  In  the  physical 
chemistry  of  sea-salt  nuclei  under  conditions  of  low  relative  humidity. 

It  Is  possible  to  compute  the  van't  Hoff  factor  as  a function  of  molal- 
ity for  any  given  solute  from  the  variation  of  the  activity  coefficient  of 
that  solute  vlth  molality  (Moore,  1950).  The  great  vealtb  of  data  on  activ- 
ity coefficients  vouid  Justify  this  approach  ones  there  Is  clear  evidence 
that  some  one  substance  Is  of  dominant  Importance  in  atmospheric  nucieatlon, 
but  the  method  Is  quite  tedious,  Involving  as  it  docs  a numerical  Integration 
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of  the  Gibbs -Duhem  equation  for  each  determination  of  1.  As  a more  direct 
approach  available  In  the  case  of  NaCl  whose  effect  on  7apor  tension  has 
been  determined  experimentally  ever  a wide  range  of  concentrations,  i -values 
have  been  computed  here  from  the  observed  magnitudes  of  the  reduced  vapor 
tensions  of  solutions  of  NaCl.  Tabulated  values  of  this  quantity  are  given 
by  Washburn  (1926)  for  a range  of  molality  of  0.1  to  6.0  (a  saturated  NaCl 
solution  is  about  6.05  molal  at  atmospheric  temperatures).  Three  additional 

i 

values  needed  for  the  very  low  concentrations  attained  In  the  later  phases 
of  drop  growth  were  taken  from  the  measurements  of  Dleterlcl  and  of  Bmlts 

| 

reported  by  Roth  and  Scheel  (1923).  i 

Table  1.  Values  of  the  van’t  Hoff  factor,  1,  for 

Aqueous  solutions  of  KaCl  of  molality  M. 

Hole -fraction 

M of  NaCl  1 


0.O+* 

0.0+ 

2.00- 

0,044 

0.00080 

1.96 

! 

1 

1 

0.070 

0.00126 

1.90 

1 

0.098 

0.001Y7 

1.86 

1 

0.1 

O.OOlo 

1.83 

4 

i 

ft 

i 

0.2 

0. 0036 

1.82 

? 

1 

3 

• • 0.4 

0.0072 

1.84 

1 

i 

0.6 

0.0108 

1.85 

\ 

0.8 

0.0144 

1.87 

1 

i 

1.0 

0.0180 

1.89 

2.0  , 

0.0361 

2.04 

1 

i 

2.8 

0.0505 

2.19 

• 

i 

1 

5.0 

0.0902 

2.66 

i 

1 

* 

6.0 

0.1083 

2.Q1 

* 0.0+  is  used  here  to  denote  the  limit  of  Infinite  dilution. 

Table  1 shows  the  values  of  i computed  from  these  observed  data.  It 
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la  to  be  noted  that  even  at  the  minimum  point,  1 equals  1,65,  so  the  reduc- 
tion of  the  vapor  tension  of  a droplet  containing  a nucleus  of  NaCl  Is  at 
least  65  per  cent  greater  than  the  value  Implied  by  the  erroneous  form  of 
Raoult’s  lav  given  by  Haurvltz  (19^1).  Lovell  (19**5)>  Neuberger  (1951)»  and 
List  (1951).  Note  also  that  near  saturation,  1 has  become  so  large  (due  to 
solvation  of  the  Iona)  that  an  error  of  191  per  cent  Is  made  If  these  Ionic 
effects  are  neglected  In  calculating  the  Raoult  effect.  For  the  frequently 
accepted  case  of  sea  salt  nuclei,  the  corresponding  error  vould  be  notice- 
ably larger  at  this  high  concentration  because  of  the  much  more  pronounced 
hydration  of  the  magnesium  lens. 

4.  Kohler  curves  for  NaCl  nuclei 

That  the  forms  of  Raoult ' s lav  appearing  In  the  literature  cited  are  In 
need  of  correction  Is  Indicated  by  the  fact  that  this  error  has  been  carried 
Into  the  very  valuable  vork  of  Howell  (19^9)  on  cloud  drop  growth.  Curi- 
ously, the  neglect  of  Ionic  effects  seems  to  be  an  error  appearing  In  the  me- 
terorologlcal  literature  of  Just  this  country.  Kohler  (1936)  In  his  vork  on 
the  vapor  tension  of  droplets  assumed  a constant  ' 'd lsaoclatlon  factor"  (as 
defined  for  the  now  abandoned  Arrhenius  theory)  of  0.75  for  his  assumed  NaCl 
nuclei.  Wright  (1936)  considered  several  chemical  possibilities  for  his  nu- 
clei and  used  constant  factors  to  correct  for  the  effects  of  Ionization, 

Hla  statement  that,  "The  hygroscopic  factor... may  vary  vlth  the  concentration 
of  the  solution,  particularly  when  the  solution  becomes  very  dilute  and  the 
phenomenon  of  dleaoclatlon  occurs",  raises  some  doubt,  however,  as  to  vhether 
he  held  a correct  view  of  the  underlying  physical  chemistry  of  Raoult'o  lav. 
Beet  (195l)>  in  a recent  paper  on  drop  growth,  merely  employs  Wright’s  con- 
stant correction  for  solute  dissociation. 

Using  the  values  of  the  van't  Hoff  factor  given  In  Table  1,  (l.e.,  ef- 
fectively using  the  experimental  values  from  which  Table  1 was  derived),  the 
equilibrium  supersaturation  over  droplets  containing  various  specified  numbers 
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of  moles  of  NaCl  were  computed  for  comparison  with  Howell's  values  (present- 
ed by  him  as  his  Figure  1).  Each  Kohler  curve  corresponding  to  a siren  nu- 
clear site  was  pushed  down  as  a result  of  the  correction  of  the  Raoult's- 
law  error  made  by  Howell,  and  its  peak  value  of  critical  supersaturation 
was  moved  towards  a larger  drop-site.  The  corrected  values  of  peak  super- 
saturation  are  only  about  three -fourths  as  large  as  those  found  by  Howell 
(not  to  be  confused,  however,  with  the  supersaturations  found  by  Howell  in 
the  course  of  his  numerical  Integrations  of  the  drop-growth  equation,  for 
comments  on  which  see  below). 

The  Important  question  arises:  What  effects  on  Howell’s  computed  growth 
rates  would  follow  from  a corrected  treatment  of  Raoult's  law?  In  an  attempt 
to  answer  this,  the  writer  calculated,  from  Howell’s  paper,  the  parameters 
needed  to  evaluate  the  growth  rates  for  several  different  nuclear  slues  and 
then  sought  to  compare  these  rates  with  those  obtained  with  the  Incorrect 
form  of  Raoult's  law.1 

The  degree  of  supersaturation,  3 In  Howell's  nomenclature,  Is  a quantity 
which  must  be  known  before  a calculation  of  growth  rate  can  be  made  In  any 
given  case.  By  reading  off  the  S-values  from  Howell's  growth  curves  for  var- 
ious drop  sizes  and  nuclear  masses  of  Interest,  and  using  these  to  calculate 
corrected  growth  rates,  apparent  errors  of  from  25  per  cent  to  50  per  cent  In 
rates  (at  the  points  of  maximum  error)  were  found.  However,  this  was  a fal- 
lacious approach  which  overlooked  the  principal  contribution  of  Howell's  en- 
tire analysis,  namely  the  treatment  of  S as  a dependent  variable  In  the  growth 
equation.  The  degree  of  supersaturation  Is  Itself  strongly  Influenced  by  the 
vapor  tension  of  the  drops  (In  reality,  and  also  In  Howell's  admirable  anal- 
yals  of  the  growth  problem),  so  to  use  Howell's  supe saturations  along  with 

1 It  may  be  worth  noting  here  for  the  benefit  of  readers  concerned  with  the 
theory  that  there  appears  to  be  an  error  of  algebraic  sign  In  Howell's  form- 
ulation of  the  compensated  diffusion  coefficient  as  well  as  a discrepant 
location  of  a factor  of  2 77  In  the  asymptotic  form  of  this  coefficient  for 
very  snail  drops. 
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the  corrected  vapor  tensions  is  not  permissible  since  this  vould  appreciably 
overestimate  the  error  in  the  growth  rates. 

There  seems  to  be  no  other  way  to  determine  exactly  how  far  off  Howell's 
growth  rates  may  be  than  to  redo  his  Integrations.  This  vould  entail  an 
amount  of  labor  which  vould  carry  this  note  considerably  beyond  its  intended 
scope,  so  no  final  assessment  of  the  Raoult's  lav  error  in  Howell's  work  is 
given  here.  However,  from  certain  numerical-physical  arguments,  the  writer 
is  led  to  suspect  that  in  the  regions  of  most  rapid  drop  growth,  Howell's 
rates  may  be  too  low  by  about  10-15  per  cent.  If  so,  the  later  history  of 
drop  growth  might  be  sensibly  affected,  since  it  is  exactly  at  this  stage  of 
cloud  formation  that  the  peak  supersaturation  is  determining  the  lower  limit- 
ing site  of  nuclei  (and  hence  the  total  number  of  nuclei)  to  be  activated 
for  rapid  growth.  For  the  region  of  grow th  beyond  a few  microns,  dilution 
has  proceeded  so  far  that  the  entire  Baoult  effect  vanishes;  hence  the  latter 
phases  of  Howell's  theoretical  growth  histories  are  not  here  in  question  ex- 
cept Inasmuch  as  they  may  be  affected  by  the  peak  superaaturation. 
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17.  The  Shape  and  Aerodynamics  of  Large  Raindrops 

Abstract --The  physical  factors  which  might  be  expected  to 
control  the  shape  of  large  raindrops  are  surface  tension,  hydro- 
static pressure,  external  aerodynamic  pressure,  electrostatic 
charge,  and  internal  circulation.  Each  of  these  is  examined 
quantitatively  and  the  conclusion  reached  that  only  the  first 
three  play  important  roles  in  producing  the  deformation  charac- 
teristic of  large  raindrops.  By  analyzing  an  actual  drop  photo- 
graph, the  distribution  of  aerodynamic  pressures  is  deduced  and 
is  shown  to  imply  that  separation  in  the  airflow  about  a rain- 
drop has  significant  effects  on  drop  shape  and  on  a number  of 
physical  processes  occurring  at  the  surfaces  of  falling  rain- 
drops. 

1.  Introduction 

It  has  been  known  for  over  half  a century  that  large  raindrops  do  not 
possess  the  streamlined  form  popularly  described  as  the  "teardrop"  shape. 
High-speed  photographs  (Flower,  1928;  Edgerton,  1S59;  Blanchard  1950)  reveal, 
instead,  that  a drop  falling  through  the  air  exhibits  a marked  flattening  on 
its  lower  surface  and  smoothly  rounded  curvature  rather  than  conical  taper 
on  its  upper  surface.  For  an  example,  see  Figure  1 here.  This  long-recog- 
nized peculiarity  of  large  drops  has  never  been  adequately  explained,  and  only 
very  few  attempts  to  elucidate  this  matter  have  even  been  undertaken. 

Although  J.  J.  Thomson  made,  in  1885;  some  observations  on  the  shape  of 
liquid  drops  moving  through  various  fluids,  the  first  serious  attempt  to  ex- 
amine the  meteorological  problem  of  the  shape  of  large  raindrops  appears  to 
have  been  made  by  Lenard  (1901*).  Using  a vertical  airstream  with  water  drops 
suspended  freely  therein,  Lenard  carried  out  a number  of  experiments  on  ter- 
minal velocities,  deformation,  and  breakup.  He  noted  that  a finite  time  , 
somewhat  greater  than  a tenth  of  a second,  was  required  for  a large  drop  to 
attain  its  equilibrium  degree  of  deformation  and  suggested  that  this  might  be 
due  to  centrifugal  distortion  sat  up  by  internal  circulations  which,  for  in- 
ertial reasons,  took  a measurable  amount  of  time  to  become  established  by  the 
surface  friction  of  the  air  rushing  past  the  drop.  To  the  preeent  writer's 


Figure  1.  High-speed  photograph  of  a large 
vater  drop  falling  at  terminal  velocity.  Dia- 
meter of  the  drop  calculated  as  a sphere  is 
6 ima,  velocity  of  fall  8.8  m sec"*,  height 
of  fall  12  m.  Photograph  by  Dr.  C.  Magono, 
Hokkaido  Imperial  Uhlverslty;  provided  through 
the  courteoy  of  Dr.  U.  Nakaya. 
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know  ledge,  no  extension  of  this  interesting  beginning  of  the  stud;  of  the  drop- 

• 

shape  problem  was  made  during  more  than  fort;  ;ears  following  Lenard's  work, 
flower  (1928),  in  a stud;  of  falling  speeds  had  found  it  necessary  to  appeal 
to  Lenard's  theory  of  centrifugal  distortion;  and  even  when  Lavs  (19UI)  car- 
ried out  his  very  extensile  measurements  of  the  terminal  velocity  of  water 
drops,  no  other  theory  of  drop  deformation  was  available  to  be  Invoked  to  ac~ 
count  for  the  distortion  of  shape,  and  this  despite  the  fact  the  reality  of 
Lenai-d's  postulated  circulations  had  never  been  demonstrated. 

Spllhaus  (19^8),  In  a short  paper  on  raindrop  shape  and  falling  speed  has 
made  the  only  other  contribution  to  this  problem  that  has  come  to  the  writer's 
attention.  Spllhaus  suggested  that  the  vertical  flattening  of  large  drops  Is 
due  to  the  combined  action  of  surface  tension  and  aerodynamic  pressures.  Due 
to  the  deficit  of  external  pressure  around  the  waist  of  a drop,  "the  drop 
must  deform  so  as  to  reduce  the  ratio  of  Its  area  of  cross  section  to  perim- 
eter In  the  vertical  plane"  In  order  to  give  the  surface  tension  an  opportun- 
ity to  equilibrate  the  aerodynamic  forces.  Spllhaus  say  not  have  been  aware 
of  Lenard's  earlier  work,  for  he  neither  mentions  It  specifically  nor  gives 
any  consideration  to  the  centrifugal  effects  which  Lenard  held  to  be  solely 
responsible  for  producing  drop  deformation.  Neither  Spllhaus  nor  Lenard  of- 
fered any  explanation  of  why  large  drops  are  not  symmetrical  about  horlrontal 
planes  through  their  centers,  and  Spllhaus  explicitly  omitted  this  asymmetry 
from  hla  theory  In  order  to  be  able  to  use  experimental  data  on  drag  coeffi- 
cients of  oblate  ellipsoids.  As  will  be  pointed  out  belot1,  Spllhaus  used  an 
Incorrect  relationship  for  determining  the  surface  pressure  increment  due  to 
surface  tension  and  elso  erred  In  Ignoring  the  significant  effects  of  Internal 
hydrostatic  pressure  gradients  present  In  a drop  falling  at  terminal  velocity. 
With  this  background  to  the  present  problem,  It  seems  appropriate  to  conclude 
that  the  Issues  Involved  are  far  from  settled.  In  the  present  paper,  some 
further  contributions  to  this  problem  will  be  made  and  the  resulte  used  to 
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gain  an  improved  understanding  of  some  important  features  of  the  airflow 
about  large  raindrops. 

2.  Factors  controlling  raindrop  shape 

In  the  course  of  the  present  stud;  an  effort  has  been  made  not  only  to 
gain  a better  appreciation  of  the  role  of  centrifugal  distortion,  surface 
tension,  and  aerodynamic  pressures,  but  also  to  obtain  a clearer  recognition 
of  the  possible  importance  of  electrostatic  charges  and  internal  hydrostatic 
pressure  gradients. 

The  only  reason  that  water  drops  can  exist  at  all  as  mechanically  stable 
systems  is  that  surface  forces  at  the  water -air  interface  continually  try  to 
minimize  the  interfacial  energy  by  tending  to  minimize  the  interfacial  area. 
When  this  effect  of  surface  tension  acts  alone,  or  nearly  so,  as  in  the  case 
of  cloud,  drizzle,  and  even  small  raindrops,  it  succeeds  in  molding  a drop 
into  the  shape  characterized  by  minimum  surface -to-volume  ratio,  l.e,,  a 
sphere.  When,  however,  other  factors  than  surface  energy  contribute  signifi- 
cantly to  the  total  energy  of  the  drop,  minimum  total  energy  may  be,  and  in 
fact  is,  inconsistent  with  perfectly  spherical  shape.  One  might  hope  to  as- 
semble all  of  these  other  energy  factors,  express  mathematically  their  contrl 
buttons  to  the  total  drop  energy  and  then  determine  the  equilibrium  shape  by 
minimizing  the  total  energy  with  respect  to  some  suitable  shape  parameter  or 
parameters.  If  the  gravitational  effects  (hydrostatic  pressures)  were  the 
only  additional  energy  factor,  thle  might  be  done  here  .Just  as  it  has  been 
done  (by  tedious  numerical  processes)  for  the  case  of  pendent  and  sessile 
drops  (Aaam  , 19**9).  However,  anyone  fully  cognizant  of  the  difficulty  of  in 
corporatlug  the  aerodynamic  factor  into  this  type  of  approach  will  understand 
that  the  raindrop  shape  problem  will  probably  never  yield  to  any  analysis 
which  treats  it  as  a classical  minimal  problem.  Certainly  this  is  eo  If  that 
analysis  is  to  be  carried  out  by  manual  rather  than  electronic-computational 
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means, 

Recognli Ing  this,  the  writer  has  sought  to  approach  the  problem  almul- 
tuneously  from  two  directions  In  order  to  converge  ultimately  upon  a result 
(deduced  aerodynamic  pressure  distribution)  whose  correctness  may  be  Judged 
tolerably  well  by  comparison  with  certain  experimental  results  In  the  fields 
of  fluid  dynamics  and  cloud  physics.  The  central  Idea  In  this  analysis  has 
been  to  evaluate  all  of  the  factors  controlling  the  pressure  distribution 
Ins lde  a large  drop  and  then,  using  certain  surface  physical  concepts,  to 
determine  the  surface  pressure  prevailing  In  the  boundary  layer  Just  outside 
the  drop  surface.  If  the  surface  pressure  pattern  thus  deduced  Is  found  to 
be  In  reasonably  good  agreement  with  aerodynamic  principles  (as  will  be  shown 
to  be  the  case),  then  some  confidence  may  be  placed  In  the  theory  of  drop 
shape  on  which  the  calculations  have  been  based.  The  logic  of  this  approach 
will  be  further  elaborated  below. 

Surface  tension.  Surface  tension  holds  a rain  drop  together  In  the  face 
of  a number  of  tendencies  to  disperse  the  water  contained  In  the  drop.  A con- 
sequence of  the  net  Inward  attraction  exerted  on  a surface  molecule  by  the 
molecules  lying  deeper  within  the  drop,  this  surface  tension  also  produces 
an  Increase  of  pressure  within  the  drop  over  and  above  that  prevailing  In  the 
air  outside.  This  Increment  In  pressure,  Apfl#  at  a given  point  on  the  drop 
surface  Is  given,  In  general,  by 

APs  » ft'  U/Rj.  + l/Rs)  (1) 

where  Y'  Is  the  surface  tension  of  the  vater-air  Interface  and  R^  and  Rg  are 
the  principal  radii  of  surface  curvature  at  the  point  In  question  (Adam,  19**9)« 
The  quantity  A Pfl  can  be  either  positive  or  negative  If  one  admits  sufficient- 
ly arbitrary  surface  geometry , A principal  radius  will  here  be  regarded  as 
positive  for  the  case  where  tue  water-air  Interface  Is  convex  as  viewed  from 
the  air  and  APs  then  becomes  the  difference  between  the  water  pressure  Just 
Inside  the  drop  minus  the  (aerodynamlcally  controlled)  air  pressure  Just  out- 
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slde  the  Interface  at  the  point  in  question.  It  la  to  be  noted  that  with 
these  sign  conventions,  the  principal  radii  are  each  everywhere  positive 
(but  not  constant),  and  ^pfl  Is  everywhere  positive  (but  not  constant)  for 
a stable  drop  such  as  that  pictured  here  In  Figure  1.  When,  as  here,  both 
radii  are  of  the  same  sign  at  each  point,  the  surface  Is  said  to  be  every- 
where synclastlc.  At  a point  where  a general  surface  has  radii  of  opposite 
signs.  It  Is  there  antlclastlc.  During  the  processes  of  breakup  of  large, 
unstable  drops  the  surface  passes  from  the  wholly  synclastlc  over  Into  a 
partially  antlclastlc  form,  as  can  be  seen  In  photographs  of  artificially 
Induced  breakup  taken  by  Blanchard  (1950).  An  antlclastlc  surface  can  be- 
come dynamically  unstable  under  certain  conditions,  but  a synclastlc  surface 
cannot,  as  has  been  shown  by  Rayleigh  and  others  (see,  for  example,  Champion 
and  Davy,  1936). 1 

In  the  special  case  of  a spherical  drop,  R^  * % * T>  vllere  r 18  tbe 
drop  radius,  and  then 

Pfl  “ (2) 

r 

This  Is  the  equation  that  was  Incorrectly  applied  by  Spllbaus  to  his  assumed 
ellipsoidal  raindrop,  using  for  r the  radius  of  the  circular  cross  section 
In  a horizontal  symmetry  plane.  Since  thin  radius  Is  only  one  of  the  two 
principal  radii  of  curvature  at  a point  on  the  waist  of  such  a drop,  and 
since  the  second  principal  radius  Is  there  smaller  than  the  first,  Spllhaus 
underestimated  the  pressure  Increment,  particularly  for  his  very  large,  and 
hence  very  much  flattened  drops.  At  the  same  time  he  neglected  to  consider 
the  fact  that  the  surface  pressure  Increment  Is  different,  In  general,  at 
each  different  point  of  the  drop  surface,  so  bis  treatment  of  surface  ten- 
sion effects  was  doubly  Invalid  In  Its  details,  even  though  acceptable  in  a 
qualitative  sense. 

Since  ^ for  a water-air  Interface  et  OC  Is  75  dyne  cm"^,  (2)  reveals 
that  a spherical  drop  of  one  millimeter  radius  must  have,  at  that  temperature, 


Jp 

an  Internal  pressure  that  Is  about  1500  dyne  cm  a core  the  external  air  pres- 
sure. Tor  a drop  of  fire  millimeter  radius,  this  Increment  would  be  only  300 
dyne  cm"^  If  the  drop  could  somehow  remain  spherical;  and  at  the  other  extreme, 
for  e cloud  drop  of  five  micron  radius  the  internal  pressure  is  some  300  mil- 
libars above  the  external  air  pressure.  In  the  case  of  the  cloud  droplet,  the 
pressure  Increment  Is  so  very  large  compared  to  hydrostatic  pressure  differ- 
ences within  the  drop  and  to  the  minute  aerodynamic  pressures  established  at 
terminal  velocity  that  each  of  these  factors  (and  also  all  others)  may  safely 
be  neglected  in  discussing  drop  shape.  Hence  cloud  drops  do  simply  assume 
the  shape  Implying  minimum  surface  free  energy,  thus  accounting  for  their 
well  known  spherical  form.  But  in  the  case  of  a raindrop  at  the  upper  end  of 
the  observed  drop-size  distribution,  the  surface  pressure  Increments  are  only 
of  the  same  order  of  magnitude  as  the  pressure  effects  due  to  gravity  and  aero- 
dynamic factors,  so  for  this  case,  one  must  examine  the  shape  problem  more 
thoroughly. 

The  technique  for  determining,  in  general,  and  ft,  from  a photograph  of 
a falling  drop  will  be  explained  below  in  Section  3. 

Internal  hydrostatic  pressure.  As  a drop  falls  at  its  particular  terminal 
velocity  it  is,  by  definition,  no  longer  accelerating  in  the  gravitational 
field.  In  a coordinate  system  moving  with  that  falling  drop,  an  observer 
would  regard  the  drop  as  being  Just  supported  against  gravity  by  the  vertical 
components  of  the  aerodynamic  nornal  pressure  forces  and  the  surface  shear 
stresses  due  to  the  apparently  upward -rushing  air.  Consequently,  there  must 
exist  within  the  drop  a vertical  pressure  gradient  of  exactly  the  sort  found  in 
any  mass  of  fl»ld  at  rest  in  a gravitational  field.  This  hydrostatic  pressure 
gradient  appears  to  have  been  completely  overlooked  by  both  Lenard  and  Spllhaus, 
yet  la  the  limit  of  very  large  raindrops  the  difference  in  hydrostatic  prespur* 
between  top  and  bottom  of  a drop  becomes  quite  Important  In  controlling  drop 
shape.  Thus,  for  a drop  of  5 na  radius  (considered  by  Spllhaus,  thougn,  of 


course  atypical  of  natural  rain),  this  difference  In  hydrostatic  pressure  Is 
about  1000  dyne  cm"^,  or  about  three  times  larger  than  the  surface  pressure 
increment  for  a hypothetically  spherical  drop  of  this  same  radius.  In  the 
range  of  normal  drop  sizes,  the  hydrostatic  effect  Is  of  course  smaller,  while 
the  surface  pressure  contribution  becomes  larger:  A drop  of  1 mm  radius  has 

a top-to-bottoo  hydrostatic  pressure  difference  of  200  dyne  cm"^  as  compared 

O 

with  Its  1500  dyne  cm  surface  pressure  increment. 

If  one  could  show  that  there  exist  appreciable  Internal  circulations  In* 
side  raindrops,  then  It  would  be  necessary  to  take  account  of  the  dynamic  pres- 
sure gradients  that  would  Inevitably  be  associated  with  these  circulations. 

This  point  will  be  considered  below.  Here  It  will  merely  be  noted,  as  a rather 
Interesting  point,  that  If  raindrop  shape  were  Influenced  only  by  surface  ten- 
sion and  hydrostatic  effects  (with  uniform  external  pressure),  then  the  equi- 
librium shape  would  be  one  that  was  flattened  on  top  and  smoothly  rounded 
below,  l.e.,  Just  the  reverse  of  the  relative  curvature  observed  In  actual 
raindrops.  This  conclusion  follows  from  the  fact  that  the  drop  could 
than  only  be  In  equilibrium  (internal  pressures  in  hydrostatic  balance)  If  the 
surface  curvature  were  largest  near  the  base  and  very  small  near  the  top,  as 
suggested  In  Figure  2.  In  several  photographs  presented  by  Spells  (1952), 
liquid  drops  falling  very  slowly  through  less  dense  liquids  may  be  seen  to 
possess  exactly  this  sort  of  meridional  profiles.  The  reason  for  this  Is,  as 

will  be  shown  later  here,  that  at  the  low  Reynolds  numbers  at  which  Spells' 
drops  were  falling,  separation  does  not  occur  in  the  boundary  layer  and  hence 
no  dynamically  low  pressure  can  develop  over  the  upper  surface.  The  conclusion 
drawn  here  as  to  the  shape  a drop  would  possess  if  only  surface  tension  and 
gravitational  effects  were  involved,  Is  sufficiently  contradictory  to  the  ob- 
served shape  of  large  drops  to  Imply  clearly  that  other  physical  factors  must 
play  an  important  part  In  the  morphology  of  raindrops. 
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Flgure  2.  Equilibrium  shape  of  a falling  liquid  drop  If  only 
surface  tension  and  hydrostatic  pressure  effects  were  significant. 
Observed  shape  of  large  raindrops  Is  approximately  that  ob- 
talned  by  turning  above  shape  upside  down. 

The  vay  In  which  one  may  consider  hydrostatic  effects  quantitatively  In 
analyzing  drop  shape  vlll  be  discussed  below  In  Section  ?. 

Electrostatic  charges.  Since  it  Is  known  that  hydrometeors  of  all  sizes 
ranging  from  cloud  droplets  up  to  the  largest  raindrops  may  carry  electric 
charge , It  Is  necessary  to  consider  the  possibility  that  the  drop-shape  prob- 
lem might  be  sensibly  affected  by  this  factor. 

By  electrostatic  standards  the  water  in  a natural  raindrop  Is  a good 
conductor.  It  can  be  shown  (e.g. , Jeans,  19**1,  p.  79)  that  a conductor  car- 
rying a local  surface  charge  density  (j" experiences  an  outward -d lrected  ten- 
sion (negative  pressure)  whose  magnitude  per  unit  area  Is  given  by 

T = 27TCT2.  (5) 

This  electrostatic  tension  opposes  the  surfuco  tcncicn  thus  constitutes 
one  of  the  several  destabilizing  factors  that  control  drop  morphology.  This 
point  assumes  real  Interest  as  soon  as  one  notes  that  on  any  conductor  of 
verlable  surface  curvature  there  Is  a tendency  (counteracted  only  by  external 
fields  due  to  neighboring  charged  bodies)  for  the  charge  to  distribute  Itself 
In  such  a manner  that  O'  becomes  largest  where  the  surface  curvature  is 
largest.  Hence  aa  a charged  raindrop  began  to  flatten  out  due  to  aerodynamic 
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effects,  more  of  the  total  available  charge  would  migrate  towards  the  waist 
of  the  drop  and  would  there  produce  a locally  exaggerated  suppression  of  the 
surface  tension  effect  which  serves  to  oppose  the  effect  of  low  aerodynamic 
pressures  near  the  waist.  Consequently  the  drop  would  have  to  deform  still 
more  to  increase  the  quantity  (l/R^  ♦ l/Rg)  around  the  waist  in  an  effort  to 
attain  pressure  equilibrium.  This  further  flattening  would  then  not  only 
further  decrease  the  external  air  pressure  at  the  waist,  due  to  continuity- 
and  Bernoulli -effects  in  the  airflow,  but  would  at  the  sane  time  have  the  ad- 
ditionally unfavorable  effect  of  calling  for  still  further  buildup  of  surface 
charge  density  near  the  increasingly  sharply  curving  waist.  This  would,  in 
turn,  oppose  even  more  strongly  the  surface  tension  effects  that  are  trying 
to  hold  the  drop  together,  and  so,  until  the  drop  became  so  flattened  as  to 
be  torn  apart  by  aerodynamic  forces.  This  qualitative  picture  suggests  so 
vividly  and  plausibly  a mechanism  for  the  breakup  of  large  raindrops  in  thun- 
derstorm precipitation  currents  that  it  is  perhaps  regrettable  that  it  must 
next  be  shown  that  this  Interesting  hypothesis  is  quant  it  at.  ire  ly  tenable  only 
for  quite  a buOx  UJGL  1 degrees  of  charging  of  the  drops. 

In  the  region  Just  outside  a point  on  the  surface  of  a raindrop  having 
local  surface  charge  density  S , the  electric  field  intensity  is 

E - UiTtr  . (4) 

Now  the  greatest  possible  value  that  <5*  can  assume  is  given  by  (4)  when  E is 
set  equal  to  E<j,  the  dielectric  strength  of  the  surrounding  air.  Any  greater 
surface  charge  will  induce  corona  discharge  that  will  reduce  E to  Eg.  is 
pressure-dependent  and  is  also  sensitive  to  the  geometry  of  the  charged  con- 
ductors involved,  but  for  the  cloud  physical  problem  at  band,  one  will  be  con- 
servative (in  the  sense  of  admitting  rather  high  values  of  *y'  and  hence  of  T) 
if  he  pule  E^  = 39,000  volt  cm”^  in  (4),  solves  for  the  implied  surface  charge 
density,  and  then  inserts  this  into  (3)  to  determine  the  greatest  value  T can 
assume  before  corona  discharge  sets  in.  The  result  is 
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Tj  * !LlLElr>  * ^ s IL  * 400  dyne  cm"^ 

TmFT2  8 If  87T 

for  the  sea-level  value  Eg  * 50,000  volt  cm"1  * 100  e.s.u.  At  the  500  oh 
level,  where  Eg  falls  to  about  15,0CC  volt  em“*,  Tg  could  be  no  greater  than 
100  dyne  cm  . 

These  electrostatic  effects  (decrements  of  pressure  on  passing  from  the 
air  side  over  Into  the  vat er  side  of  the  drop  surface)  are  of  the  order  of 
magnitude  of  £^p8  for  large  drops  as  calculated  roughly  above,  but  they  have 
been  computed  here  for  values  of  CT  that  appear  to  be  substantially  larger 
than  any  yet  observed  directly  or  Indirectly.  Thus  a drop  of  5 radius 
vlth  so  high  a surface  charge  density  as  to  imply  E * Eg  at  Its  surface  vould 
bear  a total  charge  of  qg  « r^Eg  « 0.25x100  - 25  e.s.u.,  at  the  earth's  sur- 
face or  about  12  e.s.u.  at  the  500  nb  level.  Gunn  (1947  1950),  however, 

found  by  direct  measurement  from  aircraft  flying  through  precipitating  clouds 
that  drops  seldom  bear  charges  In  excess  of  0.1  e.w.u.  Furthermore,  his  val- 
ues wore  notable  in  that  they  are  almost  an  order  of  magnitude  greater  than 
those  previously  reported  for  raindrop  charges  as  measured  e.t  the  earth's  sur- 
face (Chalmers,  1949).  Gunn  (1949)  has  given  an  interesting  possible  explana- 
tion of  vhy  raindrops  may  not  be  able  to  accumulate  charge  to  such  a degree  as 
to  produce  corona  discharge  of  the  type  here  tacitly  assumed  (isolated  drops 
discharging  Into  the  air). 

It  vould  seem,  then,  that  attractive  as  is  the  electrostatic  factor  In 
explaining  drop  shape  and  breakup,  such  observations  as  do  sxlst  emphasise  Its 
slight  quantitative  Importance  under  normal  circumstances.  The  electrostatic 
factor  should,  perhaps,  be  held  In  mind  as  possessing  possible  significance  In 
processes  occurring  In  the  regions  of  highest  electrical  activity  of  thunder- 
storms, but  elsevhere  It  may  be  Ignored  on  the  basis  of  existing  drop-electric- 
al measurements.  It  may  be  noted  that  in  the  regions  of  peak  field  strength 
within  thunderclouds,  the  shape  problem  must  also  be  considered  in  relation 


to  the  Micky  effect  due  to  external  electric  fields  (Macky,  1931) » but  this 
case  will  not  be  treated  In  the  present  paper. 

Internal  circulations.  Real  fluid  flow  la  fundamentally  distinguished 
from  perfect  fluid  flow  by  the  fact  that  the  former,  but  not  theslatter,  Is 
characterised  by  the  "condition  of  no  slip"  at  surfaces  bounding  the  region 
of  flow.  When  the  boundary  lp  a solid,  as  In  the  case  of  flow  orer  an  air- 
foil, this  condition  Implies  that  the  air  In  contact  vlth  the  boundary  Is  at 
rest  relative  to  that  boundary.  But  when  the  boundary  Is  the  surface  of  a 
liquid,  as  in  the  case  of  a raindrop  falling  through  the  air,  this  condition 
can  be  satisfied  even  If  the  surface  layer  of  air  Is  slowly  moving,  for  the 
lntferfaclal  liquid  may  be  drifting  downstream  at  some  slow  rate.  In  the 
case  of  a raindrop  moving  downward  through  air,  any  such  surface  circulation 
Induced  by  shear  stresses  exerted  by  the  ambient  air  would  In  turn  Induce 
some  sort  of  ax lsymmetr lc  Internal  circulation. 

Lenard  (1904)  postulated  the  existence  of  such  internal  motions  and  re- 
garded them  as  capable  of  accounting  for  the  drop  deformations  which  he  had 
observed;  but  he  made  no  attempt  to  demonstrate  their  reality  experimentally 
nor  to  predict  their  Intensity  theoretically.  Qualitatively,  one  can  say 
that  at  least  a very  slow  Internal  circulation  Is  almost  Inevitable;  for  the 
dynamic  boundary  condition  pertinent  here  Is  that  of  continuity  of  tangen- 
tial shear  stress  across  the  water-air  Interface,  and  since  water’s  viscosity 
Is  not  Infinite,  at  least  a slight  amount  of  Internal  motion  seems  certain 
to  develop.  In  the  analysis  that  will  be  given  below  In  8eetlon  3 It  will  be 
necessary  to  know  whether  the  pressure  at  a point  on  the  vertical  axis  of  a 
drop  la  equal  to  that  at  a point  at  the  same  height  above  the  base  of  the 
drop  but  lying  Just  Inside  the  drop  surface.  In  view  of  the  relatively  small 
radii  of  curvature  of  the  water-particle  trajectories  In  any  internal  circu- 
lations, one  can  show  that  horizontal  uniformity  of  preBBure  would  be  notice- 
ably altered  if  surface  water  velocities  of  much  over  one-tenth  the  drop’s 
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f ailing  speed  can  be  developed.  Consequently  It  becomes  indispensable  bare 
to  obtain  an  estimate  of  the  Intensity  of  the  internal  circulation. 

Bond  (1927)  has  examined  theoretically  the  problem  of  internal  circula- 
tions for  the  case  of  a liquid  sphere  moving  through  a dissimilar  liquid 

1 

medium  for  the  case  of  fcynolda  numbers  in  the  Stokes  lav  range,  and  has 

I 

shown  that  a dimensionless  quantity, 

k « 2 A*  + 

j/*  + 3>'/ 

is  a measure  of  the  extent  to  which  the  drag  law  for  a drop  of  liquid  of 

j 

dynamic  viscosity  p moving  through  a medium  of  viscosity  /*  departs  from  the 
drag  law  for  a rigid  sphere  moving  through  the  same  medium.  Since  for  air 

j.  / j. 

is  about  1.7X10  poise  at  OC,  while  /A  for  water  Is  about  170x10  poise  at 
the  same  temperature,  k Is  about  302/303  * 0. 997  for  the  raindrop  case,  where- 
as It  would  be  1.0C0  for  the  case  of  a falling  solid  sphere.  Eence  for  drop 
sites  whose  terminal  velocities  lie  within  the  Stokes  law  range,  It  appears 
that  internal  circulations  must  be  relatively  Insignificant.  However,  for  the 
sort  of  drops  of  primary  interest  here,  the  Reynolds  numbers  lie  in  the  range 
from  about  1()2  to  10^  so  Bond’s  criterion  Is  not  rigorously  applicable  to  the 
present  problem.  Nonetheless,  this  part  of  Bond's  analysis  Is  of  qualitative 
value  here  in  that  it  calls  attention  to  the  fact  that  the  development  of  in- 
ternal circulations  does  hot  depend  on  Just  the  external  Reynolds  number,  but 
rather  on  the  relative  viscosity  of  the  Interior  hnd  the  eiterior  fluids. 

It  Is  particularly  necessary  to  keep  this  latter  Implication  of  Bond's 
work  clearly  in  mind  In  examining  a number  of  experimental  studies  recently 
carried  out  to  determine,  for  chemical  engineering  reasons,  the  nature  of  the 
Internal  circulations  occurring  Inside  a heavy  liquid  drop  falling  through  a 
lighter  liquid  (Garner,  1950;  Spells,  1952).  Circulations  were  readily  ob- 
servable photographically  in  these  stuuloo,  but  they  all  concerned  situations 
■«here In  there  was  only  a very  alight  difference  in  viscosity  between  the  drop 
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material  and  the  asternal  medium,  while  in  the  cnee  of  a raindrop  falling 
through  air  the  viscosity  of  the  drop  exceeds  by  two  orders  of  magnitude  that 
of  the  surrounding  medium.  Hence,  one  must  note  veil  that  these  recent  studies 
do  not  bear  very  close  relation  to  the  raindrop  problem. 

Richardson  (1950)  has  made  some  experimental  studies  of  the  breakup  of 
vater  drops  falling  from  a tover  125  feet  high  and  has  suggested  that  their 
breakup  is  due  to  the  effects  of  internal  circulation;  but  he  appears  to  have 
obtained  no  direct  evidence  to  support  this  viev.  He  did  observe  that  drops 
of  a very  viscous  liquid  (thickened  methyl  salicylate)  resisted  breakup  far 
more  effectively  than  did  vater  drops,  and  argued  that  this  vas  due  to  their 
resistance  to  the  development  of  internal  circulations.  On*  must  question 
this  interpretation  on  the  ground  that  the  increased  viscosity  can,  and  in- 
deed must,  also  play  an  important  role  in  suppressing  breakup  by  inhibiting 
the  rapid  oscillations  that  probably  initiate  breakup  in  oversize  drops 
(Blanchard,  1950). 

Blanchard  (1949)  attempted  to  observe  internal  circulations  In  vater  drops 
suspended  in  a vertical  vlnd  tunnel.  He  introduced  fine  particles  of  alumina 
into  the  drops  to  serve  as  tracers  in  revealing  Internal  motions,  but  reported 
no  evidence  for  any  circulations.  His  observations  may  be  Inconclusive  due  to 
the  fact  that  his  tracer  particles  might  have  been  too  large  ("300  microns  and 
under")  to  be  carried  along  in  currents  of  the  order  of  centimeters  per  second. 
Kinzer  (unpublished)  has  observed  very  slow  rolling  motions  inside  drops  con- 
taining fine  talc  particles,  but  he  estimated  the  velocities  involved  to  be 
less  than  a centimeter  per  second,  vhlch  can  be  shown  to  be  too  6lcv  to  have 
appreciable  centrifugal  effect  on  drop  shape. 

Because  of  the  importance  of  settling  the  question  of  whether  internal 
circulations  should  be  important  in  the  drop-shape  problem,  the  writer  has 
sought  a theoretical  basis  for  estimating  the  upper  limit  to  the  circulation 
velocities  cnat  could  develop  within  a large  drop  falling  at  the  ex  per  mentally 
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establlshed  limiting  speed  of  about  8 a »«c*^.  The  analysis  Is  based  upon 
the  following  assumptions j 

(1)  The  shear  stress  Is  continuous  across  the  water-air  Interface. 

(2)  The  condition  of  no  slip  holds,  but  the  Interfacial  molecules  of 
water  and  air  drift  together  so  slowly  downstream  that  as  far  as 
the  external  flow  dynamics  are  concerned,  the  surface  air  molecules' 
speed  may  be  aseumed  to  be  negligibly  small  compared  to  the  relative 
airspeed  at  Infinity  (8  m sec-^). 

(3)  The  drop  remains  spherical.  A corollary  to  this  assumption  Is  the 
Implication  that  the  stream  function  Inside  the  drop  corresponds  to 
that  of  a spherical  vortex. 

Any  exception  to  assumption  (1)  would  Imply  finite  shear  stress  acting  on  an 
Infinitesimal  lamina  of  fluid,  l.e.,  Infinite  accelerations  would  result 
from  any  failure  for  this  condition  to  hold,  so  this  firet  assumption  Is 
above  reproach.  That  the  surface  air  velocity  at  the  bese  of  the  boundary 
•layer  may  be  regarded  as  zero  compared  to  the  falling  speed  of  the  drop 
(assumption  (2))  Is  an  assumption  that  can  only  be  tested  a posteriori.  It 
will  be  shown  below  to  be  admissible.  Assumption  (3)  1b  Introduced  to  sim- 
plify the  analysis  even  though  It  falls  to  hold  In  the  case  of  Interest  here. 
However,  departure  from  sphericity  ought  not  have  any  large  effect  on  the 
chances  for  development  of  Internal  circulation  since  the  experimental  work 
of  Gamer  (1950)  and  Spells  (1952)  demonstrates  that  the  type  of  vortices 
predicted  by  Bond  (1927)  appear  In  highly  deformed  drops  of  liquid  In  liquid. 
Assumption  (3)  simplifies  the  analysis  primarily  because  of  the  corollary 
Implication  concerning  the  nature  of  the  internal  stream  function.  For  a 
spherical  drop  the  well  known  stream  function  for  a spherics!  vortex  may  be 
used  to  determine  the  motion, 

A stream  function  exists  for  simply  kinematics!  reasons  by  virtue 
of  the  continuity  equation,  and  so  holds  whether  the  motion  is  Irrotatlonal 


or,  as  here,  rotational.  The  general  fora  of  the  streaa  function  f for  any 
axleymmetrlc  flow,  when  specialized  to  the  case  of  motion  within  a sphere, 
reduces  to  the  form 


p = A(rz- a2)r"  si»*& 


where  A Is  a constant,  a Is  the  radius  of  the  sphere,  r Is  the  radial 
coordinate  measured  from  the  center  of  the  sphere,  and  6 Is  the  meridional 
coordinate  measured  from  a polar  axis  directed  towards  the  pole  of  the 
sphere  towards  which  fluid  moves  along  the  polar  axis,  see  Figure  3 and 


— Stagnation 
point 


Representative  external 
streamline 


Figure  3.  Definition  sketch  of  vortical  circulation 
Inside  a spherical  drop. 


The  constant  in  (3)  assumes  the  value  3U0/4a^  for  a spherical  vortex  In  a 
perfect  fluid,  where  UQ  is  the  translation  speed  of  the  vortex  center  rela- 
tive to  the  surrounding  fluid,  and  assumes  the  value^I^/ta^^'  +/<  ) for 
a spherical  vortex  composed  of  a fluid  of  dynamic  viscosity/^  moving  very 
slowly  through  a fluid  of  dynamic  viscosity  (Bond,  1927).  Due  to  the  par- 
ticular way  in  which  (5)  will  be  used  here,  It  will  be  unnecessary  to  specify 
A (which  is  comforting  Inasmuch  us  the  raindrop  case  corresponds  neither  to 


the  perfect  fluid  Unit  nor  to  the  low  Reynolds  numbers  treated  by  Bond) 


hence 


The  meridional  Telocity  u*  at  a general  p^lnt  r,  9 la  given  by 

ji'(n  e)~±  if\  ---  A (*  r*-24'^(6) 

' 'f,  6 


- 2 a2  A . 


(7) 


(8) 


(9) 


Also,  the  velocity  shear  Just  inside  the  surface  of  the  sphere  at  Its 
waist  is,  from  (6), 

:■?"'!  = h 

\ir  >*.•/* 

so,  combining  (7)  and  (8)  to  eliminate  A,  It  follows  that 

y (3'  wfz)  ~ -v"(  Tr)a,*/£ 

and  the  problem  of  determining  u’  (a,  7 T /2)  becomes  that  of  finding  ^u’/<^/^ 
at  the  same  point.  This  can  next  be  done  with  the  aid  of  assumption  (1). 
Continuity  of  shear  stress  across  the  water-ai'-  interface  requires 

where  the  primed  quantities  refer,  as  before , to  the  water  and  the  unprimed 

t O 

to  air.  Since  /A  £ 1C  , one  has 


-2 


( r)  t # 


(10) 


(sL"/'  ^ to 

x 4 r /A 

The  velocity  shear  In  the  boundary  layer  of  air  Just  outside  the  drop  Is  of 
the  order  of  u^/ £ where  u^  is  the  local  air  speed  relative  to  the  sphere 
at  the  outer  limit  of  the  laminar  boundary  layer  and  ^ Is  that  layer's 
radial  thickness.  Since  outside  of  the  boundary  layer  the  flew  at  raindrop 
Reyn0j.u9  numbers  will  be  essentially  potential  flow  at  least  up  to  about 
the  waist  of  the  drop,  and  furthermore  since  to  assume  potential  flow  Is  to 
be  conservative  here  In  the  sense  of  admitting  rather  large  waist  velocities 
In  the  airflow,  ui  is  taken  as  1.5U0,  where  U0  In  the  rc lndrop  case  Is  the 
terminal  falling  velocity.  On  the  other  hand,  to  evaluate  £ , the  boundary 
layer  theory  of  Tomotlka  (1955)  will  be  utilized,  Tomotika  showed  that  for 
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a point  at  about  80°  from  the  forvard  stagnation  point  of  a sphere, 

6.1?  iT; 


£ 


2 Vo  V>3l 


(ii) 


vhere  "V*  Is  the  kinematic  viscosity  of  the  air  and  a la  the  radius  of  the 
sphere.  For  UQ  « S a sec  “*  and  a * 0.5  cm,  £ la  found  from  (6)  to  be 
about  0.05  cm.  Hence,  near  the  valst  of  such  a large  drop,  9u/3r  ± Uj/ £ 
1200/0.05  sec"*  ■ 2.UX1014  see"*,  and  then  from  (10)  Bu ' /^t  sec"*. 

Combining  this  last  result  vlth  (9),  one  finds  that  u*(a,  Jf /2)  *^30 


cm  sec"*,  It  Is  to  be  noted  lmned lately  that  since  this  speed  Is  less  than 
5 per  cent  of  the  air  speed  (12  a see”*)  at  the  outer  limit  of  the  boundary 
layer  near  the  valst  of  a drop  falling  at  the  maximum  speed  of  about  8 m 
sec'*,  assumption  (2)  is  rendered  quite  plausible  a posteriori.  If  a simi- 
lar calculation  Is  carried  out  for  the  more  probable  value  of  2.5  mm  for 
the  radius  of  a "large*  drop,  the  circulation  velocity  at  the  valst  Is  found 
to  be  only  about  20  cm  sec"*. 

The  circulation  Intensity  predicted  by  the  above  argument  Is  substantial- 
ly higher  than  that  observed  by  Klnzer  (In  smaller  drops),  and  It  vlll  become 
possible  later  In  this  paper  to  point  out  a very  good  reason  vhy  the  actual 
circulations  fall  to  reach  the  Intensity  Just  predicted.  Hovever,  since  this 
point  vlll  depend  upon  a deduction  vhleh  hinges  In  part  upon  the  negligibility 
of  Internal  circulations,  It  Is  here  neces£«ry  to  proceed  to  shov  that  even  a 
surface  motion  of  50  cm  sec"*  In  a large  drop  vlll  not  lead  to  Internal  pres- 
sure gradients  large  enough  to  alter  seriously  the  hydrostatic  balance  Inside 
the  drop.  Consider  the  radially  lnvard  force  that  would  develop  as  a cen- 
trifugal reaction  to  the  here  predicted  vortical  motion.  The  radii  of  curv- 
ature of  the  trajectories  of  particles  moving  Just  Inside  the  boundary  of  a 
spherical  vortex  near  its  valst  are  nearly  identical  vlth  the  radius,  a,  of 
the  vortex  boundary,  as  may  be  seen  in  Figure  5*  Since  the  speed  of  circula- 
tion falls  off  rapidly  Inward  ( XL*  oCe  ) as  one  approaches  the  ring  of 
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stagnation  points  lying  In  an  equatorial  circle  at  0.71a  from  the  axle,  and 
since  furthermore  the  centrifugal  reaction  varies  as  the  square  of  the  cir- 
culation velocity,  one  will  obtain  a reasonable  order  of  magnitude  estimate 
of  the  centripetal  force  Involved  If  he  assumes  that  a lamina  of  radial 
thickness  equal  to  about  one-tenth  of  the  distance  from  the  surface  to  the 
stagnation  point  moves  meridionally  near  the  valst  with  the  predicted  speed 
of  30  cm  sec"1  In  a very  large  drop  of  0.5  cm  radius.  The  centripetal 
force  per  unit  aree  acting  on  this  lamina  Is  found  to  be  about  30  dyne  cm"^. 
This  contribution  to  the  Internal  pressure  field  Is  only  10  per  cent  of  the 
surface  tension  contribution  of  about  300  dyne  cm for  a drop  of  this  site, 
and  Is  a still  smaller  fraction  of  the  hydrostatic  pressure  difference  from 
top  to  bottom  of  such  a drop.  For  a drop  of  0.25  cm  radius  one  finds  that 
the  centrifugal  pressure  effect  represents  a much  smaller  fraction  of  the 
surface  tension  Incremental  pressure.  Hence  It  is  concluded  here  that  one 
may  neglect  Internal  circulation  within  falling  raindrops  as  a first  approx- 
imation, at  least  as  far  as  such  circulations  might  affect  the  shape  of  drops 
In  the  range  of  sires  now  known  to  occur  In  natural  rain.  This  conclusion 
■will  be  strengthened  in  Section  U,  when  the  phenomenon  of  separation  Is  dis- 
cussed. 

Aerodynamic  pressure  distribution.  In  the  writer's  opinion,  all  of  the 
factors  capable  of  Influencing  drop  shape,  with  the  Important  exception  of 
the  aerodynamic  factor,  have  now  been  considered.  The  next  logical  step  for 
completing  the  drop-shape  theory  should  thus  be  a direct  evaluation  of  there 
aerodynamic  pressure  effects.  Unfortunately,  to  carry  out  this  last  step 
would  be  extremely  difficult.  Without  belaboring  this  fairly  obvious  point, 
it  may  be  noted  that  the  aerodynamic  pressure  distribution  over  the  surface 
of  a drop  falling  through  the  air  Is  Itself  determined  by  the  very  shape  one 
wishes  to  deduce.  One  could  only  proceed  here,  in  principle,  by  some  method 
of  successive  approximations  In  which  each  aerodynamic  pressure  calculation 
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( based  on  the  previous  Iterative  approximation  to  the  equilibrium  shape) 

\’ould  be  used  to  deduce  a modified  shape  consistent  with  the  surface  ten- 
sion and  hydrostatic  pressure  requirements,  and  then  this  new  shape  wov.ld 
have  to  be  used  in  the  next  Iteration  of  the  aerodynamic  calculations,  and 
so  on.  Difficult  as  this  would  be,  one  might  be  ready  to  attempt  it  were 
It  not  for  the  fact  that  there  exists  no  general  method  for  calculating 
analytically  the  pressure  pattern  about  an  arbitrary  surface.  The  method 
of  superposition  of  a suitable  array  of  sources  and  sinks  which  Is  some- 
times useful  in  treating  the  flow  around  revolutes  le  only  sufficiently 
convergent  to  be  practicable  In  the  limit  of  very  elongated  revolutes  such 
as  dirigibles.  To  proceed  with  an  Iterative  method  In  which  each  aerodynamic 
calculation  had  to  be  performed  graphically  or  numerically  was  one  course 
open  to  the  writer;  but  he  has  Instead  chosen  to  proceed  at  thla  point  upon 
a different  tackj  namely,  to  analyze  a photograph  of  an  actual  drop  of 
known  size  and  falling  speed  In  order  to  deduce  these  aerodynamic  pressures 
which  would  be  so  difficult  to  calculate  directly. 

3.  Calculation  of  aerodynamic  surface  pressures 

On  the  'oasis  of  the  discussions  of  drop  morphology  given  in  the  preced- 
ing Section,  the  writer  adopts  the  following  hypothesis:  The  equilibrium 

shape  of  a large  drop,  bearing  at  most  a charge  small  compared  to  the  limit- 
ing v.nlue  Imposed  by  the  dielectric  strength  of  air  and  falling  at  terminal 
velocity.  Is  that  particular  shape  for  which  the  Joint  action  of  the  external 
aerodynamic  pressures  and  the  surface  pressure  Increments  Just  produce  an  in- 
ternal pressure  distribution  that  satisfies  the  hydrostatic  equation  within 
the  drop. 

This  hypothesis  has  led  the  writer  to  employ  the  following  method  for 
calculating  the  aerodynamically  developed  external  surface  pressures:  Given 

a photograph  of  a drop  of  known  size  and  falling  speed,  one  first  calculates 


-91- 


the  stagnation  pressure  developed  at  the  lover  pole  of  the  drop.  It  Is 
fortunately  one  of  the  veil  established  facts  of  fluid  dynamics  that  re- 
gardless of  almost  all  peculiarities  of  a given  flov  pattern  about  an  ob- 
ject immersed  in  a fluid  stream,  the  excess  pressure  developed  at  the 
leading  stagnation  point  is  v*\.  vhere  ^ is  the  fluid  density  and  v is 
the  speed  of  the  fluid  far  from  the  object,  measured  relative  to  that  object. 

Hence  this  first  step  Involves  no  approximations.  Next  one  measures,  on 
the  available  photograph,  the  radius  of  curvature  Rq  of  the  drop  surface 
profile  at  the  lover  stagnation  point.  For  reasons  of  axial  symmetry, 

Rl  « Rg  — Ro  at  this  point  of  the  drop,  i.e.,  the  drop  surface  is  locally 
a portion  of  a sphere  of  radius  Eg  et  the  stagnation  point.  Using  this 
measured  radius  in  (1)  to  compute  £> pB  at  the  stagnation  point,  and  adding 
the  result  to  the  computed  stagnation  pressure,  one  obtains  the  pressure 
prevailing  just  lnaide  the  drop  at  its  lover  pole.  Next,  using  the  hydro- 
static equation,  one  may  quickly  determine  the  Internal  pressure  p^(z)  at 
any  height  z measured  upwards  along  the  vertical  axis  of  symmetry  from  tero 
at  the  lover  pole.  Ilien,  as  long  as  internal  circulations  produce  only  neg- 
ligible internal  pressure  gradients,  the  pressure  just  inside  the  drop  sur- 
face at  height  z is  equal  to  that  already  determined  for  the  point  along  the 

axis  at  that  height.  Finally,  if  one  can  determine  from  the  drop  photograph 

\ 

the  values  R^(z)  ana  Rg(z),  then  (l)  can  be  used  to  compute  &Ps(r)>  an<*  ex- 
tracting this  from  Pi(z)  yields  pe(z),  the  external  aerodynamics lly  Induced 
pressure  at  height  z.* 

From  the  preceding  discussion  it  can  be  seen  that  the  success  of  this 

method  hinges  upon  being  able  to  determine  R^(z)  and  R2(z)  from  a single  slde- 

vlev  photograph  of  a given  drop.  The  technique  for  doing  this  turns  out  to 

%ote  that  all  pressures  represent  algebraic  excesses  over  the  prevailing 
barometric  pressure,  and  that  the  slight  variation  in  the  latter  through  the 
height  Interval  spanned  by  the  drop  at  any  instant  Is  lgnorable  because  it  is 
only  of  the  order  of  10-3  times  the  internal  hydrostatic  pressure  variation 
in  the  same  Interval. 
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p 

b«  quite  stra lghtf  orvard . 

One  of  the  tvo  principal  radii!  say  R^( z ) , la  simply  the  radius  of 
curvature  of  the  meridional  profile  at  the  height  i.  This  can  be  aeasured 
on  a tracing  made  fro*  an  enlargement  of  a photograph  of  the  vertical  cross 
section,  such  as  the  one  shovn  here  In  Figure  1.  One  constructs  normals  to 
the  profile  curve  at  each  of  a fairly  dense  series  of  points  spaced  regularly 
along  the  profile  (the  writer  used  a 4 50 -90°  prism  tangsntometer  In  the  deter- 
mination of  these  normals)  and  from  these  the  values  of  R^(z)  can  be  deter- 
mined by  measuring  the  distance  along  the  normal  at  z to  the  point  of  inter- 
section with  the  normal  drawn  from  the  next  adjoining  point  on  the  profile. 

In  Figure  4 the  distance  for  point  P Is  shovn  as  PCi.  By  using  an  en- 
largement factor  of  about  20  and  by  using  a fairly  dense  set  of  points  (fif- 
teen in  all)  along  the  profile,  it  was  found  to  be  easy  to  achieve  precision 
of  about  5 per  cent  in  R1(z),  as  ascertained  from  repeated  trials. 

Rp(z)  Is  even  more  easily  determined  sines  this  second  principal  radius, 

p 

for  a surface  of  revolution,  can  be  shovn  to  be  simply  the  distance  from 
the  profile  point  at  z to  the  axis  of  revolution  aeasured  along  the  local 
normal  to  the  profile  et  z.  In  Figure  4,  Rg  Is  shown  as  the  distance  PCg. 

The  normals  already  constructed  in  the  process  of  finding  R^(z)  facilitate 
rapid  determination  of  Rg  (z). 


p 

The  writer  Is  Indebted  to  Dr.  J.  M.  Keller  of  the  Department  of  Physics, 
Iova  Stats  Collage,  for  examining  and  solving  the  problem  In  differential 
geometry  that  underlies  the  technique  for  determining  the  principal  radii. 
Hla  solution  is  too  lengthy  to  be  reproduced  here  and  the  writer  knows  of 
no  published  eolutlon.  However,  subsequent  to  Dr.  Keller's  solving  this 
problem,  the  writer  came  across  a brief  statement  (Adam,  1949*  p.  366)  which 
agrees  with  hi?  result  and  Implies  that  the  same  problem  has  been  treated 
previously. 
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Ftgure  4.  Determination  of  R^(t)  and  Rg(r).  /or  point  P, 

*l(»)  - PC^  and  RgU)  ■ PCj>. 

Having  established  a method  for  determining  pe(z),  there  remains  only 
the  question  of  availability  of  suitable  photographs.  Since  It  Is  obvious 
that  satisfactory  photographic  resulte  depend  on  using  extremely  high  lllua- 
Inatlona  to  permit  exceedingly  short  exposure  times  (order  of  10"^  sec),  it 
la  not  entirely  surprising  that  the  vrlter  has  found  virtually  no  existing 
photographs  of  either  artificial  or  natural  drops  for  vhlch  the  falling  speed 
and  drop  site  are  accurately  determined.  In  fact,  the  vrlter  has  succeeded 
In  obtaining  only  a single  photograph  fulfilling  these  requirements.  This 
photograph,  reproducsd  here  as  Figure  1,  vas  taken  by  Dr.  Chojl  Magoao  of 
the  Hokkaido  Imperial  University,  under  conditions  vhere  both  terminal  vel- 
ocity and  drop  volume  could  be  measured.  All  of  the  remaining  discussion  Is 
necessarily  based  upon  this  one  photograph. 

Since  only  the  equivalent  spherical  diameter  of  the  drop  in  Figure  1 
vas  knovn,  the  exact  distance  seals  on  the  photograph  had  to  be  determined 
preliminarily.  On  a tracing  of  an  enlargement  of  the  photograph  the  cross- 
sectional  area  of  the  drop  vas  divided  Into  a large  number  of  horltontal 
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•trlpa  of  equal  vertical  width &z.  Next,  for  each  each  narrow  strip,  the 
radial  distance  rt  was  measured  from  a point  on  the  axis  of  revolution  (mid- 
way between  the  upper  and  lower  edges  of  the  strip)  out  to  toe  meridional 

o 

profile  curve.  Then  the  quantities  7^  r^  ^z  were  computed  for  each  of  the 
entire  series  of  strips  and  their  sum  equated  to  the  actual  drop  volume  as 
given  by  Magono.  From  this  equality,  the  scale-factor  of  the  enlargement 
was  ascertained  for  use  in  converting  all  subsequently  measured  distances 
on  the  photograph  to  true  distances. 

Proceeding  In  the  manner  now  fully  outlined,  the  writer  determined  p#(i ) 
for  tfcgono'e  drop.  The  results  are  presented  In  Figure  5«  To  aid  in  Identi- 
fying the  positions  along  the  meridional  profile  where  pe  assumes  certain 
values  of  particular  Interest,  five  points  (A  to  are  labeled  on  the  pres- 
sure curve  and  the  locations  of  these  points  are  Indicated  on  the  Inset 
sketch  of  the  drop  profile.  The  lower  stagnation  point  la  subjected  to  an 
external  pressure  of  460  dyne  cm-2  In  excess  of  barometric  pressure,  but  as 
the  air  accelerates  in  sweeping  up  auu  around  to*  drop,  the  surface  pressure 
falls,  reaching  a minimum  of  -590  dyne  cm-2  near  point  C,  Just  below  the 
point  of  maximum  horizontal  cross-sect Ion.  The  pressures  In  the  profile 
Interval  C-D  are  less  accurately  determined  than  those  for  other  loci  since 
here  the  profile's  radius  of  curvature  Is  changing  quite  rapidly  and  Is  rtry 
small.  By  performing  the  measurements  in  this  region  five  times  and  averag- 
ing the  results,  the  accepted  values  were  found  to  deviate  by  about  10  per 
cent  from  their  respective  extremal  values.  As  the  air  passes  C the  pressure 
very  suddenly  rises  to  a local  maximum  at  D and  then  falls  off  again  very 
slowly  towards  the  upper  pole  at  S,  with  a alight  (and  somewhat  questionable) 
rise  right  at  the  pole. 

If  the  airflow  around  this  drop  had  been  potential  flow,  the  stream- 
lines would  have  closed  In  above  the  drop  In  such  » vay  that  the  pressure 
cjrve  would  have  returned  to  positive  values  on  the  upper  surface,  ultimately 
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Figure  5.  Aerodynamic  preset* re  distribution  deduced  for  Magono'e 


drop  (Figure  1).  Letters  refer  to  Inset  sketch. 


9 (degrees)  > 

Figure  6.  Aerodynamic  pressure  distribution  around  a sphere  at 
Re^l.fixlO^  and  2,5*i(y  after  Jhge.  Dotted  curve  represents  profile 
for  the  case  of  potential  flov.  Abscissa  plotted  in  units  of  v^. 


-96- 


attaining  a Talus  of  k60  dyne  cm  at  an  upper  stagnation  point  at  E,  That 
the  pressures  are,  Instead,  found  to  remain  negatlre  ( l.e. , pressure  less 
than  barometric  pressure)  will  be  shown,  In  the  next  section,  to  Imply  that 
separation  of  the  laminar  boundary  layer  occurs  In  the  flow  around  large 
raindrops. 

1*.  Discussion  of  results 

To  review  the  logic  of  the  present  study,  It  must  be  noted  that  It  has 
not  been  found  possible  here  to  predict  the  aerodynamic  pressure  effects  on 
drop  shape  by  proceeding  directly  from  fundamental  principles  of  fluid  dynam- 
ics; Instead  the  aerodynamic  pressures  have  been  deduced  from  the  writer's 
hypothesis  of  drop  shape.  Therefore  It  follows  that  the  test  of  the  correct- 
ness of  the  shape  hypothesis  must  Involve  a comparison  between  the  results 
shown  In  Figure  5 and  any  available  and  pertinent  experimental  data. 

It  is  understandable  that  the  literature  of  experimental  aerodynamics 
does  not  contain  any  data  on  observed  pressure  patterns  around  objects  of 
raindrop  shape  since  this  Is  a quite  -nususl  shape.  Furthermore,  It  would 
not  be  easy  to  set  about  securing  precisely  the  data  desired  because  of  the 
following  considerations:  If  a "model"  of  a large  raindrop  were  to  be  fab- 

ricated with  an  array  of  static  pressure  orifices  distributed  over  Its  sur- 
face, the  necessary  internal  tubing  would  probably  preclude  use  of  a model 
smaller  than  four  or  five  centimeters  In  cross  section,  l.e.,  the  model 
would  have  to  be  six  or  seven  times  larger  than  the  prototype.  To  preserve 
dynamic  similarity  between  model  and  prototype  in  the  Reynolds  sense,  the  ex- 
perimental velocity  could  be  only  one -sixth  or  ons -seventh  of  the  prototype 
velocity  of  fall  (about  8 m sec”*)  if  air  were  used  as  the  model  medium;  but 
at  such  low  airspeeds,  the  uncertainties  in  measurement  of  the  slight  static 
pressure  changes  would  pose  rather  serious  obstacles  to  the  success  of  the 
experiment.  To  use  water  as  a medium  can  be  seen  to  provide  some  little 
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advantage,  since  the  kinematic  Tlacoalty  of  vater  le  only  about  one-tenth 
that  of  air  but  the  density  Is  about  a thousand  tinea  that  of  air,  so  for 
given  Reynolds  number,  the  pressure  effects  vould  be  ten  tines  greater  In 
vater  than  in  air.  If  any  such  model  studies  vere  to  be  undertaken,  one 
vould  probably  choose  to  use  air,  employ  a conveniently  large  drop  model, 
and  then  sacrifice  dynamic  similitude  In  order  to  get  at  least  a rough  no- 
tion of  the  behavior  of  the  flov.  The  vrlter  has  not  contemplated  doing 
this  himself,  but  the  results  vould  be  of  considerable  Interest. 

In  vlev  of  the  lack  of  precisely  the  type  cf  observational  data  needed 
fbr  comparison  vlth  the  present  results,  It  becomes  a next-best  substitute 
to  use  data  on  spheres . Even  for  spheres  one  finds  gap  In  the  experimental 
results  for  Reynolds  numbers  lying  betveen  the  upper  limit  of  the  Stokes  lav 
range  and  the  lover  limit  of  the  region  of  critical  Reynolds  numbers  for 
transition  to  a turbulent  boundary  layer.  Drag  data  vere  found  to  be  abund- 
ant for  the  entire  range,  but  not  surface  pressure  data.  Thy*  It  finally 
became  Decease ry  merely  to  look  at  pressure  profiles  for  spheres  at  Reynolds 
numbers  of  the  order  of  10“  (Fage,  1937)*  despite  the  fact  that  Magono's  drop 
fell  vlth  a Reynolds  number  of  only  about  4xl(>5. 

Two  of  Page's  curves  for  the  profile  around  a sphere  (for  Re£2.5xlC)5 
and  Res}1.6xl0^)  are  plotted  here  In  Figure  6.  The  similarities  betveen 
these  curves  and  that  of  Figure  5 are  sufficient  to  provide  considerable  as- 
surance that  the  present  drop-shape  hypothesis  is  at  least  fairly  close  to 
the  truth.  The  pressures  around  the  dovnstream  hemisphere  of  Fage's  sphere, 
like  those  deduced  here  for  the  upper  surface  of  Magono's  drop,  fall  to  re- 
turn to  positive  values  (as  do  the  pressures  in  potential  flov,  dotted  curve 
of  Figure  6)  after  starting  to  rise  fairly  rapidly  Just  ahead  of  the  region 
of  maximum  cross  section.  Now  J.t  Is  veil  knovn  that  this  feature,  In  curves 
such  as  those  of  Figure  6;  Is  due  to  the  occurrence  of  separation  In  the 
boundary  layer.  At  Reynolds  numbers  be lev  those  for  vhlch  the  boundary  layer 
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becomes  turbulent  (Rsc*4xlo5)  but  above  the  Stokes  la*  rang*  (As&l),  the 
streamlines  full  to  oloaa  In  downstream  from  the  apbera  and  by  thla  failure 
to  create  a downstream  stagnation  point  the  preaeuraa  remain  low  downstream 
from  the  zone  of  separation.  Separation  may  be  Inhibited  by  fairing  out  the 
downstream  portion  of  a body  Into  a gently  tapering  form,  but  in  bod  lea  of 
large  curvature  aft  (e.g.,  a deformed  raindrop)  separation  Is  known  to  be 
easily  established.  The  fact  that  the  pressures  deduced  hare  for  Magono's 
drop  do  exhibit  this  behavior  nay  be  regarded,  then,  as  constituting  fairly 
strong  evidence  In  favor  of  the  shape  hypothesis  used  to  determine  these 
pressures.  In  still  further  support  of  this  conclusion  there  exist  two  ex* 
perlmental  observations  on  actual  water  drops  which  will  next  be  shown  to 
be  quite  consistent  with  the  present  deduction  of  separation  In  the  airflow 
around  large  raindrops. 

First,  Gunn  (191*9)  has  reported  a curious  tendency  for  drops  of  one 
certain  size  (about  0.5  *0  radius)  to  undergo  marked  sideslipping  as  they 
fall.  Gunn  has  shown  very  convincingly  that  this  must  be  due  to  a resonance 
phenomenon  Involving  the  natural  frequency  of  mechanical  oscillation  of  the 
drops  and  the  frequency  with  which  eddies  are  shed  from  the  upper  surface  of 
tha  falling  drop.  Gunn  Invoked  Moeller's  (1938)  extensive  results  on  eddy 
frequencies  for  spheres  to  show  that  a sphere  of  0.5  mm  radius  sheds  eddies 
at  a frequency  of  about  500  cps,  almost  exactly  the  natural  frequency  of 
ellipsoidal  vlbratloa  of  a water  sphere  of  that  size.  Gunn's  detection  of 
thla  eddying  phenomenon  constitutes  clear  evidence  of  separation,  since 
eddies  can  only  be  shed  from  a "deadwater"  region  bounded  by  a separating 
stream  surface  extending  downstream  from  the  given  object.  Furthermore, 
Gunn's  observations  concerned  a drop-size  much  smaller  than  the  one  here 
analyzed;  so  at  the  much  higher  Reynolds  number  at  which  Magono'a  drop  fell, 
separation  should  almost  certainly  be  expected  (and  this  the  more  ao  because 
of  the  sharper  curvature  on  the  upper  surface  of  the  larger  and  more  flat- 


tened  drop). 

Second,  Blanchard  (1950)  reported  that  when  ora  drop  of  appropriate 
site  Is  inserted  into  the  alrstream  above  another  drop  which  la  aerodynam- 
lcally  suspended  In  a vertical  slrstresm,  the  upper  descends  upon  the  lower 
along  a peculiar  spiral  path.  Wlnny  (1952)  has,  by  means  of  photographs  of 
the  flow  behind  spheres  at  which  separation  was  occurring.  Illustrated  what 
appeared  to  be  a spiral  eddy  pattern  In  the  Interval  Re£2xlo3  to  6x10^,  so 
Blanchard's  observations  may  be  taken  as  further  evidence  for  separation  In 
the  flow  around  water  drops.  A very  complete  discussion  of  the  sake  phenom- 
ena behind  solid  spheres  has  been  given  by  Moeller  (1936),  vbo  concluded  that 
separation  first  appears  at  about  As  = 150  and  that  periodic  eddy  detachment 
begins  at  about  Re  = 4^0  for  solid  spheres.  For  a somewhat  deformed  drops, 
these  phenomena  would  be  expected  to  appear  at  somewhat  lower  Reynolds  num- 
bers. 

In  all,  there  seems  to  be  very  good  reason  for  believing  that  separation 
occurs  In  the  airflow  around  all  raindrops  with  diameters  greater  than  about 
0,5  - 1,0  mm.  Since  the  present  shape  hypothesis  has  led  to  the  deduction  of 
a pressure  profile  of  a type  entirely  different  from  that  obtained  with  po- 
tential flow  (see  dotted  curve,  Figure  6)  but  quite  similar  to  that  character- 
istic of  viscous  flow  at  high  Reynolds  number  (separating  boundary  layer), 
there  would  appear  to  be  sufficient  ground  for  accepting  the  shape  hypothesis 
adopted  here . 

5*  Implications  of  separation 

It  should  be  emphaslted  that  the  calculation  of  pe(i)  for  Magono'a  drop 
has  Interest  not  only  In  that  It  provides  a check  on  the  present  theory  of 
drop  shape,  but  also  In  that  It  focuses  attention  on  the  phenomenon  of  separ- 
ation Itself.  Gunn's  (194$)  observations  Imply  that  this  phenomenon  Is  al- 
ready veil  established  In  flow  about  raindrops  of  acXidl  else  (1  era  dlaneler) 


-100- 


falling  at  Reynolds  numbers  of  laaa  than  300,  which  finding  la  In  good  agree - 
nt  with  Moeller' e work  on  solid  spheres.  Bence  separation  appears  to  be  a 
phenomenon  to  be  reckoned  with  In  any  theory  of  microphyslcal  processes  In- 
volving raindrops.  The  kinematics  of  Ion  deposition  on  falling  raindrops, 
for  example,  must  be  admitted  to  be  affected  to  some  extent  by  the  failure 
for  simple  potential  flcv  to  occur  on  the  upper  surfaces  of  drope.  Anal- 
yses of  heat  and  vater  vapor  transfer  to  or  from  vater  drops  must  take  sep- 
aration Into  account.  As  Blanchard  has  found  experimentally,  the  dynamlos 
of  a "fall-on"  collision  between  two  drops  will  be  sensibly  affected  by  the 
existence  of  a pulsatory  wake  set  up  behind  a drop  about  which  the  flow' Is 
undergoing  separation,  finally,  the  else -dependence  of  the  terminal  velocity 
of  raindrops  must  be  largely  due  to  the  well  established  effect  of  separation 
on  form  drag.  For  Reynolds  numbers  in  the  approximate  Interval  101  to  Wp, 
the  drag  coefficient  of  a sphere  Is  nearly  constant  because  In  this  range 
form  drag  (due  to  the  presence  of  a deadwater  region)  rather  than  skin  fric- 
tion drag  is  of  controlling  Importance,  and  separation  produces  a deadwater 
region  whose  dlmenslc  « remain  roughly  Independent  of  Reynolds  number  until 
the  latter  reaches  the  critical  transition  value  (meteorologically  (unat- 
tainable for  raindrops,  which  break  up  well  below  Rec)  • A*  wOfiw  vbcr  y the  ef- 
fects due  to  separation,  whose  existence  In  raindrop  aerodynamics  seems  quite 
clearly  Indicated,  must  constitute  a significant  factor  In  the  physics  of  rain. 

The  detection  of  separation  sheds  further  Interesting  light  on  the  rea- 
sons for  the  probable  non-existence  of  appreciable  Internal  circulation  In- 
side raindrops.  The  skln-frlctlon  drag  at  the  surface  of  a liquid  drop  about 
which  the  flow  undergoes  uaparatlon  Is  markedly  less  favorable  to  the  estab- 
lishment of  Internal  circulation  of  the  spherical  vortex  type  than  Is  the 
skln-frlctlon  drag  in  the  non-separating  csss.  Starting  from  a zero  value 
at  the  lower  stagnation  point,  the  downstream  drag  of  air  on  vater  must 
reach  Its  maximum  value  well  forward  of  the  separation  point  and  then  must 
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fall  once  more  to  eero  at  the  separation  point  where,  by  definition,  there 
exists  a zero  radial  Telocity  shear  right  at  the  drop  surface  (aee  Goldstein, 
193d,  and  the  suggested  flow  pattern  of  Figure  7).  Then,  downstream  from 
the  separation  point, 


Figure  7.  Suggested  flow  pattern  around  a large 
raindrop  when  separation  occurs. 

there  will  tend  to  exist  a ring  vortex  of  reversed  circulation  concentric 
with  the  symmetry  axle  of  the  drop.  The  presence  of  such  a circulation  over 
the  upper  portion  of  the  drop  surface  will  serve  to  cut  down  any  Incipient 
Internal  vortical  circulations  that  do  tend  to  become  established  by  the  drag 
forces  acting  over  the  lower  surface  of  the  drop.  Hence,  once  one  has  recog - 
nlzed  that  separation  Is  characteristic  of  the  flow  arc.od  raindrops,  he  has 


available  a vary  strong  basis  for  predicting  that  nothing  like  a spherical 
vortex  of  the  Bond  type  will  occur  vithln  large  drops,  for  there  simply  does 
not  exist  a sustained  downstream  tugging  of  surface  air  on  eurfacevater  over 
the  entire  surface  of  the  drop.  This  argument  could  not,  of  course,  logically 
be  used  earlier,  since  to  assume  separation  from  the  start  would  be  to  poetu- 
late  the  very  phenomenon  whose  theoretical  deduction  has  here  constituted  the 
best  available  test  of  the  working  hypothesis  for  drop  shape.  Furthermore, 
at  the  time  that  the  writer  waa  still  striving  to  settle  the  question  of  the 
existence  or  non-existence  of  internal  circulations  he  was  not  yet  aware  of 
the  experimental  Indications  that  separation  occurs  In  the  flow  around  drops. 
In  the  analysis  of  Internal  circulations  carried  out  above  (Section  2)  a 
surface  circulation  speed  of  about  30  cm  aec~^  was  predicted.  That  this  re- 
sult Is  substantially  too  large  to  agree  with  the  observations  of  Blanchard 
and  of  K Inzer  la  now  understandable  In  retrospect,  for  In  that  analysis  it 
was  tacitly  assumed  that  the  surface  air  flow  Is  favorable  to  generation  of 
a vortical  circulation  everywhere  over  the  surface  of  a drop.  In  reality, 
the  opposing  Influences  of  frictional  drag  at  lower  and  upper  surfaces  of  a 
drop  must  be  expected  to  prevent  such  strong  circulation  from  developing. 

Finally,  it  is  to  be  noted  that  separation  appears  to  be  responsible  for 
the  asymmetry  of  a large  raindrop  with  respect  to  a horizontal  plane  torough 
its  center.  It  was  pointed  out  in  Section  2 above  that  If  only  surface  ten- 
sion and  hydrostatic  pressure  effects  controlled  drop  shape,  these  would  pro- 
duce drops  with  rounded  bottoms  and  rather  flattened  tops.  However,  because 
a stagnation  point  inevitably  occurs  at  the  lower  pole,  while  one  cannot  de- 
velop at  the  upper  pole  because  of  separation  of  the  boundary  layer,  it  fol- 
lows that  the  lower  aerodynamic  pressures  over  the  upper  surface  demand  an 
appreciable  larger  there  than  on  the  underside,  i.e.,  the  curvature  of 

the  upper  surface  must  bo  greater  in  order  to  satisfy  the  requirement  of  in- 
ternal hydrostatic  equilibrium  in  face  of  the  considerable  aerodynamic 
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pressure  difference  on  upper  and  lower  surfaces.  Hence  the  drop  becomes 
rounded  on  the  top  but  flattened  on  the  bottom,  as  revealed  by  photographs. 
This  asymmetry  of  large  drops , rather  veil  accounted  for  on  the  present  drop- 
shape  hypothesis!  vas  not  explained  by  either  Lenard  (1904)  nor  Spilhaue 
(1948)  In  their  studies  of  drop  morphology. 

6.  Summary 

A rev lev  of  previous  efforts  to  understand  the  peculiar  deformation 
characteristic  of  large  raindrops  revealed  that  little  has  been  done  In  the 
past  to  clarify  this  problem.  The  factors  of  surface  tension,  hydrostatic 
pressure  gradients,  external  aerodyuamlc  pressures,  electrostatic  charge, 
and  internal  circulations  vere  examined  quantitatively  and  only  the  first 
three  of  these  vere  found  to  be  significant  in  controlling  drop  shape. 
Adopting  the  hypothesis  that  the  equilibrium  shape  of  a raindrop  falling  at 
terminal  velocity  la  that  for  vhlch  the  aerodynamic  pressures  and  the  surface 
tension  pressure  Increments  conspire  to  Just  produce  an  Internal  pressure 
pattern  satisfying  the  hydrostatic  equation,  it  vas  found  possible,  oy  anal- 
yzing a single  drop  photograph  to  deduce  the  aerodynamic  pressure  profile 
along  a meridian  of  the  drop.  This  profile  revealed  clear  evidence  of  separ- 
ation of  the  laminar  boundary  layer  in  the  airflcv  around  the  drop.  Since 
separation  effects  could  also  be  shown  to  exist  in  certain  experimental  ob- 
sarvations,  it  has  been  concluded  that  the  vorking  hypothesis  for  drop  chape 
has  been  successful  In  explaining  the  long -re cognized  deformation  of  large 
raindrops.  A number  of  implication*  of  separation  vere  pointed  cut  qualita- 
tively, and  it  vas  noted  that  the  distribution  of  skin-friction  dreg  over  the 
surface  of  •»  raindrop  at  vhlch  separation  occurs  is  decidedly  unfavorable  to 
the  generation  of  Internal  vortical  circulation. 

7.  Suggestions  for  future  research 

That  the  writer  has  only  succeeded  in  securing  a single  drop  photograph 
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sult&ble  for  the  present  type  of  analysis  is  regrettable.  In  the  past,  pre- 
sumably, there  has  been  little  occasion  for  obtaining  such  difficult  photo- 
graphs under  conditions  where  drop  sire  and  felling  speed  could  be  deter- 
mined accurately.  However,  It  new  seems  very  desirable  for  the  experimental- 
ist to  obtain  a set  of  such  photographs  over  the  whole  urop-site  range  of 
meteorological  Interest  in  order  that  a complete  picture  of  the  flow  regions 
characteristic  of  the  various  drop  sizes  may  be  delineated  by  the  method 
developed  acre.  Such  information  should  prove  very  useful  In  extending  pres- 
ent theories  of  heat  and  vapor  transfer  to  and  from  raindrops. 

Secondly,  since  the  present  study  is  regarded  by  the  writer  as  only  a 
first  step  towards  the  more  significant  objective  of  understanding  the  breakup 
of  large  raindrops,  he  wishes  to  take  this  opportunity  to  urge  extension  of 
the  sort  of  laboratory  studies  so  well  begun  by  Blanchard  (1950).  The  wealth 
of  intriguing  questions  raised  by  Blanchard’s  work  and  the  light  that  could 
potentially  be  shed  on  the  dynamics  of  drop  breakup  by  further  studies  of 
this  sort  make  a continuation  of  the  research  most  desirable.  In  this  con- 
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work  by  McDoru  Id  (195l)>  where  a number  of  suggestions  for  improving  this  type 
of  study  have  been  offered. 

Finally,  it  would  still  be  desirable  to  have  more  conclusive  experimental 
data  on  the  intensities  of  Internal  circulation  in  water  drops  of  various 
sizes.  For  drops  Just  belcw  the  size  for  which  separation  first  appears 
(probably  around  0.5  mm  diameter)  it  may  be  possible  that  circulation  is  act- 
ually better  developed  than  for  larger  drops  cn  whose  upper  surface  a re  ■ 
versed  flow  develops. 
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